Exponenta Banner
exponenta event banner

interpolateMagneticPotential

Интерполяция магнитного потенциала в магнитостатическом результате в произвольных пространственных местоположениях

свернуть все на странице

    Синтаксис

    Aintrp = interpolatePotential (магнитостатикрезульты, xq, yq)
    Aintrp = interpolatePotential (магнитостатикрезультаты, точки запроса)

    Описание

    пример

    Aintrp = interpolateMagneticPotential(magnetostaticresults,xq,yq) возвращает интерполированные значения магнитного потенциала в точках, указанных в xq и yq.

    пример

    Aintrp = interpolateMagneticPotential(magnetostaticresults,querypoints) возвращает интерполированные значения магнитного потенциала в точках, указанных в querypoints.

    Примеры

    свернуть все

    Открыть сценарий в реальном времени

    Создайте электромагнитную модель для магнитостатического анализа.

    emagmodel = createpde('electromagnetic','magnetostatic');

    Создайте квадратную геометрию и включите ее в модель. Постройте график геометрии с метками кромок.

    R1 = [3,4,-1,1,1,-1,1,1,-1,-1]';
g = decsg(R1,'R1',('R1')');
geometryFromEdges(emagmodel,g);
pdegplot(emagmodel,'EdgeLabels','on')
xlim([-1.5 1.5])
axis equal

    Figure contains an axes. The axes contains 5 objects of type line, text.

    Укажите проницаемость для вакуума в системе единиц измерения СИ.

    emagmodel.VacuumPermeability = 1.2566370614E-6;

    Укажите относительную проницаемость материала.

    electromagneticProperties(emagmodel,'RelativePermeability',5000);

    Примените граничные условия магнитного потенциала к границам квадрата.

    electromagneticBC(emagmodel,'MagneticPotential',0,'Edge',[1 3]); 
electromagneticBC(emagmodel,'MagneticPotential',0.01,'Edge',[2 4]);

    Укажите текущую плотность для всей геометрии.

    electromagneticSource(emagmodel,'CurrentDensity',0.5);

    Создайте сетку.

    generateMesh(emagmodel);

    Решите модель и постройте график магнитного потенциала.

    R = solve(emagmodel);
pdeplot(emagmodel,'XYData',R.MagneticPotential, ...
                  'Contour','on')
axis equal

    Figure contains an axes. The axes contains 12 objects of type patch, line.

    Интерполяция результирующего магнитного потенциала в сетку, покрывающую центральную часть геометрии, для x и y от -0.5 кому 0.5.

    v = linspace(-0.5,0.5,51);
[X,Y] = meshgrid(v);
Aintrp = interpolateMagneticPotential(R,X,Y)
    Aintrp = 2601×1

    0.0056
    0.0057
    0.0058
    0.0059
    0.0060
    0.0061
    0.0062
    0.0063
    0.0064
    0.0065
      ⋮

    Измениться Aintrp и построить график результирующего магнитного потенциала.

    Aintrp = reshape(Aintrp,size(X));
figure
contourf(X,Y,Aintrp)
colormap(cool)
colorbar

    Figure contains an axes. The axes contains an object of type contour.

    Кроме того, сетку можно задать с помощью матрицы точек запроса.

    querypoints = [X(:),Y(:)]';
Aintrp = interpolateMagneticPotential(R,querypoints);

    Входные аргументы

    свернуть все

    Решение магнитостатической задачи, указанной как MagnetostaticResults объект. Создать magnetostaticresults с использованием solve функция.

    Пример: magnetostaticresults = solve(emagmodel)

    точки запроса координат x, заданные как вещественный массив. interpolateMagneticPotential оценивает магнитный потенциал в координатных точках [xq(i) yq(i)] для каждого i. Из-за этого, xq и yq должно иметь одинаковое количество записей.

    interpolateMagneticPotential преобразует точки запроса в векторы столбцов xq(:) и yq(:). Он возвращает значения магнитного потенциала в виде вектора столбца того же размера. Чтобы убедиться, что измерения возвращенного решения согласуются с измерениями исходных точек запроса, используйте reshape. Например, использовать Aintrp = reshape(Aintrp,size(xq)).

    Пример: xq = [0.5 0.5 0.75 0.75]

    Типы данных: double

    точки запроса координат y, заданные как вещественный массив. interpolateMagneticPotential оценивает магнитный потенциал в координатных точках [xq(i),yq(i)] для каждого i. Из-за этого, xq и yq должно иметь одинаковое количество записей.

    interpolateMagneticPotential преобразует точки запроса в векторы столбцов xq(:) и yq(:). Он возвращает значения магнитного потенциала в виде вектора столбца того же размера. Чтобы убедиться, что измерения возвращенного решения согласуются с измерениями исходных точек запроса, используйте reshape. Например, использовать Aintrp = reshape(Aintrp,size(xq)).

    Пример: yq = [1 2 0 0.5]

    Типы данных: double

    Точки запроса, заданные как вещественная матрица с двумя строками для 2-D геометрии. interpolateMagneticPotential оценивает магнитный потенциал в координатных точках querypoints(:,i) для каждого i, таким образом, каждый столбец querypoints содержит ровно одну 2-D точку запроса.

    Пример: querypoints = [0.5 0.5 0.75 0.75; 1 2 0 0.5]

    Типы данных: double

    Выходные аргументы

    свернуть все

    Магнитный потенциал в точках запроса, возвращенный как FEStruct объект со свойствами, представляющими пространственные компоненты магнитного потенциала в точках запроса. Для точек запроса, которые находятся вне геометрии, Aintrp(i) = NaN. Свойства FEStruct объекты доступны только для чтения.

    См. также

    | | | |

    Представлен в R2021a

    Документация по набору инструментов для дифференциальных уравнений в частных производных

    Поддержка