Магнитостатический раствор и производные количества
A MagnetostaticResults объект содержит значения магнитного потенциала, магнитного поля и плотности магнитного потока в форме, удобной для построения графика и постобработки.
Магнитный потенциал, магнитное поле и плотность магнитного потока вычисляются в узлах треугольной или тетраэдрической сетки, генерируемой 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);Создайте сетку.
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
