Постройте mesh УЧП
pdemesh(___, строит график mesh или решения с помощью любого из аргументов в предыдущих синтаксисах и одного или нескольких Name,Value)Name,Value аргументы в виде пар.
pdemesh( строит графики для данных о узле PDE или треугольнике p,e,t,u)u использование сетчатого графика. А графики функций данные узла, если u - вектор-столбец, и данные треугольника, если u является вектор-строка.
Если вы хотите иметь больше контроля над своим сетчатым графиком, используйте pdeplot или pdeplot3D вместо pdemesh.
h = pdemesh(___)
Создайте сетчатый график и отобразите метки узла и элемента mesh.
Создайте модель УЧП. Включите геометрию встроенной функции 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)
The pdemesh функция игнорирует NodeLabels и ElementLabels когда вы строите график данных решения как 3-D график.
Создайте модель УЧП, включите геометрию и создайте сетку.
model = createpde; importGeometry(model,'Plate10x10x1.stl'); generateMesh(model,'Hmax',5);
Постройте график mesh, задав прозрачность 0,5.
pdemesh(model,'FaceAlpha',0.5)
Найти элементы, связанные с геометрической областью.
Создайте модель УЧП.
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 | объекта 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 - Объект MeshMesh свойство PDEModel объект | выход generateMeshОбъект Mesh, заданный как Mesh свойство PDEModel объект или как выход generateMesh.
Пример: model.Mesh
nodes - Узловые координатыУзловые координаты, заданные как 2-байтовая NumNodes матрица для 2-D mesh и 3-байтовая NumNodes матрица для 3-D mesh. NumNodes является числом узлов.
elements - Матрица связности элемента в терминах идентификаторов узлаМатрица связности элемента в терминах идентификаторов узла, заданная как матрица NodesPerElem -by- NumElements. NodesPerElem является числом узлов по элементам. Линейные сетки содержат только угловые узлы, поэтому существует три узла на 2-D элемент и четыре узла на 3-D элемент. Квадратичные сетки содержат угловые узлы и узлы в середине каждого ребра элемента. Для квадратичных сеток существует шесть узлов на 2-D элемент и 10 узлов на 3-D элемент.
Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.
pdemesh(model,'NodeLabels','on')