Печать сетки PDE
pdemesh(___, отображает данные сетки или решения с использованием любого из аргументов в предыдущих синтаксисах и одного или нескольких Name,Value)Name,Value аргументы пары.
pdemesh( графика узла PDE или данных треугольника p,e,t,u)u с использованием графика сетки. Функция отображает данные узла, если u является вектором столбца, а данные треугольника, если u - вектор строки.
Если необходимо иметь больший контроль над графиком сетки, используйте pdeplot или pdeplot3D вместо pdemesh.
h = pdemesh(___)
Создайте график сети и отобразите метки узлов и элементов сети.
Создайте модель PDE. Включить геометрию встроенной функции lshapeg. Создание сетки геометрии.
model = createpde; geometryFromEdges(model,@lshapeg); mesh = generateMesh(model);
Постройте график сетки.
pdemesh(model)
Также можно распечатать сетку с помощью mesh в качестве входного аргумента.
pdemesh(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])
Примените граничные условия, задайте коэффициенты и решите PDE.
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);
Постройте график сетки, установив прозрачность равной 0,5.
pdemesh(model,'FaceAlpha',0.5)
Найдите элементы, связанные с геометрической областью.
Создайте модель PDE.
model = createpde;
Включить геометрию встроенной функции lshapeg. Постройте график геометрии.
geometryFromEdges(model,@lshapeg); pdegplot(model,'FaceLabels','on','EdgeLabels','on')
Создайте сетку.
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')
Постройте график сетки для геометрии Г-образной мембраны.
[p,e,t] = initmesh('lshapeg'); [p,e,t] = refinemesh('lshapeg',p,e,t); pdemesh(p,e,t)
Теперь решите уравнение Пуассона 1 по геометрии, определенной 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 - решение PDEPDE раствор, указанный как вектор или матрица.
Пример: results = solvepde(model); u = results.NodalSolution; или u = assempde(model,c,a,f);
mesh - Объект-сетьMesh свойство PDEModel object | вывод generateMeshОбъект-сеть, указанный как Mesh свойство PDEModel объект или как вывод generateMesh.
Пример: model.Mesh
nodes - Узловые координатыУзловые координаты, заданные как матрица 2-by-NumNodes для 2-D сетки и матрица 3-by-NumNodes для 3-D сетки. NumNodes - количество узлов.
elements - Матрица связности элементов в терминах идентификаторов узловМатрица связности элементов в терминах идентификаторов узлов, заданная как матрица Node PerElem-by-NumElements. Node PerElem - количество узлов на элемент. Линейные сети содержат только угловые узлы, поэтому существует три узла на элемент 2-D и четыре узла на элемент 3-D. Квадратичные сетки содержат угловые узлы и узлы в середине каждого ребра элемента. Для квадратичных сеток существует шесть узлов на элемент 2-D и 10 узлов на элемент 3-D.
Укажите дополнительные пары, разделенные запятыми Name,Value аргументы. Name является именем аргумента и Value - соответствующее значение. Name должен отображаться внутри кавычек. Можно указать несколько аргументов пары имен и значений в любом порядке как Name1,Value1,...,NameN,ValueN.
pdemesh(model,'NodeLabels','on')