Поиск элементов сетки в заданной области
elemIDs = findElements(mesh,'region',RegionType,RegionID)
elemIDs = findElements(mesh,'attached',nodeID)nodeID - идентификатор углового узла. Этот синтаксис игнорирует идентификаторы узлов, расположенных в середине ребер элемента.
Для нескольких узлов задайте nodeID как вектор.
Найти элементы, связанные с геометрической областью.
Создайте модель УЧП.
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')
Найти элементы, расположенные в указанном поле.
Создайте модель УЧП.
model = createpde;
Импортируйте и постройте график геометрии.
importGeometry(model,'PlateHolePlanar.stl');
pdegplot(model)
Сгенерируйте mesh.
mesh = generateMesh(model,'Hmax',2,'Hmin',0.4)
mesh =
FEMesh with properties:
Nodes: [2x386 double]
Elements: [6x172 double]
MaxElementSize: 2
MinElementSize: 0.4000
MeshGradation: 1.5000
GeometricOrder: 'quadratic'
Найдите элементы, расположенные в следующем поле.
Eb = findElements(mesh,'box',[5 10],[10 20]);Выделите эти элементы зеленым цветом на сетчатый график.
figure pdemesh(model) hold on pdemesh(mesh.Nodes,mesh.Elements(:,Eb),'EdgeColor','green')
Найти элементы, расположенные на указанном диске.
Создайте модель УЧП.
model = createpde;
Импортируйте и постройте график геометрии.
importGeometry(model,'PlateHolePlanar.stl');
pdegplot(model)
Сгенерируйте mesh.
mesh = generateMesh(model,'Hmax',2,'Hmin',0.4,'GeometricOrder','linear')
mesh =
FEMesh with properties:
Nodes: [2x107 double]
Elements: [3x172 double]
MaxElementSize: 2
MinElementSize: 0.4000
MeshGradation: 1.5000
GeometricOrder: 'linear'
Найдите элементы, расположенные в радиусе 2 от центра [5,10].
Er = findElements(mesh,'radius',[5 10],2);Выделите эти элементы зеленым цветом на сетчатый график.
figure pdemesh(model) hold on pdemesh(mesh.Nodes,mesh.Elements(:,Er),'EdgeColor','green')
Найдите элементы, присоединенные к указанному угловому узлу.
Создайте модель УЧП.
model = createpde;
Импортируйте и постройте график геометрии.
importGeometry(model,'PlateHolePlanar.stl');
pdegplot(model)
Сгенерируйте линейный треугольный mesh путем установки значения геометрического порядка равным linear. Этот mesh содержит только угловые узлы.
mesh = generateMesh(model,'Hmax',2,'Hmin',0.4, ... 'GeometricOrder','linear');
Найдите узел, ближайший к точке [15; 10].
N_ID = findNodes(mesh,'nearest',[15;10])N_ID = 10
Найти элементы, присоединенные к этому узлу.
En = findElements(mesh,'attached',N_ID)En = 1×3
95 97 98
Выделите узел и элементы зеленого цвета на сетчатом графике.
figure pdemesh(model) hold on plot(mesh.Nodes(1,N_ID),mesh.Nodes(2,N_ID),'or','Color','g', ... 'MarkerFaceColor','g') pdemesh(mesh.Nodes,mesh.Elements(:,En),'EdgeColor','green')
area | Свойства FEMesh | findNodes | meshQuality | volume