Постройте mesh УЧП
pdemesh(___, строит mesh или данные о решении с помощью любого из аргументов в предыдущих синтаксисах и одном или нескольких Name,Value)Name,Value парные аргументы.
pdemesh( узел УЧП графиков или треугольные данные p,e,t,u)u использование сетчатого графика. Графики функций данные об узле, если u вектор-столбец и треугольные данные если u вектор-строка.
Если вы хотите иметь больше контроля над своим сетчатым графиком, использовать pdeplot или pdeplot3D вместо pdemesh.
h = pdemesh(___)
Создайте сетчатый график и отобразите узел и метки элемента mesh.
Создайте модель PDE. Включайте геометрию встроенной функции lshapeg. Поймайте в сети геометрию.
model = createpde; geometryFromEdges(model,@lshapeg); mesh = generateMesh(model);
Постройте mesh.
pdemesh(model)
В качестве альтернативы можно построить mesh при помощи mesh как входной параметр.
pdemesh(mesh)
Другой подход должен использовать узлы и элементы mesh как входные параметры для pdemesh.
pdemesh(mesh.Nodes,mesh.Elements)
Отобразите метки узла.
pdemesh(model,'NodeLabels','on')
Используйте xlim и ylim увеличить масштаб конкретных узлов.
xlim([-0.4,0.4]) ylim([-0.4,0.4])
Метки элемента отображения.
pdemesh(model,'ElementLabels','on') xlim([-0.4,0.4]) ylim([-0.4,0.4])
Примените граничные условия, задайте коэффициенты и решите УЧП.
applyBoundaryCondition(model,'dirichlet', ... 'Edge',1:model.Geometry.NumEdges, ... 'u',0); specifyCoefficients(model,'m',0,... 'd',0,... 'c',1,... 'a',0,... 'f',1); generateMesh(model); results = solvepde(model)
results =
StationaryResults with properties:
NodalSolution: [1177x1 double]
XGradients: [1177x1 double]
YGradients: [1177x1 double]
ZGradients: []
Mesh: [1x1 FEMesh]
u = results.NodalSolution;
Постройте решение в узловых местоположениях при помощи pdemesh.
pdemesh(model,u)
pdemesh функция игнорирует NodeLabels и ElementLabels когда вы отображаете данные о решении на графике как 3-D график.
Создайте модель PDE, включайте геометрию и поймайте в сети ее.
model = createpde; importGeometry(model,'Plate10x10x1.stl'); generateMesh(model,'Hmax',5);
Постройте mesh, устанавливающую прозрачность на 0,5.
pdemesh(model,'FaceAlpha',0.5)
Найдите элементы сопоставленными с геометрической областью.
Создайте модель PDE.
model = createpde;
Включайте геометрию встроенной функции lshapeg. Постройте геометрию.
geometryFromEdges(model,@lshapeg); pdegplot(model,'FaceLabels','on','EdgeLabels','on')
Сгенерируйте mesh.
mesh = generateMesh(model,'Hmax',0.5);Найдите элементы сопоставленными с поверхностью 2.
Ef2 = findElements(mesh,'region','Face',2);
Подсветите эти элементы зеленого цвета на сетчатом графике.
figure pdemesh(mesh,'ElementLabels','on') hold on pdemesh(mesh.Nodes,mesh.Elements(:,Ef2),'EdgeColor','green')
Постройте mesh для геометрии L-образной мембраны.
[p,e,t] = initmesh('lshapeg'); [p,e,t] = refinemesh('lshapeg',p,e,t); pdemesh(p,e,t)
Теперь решите уравнение Пуассона по геометрии, заданной L-образной мембраной. Используйте граничные условия Дирихле on , и постройте результат.
u = assempde('lshapeb',p,e,t,1,0,1);
pdemesh(p,e,t,u)
model — Объект моделиPDEModel возразите | ThermalModel возразите | StructuralModel возразите | ElectromagneticModel объектОбъект модели в виде PDEModel объект, ThermalModel объект, StructuralModel объект или ElectromagneticModel объект.
Пример: model = createpde(3)
Пример: thermalmodel = createpde('thermal','steadystate')
Пример: structuralmodel = createpde('structural','static-solid')
Пример: emagmodel = createpde('electromagnetic','electrostatic')
u — Решение для УЧПРешение для УЧП в виде вектора или матрицы.
Пример: results = solvepde(model); u = results.NodalSolution; или u = assempde(model,c,a,f);
mesh — Объект ячейкиMesh свойство PDEModel возразите | выход generateMeshОбъект ячейки в виде Mesh свойство PDEModel возразите или как выход generateMesh.
Пример: model.Mesh
nodes — Узловые координатыУзловые координаты в виде 2 NumNodes матрицей для 2D mesh и 3 NumNodes матрицей для 3-D mesh. NumNodes является количеством узлов.
elements — Матрица смежности элемента в терминах идентификаторов узлаМатрица смежности элемента в терминах идентификаторов узла в виде NodesPerElem-by-NumElements матрица. NodesPerElem является количеством узлов на элемент. Линейные сетки содержат только угловые узлы, таким образом, существует три узла на 2D элемент и четыре узла на 3-D элемент. Квадратичные сетки содержат угловые узлы и узлы посреди каждого ребра элемента. Для квадратичных сеток существует шесть узлов на 2D элемент и 10 узлов на 3-D элемент.
Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.
pdemesh(model,'NodeLabels','on')