createPDEResults

Создайте объект решения

Эта страница описывает устаревший рабочий процесс. Новые возможности не могут быть совместимы с устаревшим рабочим процессом. Для соответствующего шага в рекомендуемом рабочем процессе смотрите solvepde и solvepdeeig.

Исходная (R2015b) версия createPDEResults имел только один синтаксис и создал PDEResults объект. Начиная с R2016a, вы обычно не должны использовать createPDEResults, потому что solvepde и solvepdeeig функции возвращают объекты решения. Кроме того, createPDEResults возвращает объект более нового типа, чем PDEResults. Если вы открываете существующий PDEResults объект, это преобразовано в StationaryResults объект.

Если вы используете один из более старых решателей такой как adaptmesh, затем можно использовать createPDEResults получить объект решения. Стационарные и зависящие от времени объекты решения имеют градиенты в наличии, тогда как PDEResults не включал градиенты.

Синтаксис

results = createPDEResults(model,u)

results =
createPDEResults(model,u,'stationary')

results = createPDEResults(model,u,utimes,'time-dependent')

results = createPDEResults(model,eigenvectors,eigenvalues,'eigen')

Описание

пример

results = createPDEResults(model,u) создает StationaryResults объект решения от model и его решение u.

Этот синтаксис эквивалентен results = createPDEResults(model,u,'stationary').

пример

results = createPDEResults(model,u,utimes,'time-dependent') создает TimeDependentResults объект решения от model, его решение u, и времена utimes.

пример

results = createPDEResults(model,eigenvectors,eigenvalues,'eigen') создает EigenResults объект решения от model, его решение для собственного вектора eigenvectors, и его собственные значения eigenvalues.

Примеры

свернуть все

Результаты эллиптической проблемы

Скрипт Open Live Script

Создайте StationaryResults объект от решения до эллиптической системы.

Создайте модель PDE для системы трех уравнений. Импортируйте геометрию скобки и постройте метки поверхности.

model = createpde(3);
importGeometry(model,'BracketWithHole.stl');
figure
pdegplot(model,'FaceLabels','on')
view(30,30)
title('Bracket with Face Labels')

figure
pdegplot(model,'FaceLabels','on')
view(-134,-32)
title('Bracket with Face Labels, Rear View')

Условия границы множества: столкнитесь 3, неподвижно, и в отрицательном z существует сила направление на поверхности 6.

applyBoundaryCondition(model,'dirichlet','face',4,'u',[0,0,0]);
applyBoundaryCondition(model,'neumann','face',8,'g',[0,0,-1e4]);

Установите коэффициенты, которые представляют уравнения линейной эластичности.

E = 200e9;
nu = 0.3;
c = elasticityC3D(E,nu);
a = 0;
f = [0;0;0];

Создайте mesh и решите задачу.

generateMesh(model,'Hmax',1e-2);
u = assempde(model,c,a,f);

Создайте StationaryResults объект из решения.

results = createPDEResults(model,u)

results = 
  StationaryResults with properties:

    NodalSolution: [14002x3 double]
       XGradients: [14002x3 double]
       YGradients: [14002x3 double]
       ZGradients: [14002x3 double]
             Mesh: [1x1 FEMesh]

Постройте решение для z- компонент, который является компонента 3.

pdeplot3D(model,'ColorMapData',results.NodalSolution(:,3))

Результаты зависящей от времени проблемы

Скрипт Open Live Script

Получите решение из параболической проблемы.

Проблемный тепловой поток моделей в теле.

model = createpde();
importGeometry(model,'Tetrahedron.stl');
pdegplot(model,'FaceLabels','on','FaceAlpha',0.5)
view(45,45)

Установите температуру на поверхности 2 - 100. Оставьте другие граничные условия в их значениях по умолчанию (изоляция).

applyBoundaryCondition(model,'dirichlet','face',2,'u',100);

Установите коэффициенты моделировать параболическую проблему с 0 начальными температурами.

d = 1;
c = 1;
a = 0;
f = 0;
u0 = 0;

Создайте mesh и решите УЧП в течение многих времен от 0 до 200 с шагом 10.

tlist = 0:10:200;
generateMesh(model);
u = parabolic(u0,tlist,model,c,a,f,d);

168 successful steps
0 failed attempts
324 function evaluations
1 partial derivatives
28 LU decompositions
323 solutions of linear systems

Создайте TimeDependentResults объект из решения.

results = createPDEResults(model,u,tlist,'time-dependent');

Постройте решение на поверхности геометрии во время 100.

pdeplot3D(model,'ColorMapData',results.NodalSolution(:,11))

Результаты задачи о собственных значениях

Скрипт Open Live Script

Создайте EigenResults объект от решения до задачи о собственных значениях.

Создайте геометрию и mesh для L-образной мембраны. Примените граничные условия Дирихле ко всем ребрам.

model = createpde;
geometryFromEdges(model,@lshapeg);
generateMesh(model,'Hmax',0.05,'GeometricOrder','linear');
applyBoundaryCondition(model,'dirichlet','Edge',1:model.Geometry.NumEdges,'u',0);

Решите задачу о собственных значениях для коэффициентов c = 1, a = 0, и d = 1. Получите решения для собственных значений от 0 до 100.

c = 1;
a = 0;
d = 1;
r = [0,100];
[eigenvectors,eigenvalues] = pdeeig(model,c,a,d,r);

              Basis= 10,  Time=   0.02,  New conv eig=  0
              Basis= 14,  Time=   0.04,  New conv eig=  0
              Basis= 18,  Time=   0.04,  New conv eig=  1
              Basis= 22,  Time=   0.06,  New conv eig=  2
              Basis= 26,  Time=   0.08,  New conv eig=  3
              Basis= 30,  Time=   0.09,  New conv eig=  5
              Basis= 34,  Time=   0.12,  New conv eig=  5
              Basis= 38,  Time=   0.15,  New conv eig=  7
              Basis= 42,  Time=   0.18,  New conv eig=  8
              Basis= 46,  Time=   0.22,  New conv eig= 11
              Basis= 50,  Time=   0.25,  New conv eig= 12
              Basis= 54,  Time=   0.33,  New conv eig= 14
              Basis= 58,  Time=   0.36,  New conv eig= 14
              Basis= 62,  Time=   0.41,  New conv eig= 16
              Basis= 66,  Time=   0.46,  New conv eig= 18
End of sweep: Basis= 66,  Time=   0.46,  New conv eig= 17
              Basis= 27,  Time=   0.51,  New conv eig=  0
              Basis= 31,  Time=   0.53,  New conv eig=  0
              Basis= 35,  Time=   0.55,  New conv eig=  0
End of sweep: Basis= 35,  Time=   0.57,  New conv eig=  0

Создайте EigenResults объект из решения.

results = createPDEResults(model,eigenvectors,eigenvalues,'eigen')

results = 
  EigenResults with properties:

    Eigenvectors: [1440x17 double]
     Eigenvalues: [17x1 double]
            Mesh: [1x1 FEMesh]

Постройте решение для режима 10.

pdeplot(model,'XYData',results.Eigenvectors(:,10))

Входные параметры

свернуть все

`model` — Модель PDE
`PDEModel` объект

Модель PDE в виде PDEModel объект.

Пример: model = createpde

`u` — Решение для УЧП
вектор | матрица

Решение для УЧП в виде вектора или матрицы.

Пример: u = assempde(model,c,a,f);

`utimes` — Времена для решения для УЧП
монотонный вектор

Времена для решения для УЧП в виде монотонного вектора. Эти времена должны совпасть с tlist времена, которые вы задали для решения hyperbolic или parabolic решатели.

Пример: utimes = 0:0.2:5;

`eigenvectors` — Решение для собственного вектора
матрица

Решение для собственного вектора в виде матрицы. Предположим

Np количество узлов mesh
N количество уравнений
ev количество собственных значений, заданных в eigenvalues

Затем eigenvectors имеет размер Np- N- ev. Каждый столбец eigenvectors соответствует собственным векторам одного собственного значения. В каждом столбце, первом Np элементы соответствуют собственному вектору уравнения 1, оцененного в узлах mesh, следующем Np элементы соответствуют уравнению 2 и так далее.

`eigenvalues` — Решение для собственного значения
вектор

Решение для собственного значения в виде вектора.

Выходные аргументы

свернуть все

`results` — Решение для УЧП
`StationaryResults` объект (значение по умолчанию) | `TimeDependentResults` возразите | `EigenResults` объект

Решение для УЧП в виде StationaryResults объект, TimeDependentResults объект или EigenResults объект. Создайте results использование solvepde, solvepdeeig, или createPDEResults.

Пример: results = solvepde(model)

Советы

Размерности возвращенных решений и градиентов совпадают с возвращенными solvepde и solvepdeeig. Для получения дополнительной информации смотрите Размерности Решений, Градиентов и Потоков.

Алгоритмы

Процедура для оценки градиентов в узловых местоположениях следующие:

Вычислите градиенты в точках Гаусса, расположенных в каждом элементе.
Экстраполируйте градиенты в узловых местоположениях.
Насчитайте значение градиента от всех элементов, которые встречаются в узле. Этот шаг необходим из-за межэлементного разрыва градиентов. Элементы, которые соединяются в том же узле, дают различные экстраполируемые значения градиента для точки. createPDEResults выполняет взвешенное областью усреднение для 2D сеток и взвешенное объемом усреднение для 3-D сеток.

Документация

createPDEResults

Синтаксис

Описание

Примеры

Результаты эллиптической проблемы

Результаты зависящей от времени проблемы

Результаты задачи о собственных значениях

Входные параметры

`model` — Модель PDE
`PDEModel` объект

`u` — Решение для УЧП
вектор | матрица

`utimes` — Времена для решения для УЧП
монотонный вектор

`eigenvectors` — Решение для собственного вектора
матрица

`eigenvalues` — Решение для собственного значения
вектор

Выходные аргументы

`results` — Решение для УЧП
`StationaryResults` объект (значение по умолчанию) | `TimeDependentResults` возразите | `EigenResults` объект

Советы

Алгоритмы

Смотрите также

Темы

Документация Partial Differential Equation Toolbox

Поддержка

Документация

createPDEResults

Синтаксис

Описание

Примеры

Результаты эллиптической проблемы

Результаты зависящей от времени проблемы

Результаты задачи о собственных значениях

Входные параметры

model — Модель PDE PDEModel объект

u — Решение для УЧП вектор | матрица

utimes — Времена для решения для УЧП монотонный вектор

eigenvectors — Решение для собственного вектора матрица

eigenvalues — Решение для собственного значения вектор

Выходные аргументы

results — Решение для УЧП StationaryResults объект (значение по умолчанию) | TimeDependentResults возразите | EigenResults объект

Советы

Алгоритмы

Смотрите также

Темы

Документация Partial Differential Equation Toolbox

Поддержка

`model` — Модель PDE
`PDEModel` объект

`u` — Решение для УЧП
вектор | матрица

`utimes` — Времена для решения для УЧП
монотонный вектор

`eigenvectors` — Решение для собственного вектора
матрица

`eigenvalues` — Решение для собственного значения
вектор

`results` — Решение для УЧП
`StationaryResults` объект (значение по умолчанию) | `TimeDependentResults` возразите | `EigenResults` объект