Постройте mesh УЧП
pdemesh(model)pdemesh(mesh)pdemesh(nodes,elements)pdemesh(model,u)pdemesh(___,Name,Value)pdemesh(p,e,t)pdemesh(p,e,t,u)h = pdemesh(___)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-образной мембраной. Используйте граничные условия Дирихле на , и постройте результат.
u = assempde('lshapeb',p,e,t,1,0,1);
pdemesh(p,e,t,u)
model — Объект моделиPDEModel | объект ThermalModel | объект StructuralModelОбъект модели, заданный как объект PDEModel, объект ThermalModel или объект StructuralModel.
Пример: model = createpde(1)
Пример: thermalmodel = createpde('thermal','steadystate')
Пример: structuralmodel = createpde('structural','static-solid')
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 должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.
pdemesh(model,'NodeLabels','on')