Магнитостатическое решение и выведенные величины
A MagnetostaticResults объект содержит магнитный потенциал, магнитное поле и значения плотности магнитного потока в форме, удобной для графического изображения и постобработки.
Магнитный потенциал, магнитное поле и плотность магнитного потока вычисляются в узлах треугольного или четырехгранного mesh, сгенерированной generateMesh. Значения магнитного потенциала в узлах появляются в MagneticPotential свойство. Значения магнитного поля в узлах появляются в MagneticField свойство. Магнитная плотность потока в узлах появляется в MagneticFluxDensity свойство.
Чтобы интерполировать магнитный потенциал, магнитное поле и магнитную плотность потока в пользовательскую сетку, такую как та, которая задана как meshgrid, использовать interpolateMagneticPotential, interpolateMagneticField, и interpolateMagneticFlux функций.
Решите магнитостатическую задачу, используя solve функция. Эта функция возвращает решение как MagnetostaticResults объект.
interpolateMagneticPotential | Интерполяция магнитного потенциала в магнитостатическом результате в произвольных пространственных местоположениях |
interpolateMagneticField | Интерполяция магнитного поля в магнитостатическом результате в произвольных пространственных местоположениях |
interpolateMagneticFlux | Интерполяция плотности магнитного потока в магнитостатическом результате в произвольных пространственных местоположениях |
Решите 2-D электромагнитную задачу.
Создайте электромагнитную модель для магнитостатического анализа.
emagmodel = createpde('electromagnetic','magnetostatic');
Импортируйте и постройте график геометрии, представляющей пластину с отверстием.
importGeometry(emagmodel,'PlateHolePlanar.stl')ans =
DiscreteGeometry with properties:
NumCells: 0
NumFaces: 1
NumEdges: 5
NumVertices: 5
Vertices: [5x3 double]
pdegplot(emagmodel,'EdgeLabels','on')
Задайте значение вакуумной проницаемости в системе модулей СИ.
emagmodel.VacuumPermeability = 1.2566370614E-6;
Задайте относительную проницаемость материала.
electromagneticProperties(emagmodel,'RelativePermeability',5000);Примените граничные условия магнитного потенциала к ребрам, обрамляющим прямоугольник и окружность.
electromagneticBC(emagmodel,'MagneticPotential',0,'Edge',1:4); electromagneticBC(emagmodel,'MagneticPotential',0.01,'Edge',5);
Задайте плотность тока для всей геометрии.
electromagneticSource(emagmodel,'CurrentDensity',0.5);Сгенерируйте mesh.
generateMesh(emagmodel);
Решить модель.
R = solve(emagmodel)
R =
MagnetostaticResults with properties:
MagneticPotential: [1218x1 double]
MagneticField: [1x1 FEStruct]
MagneticFluxDensity: [1x1 FEStruct]
Mesh: [1x1 FEMesh]
Постройте график магнитного потенциала и поля.
pdeplot(emagmodel,'XYData',R.MagneticPotential, ... 'FlowData',[R.MagneticField.Hx ... R.MagneticField.Hy]) axis equal
