interpolateMagneticFlux

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

    Описание

    пример

    Bintrp = interpolateMagneticFlux(magnetostaticresults,xq,yq) возвращает интерполированную плотность магнитного потока в точках, заданных в xq и yq.

    пример

    Bintrp = interpolateMagneticFlux(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);

    Сгенерируйте mesh.

    generateMesh(emagmodel);

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

    R = solve(emagmodel);
    pdeplot(emagmodel,'FlowData',[R.MagneticFluxDensity.Bx ...
                                  R.MagneticFluxDensity.By])
    axis equal

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

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

    v = linspace(-0.5,0.5,51);
    [X,Y] = meshgrid(v);
    Bintrp = interpolateMagneticFlux(R,X,Y)
    Bintrp = 
      FEStruct with properties:
    
        Bx: [2601x1 double]
        By: [2601x1 double]
    
    

    Измените Bintrp.Bx и Bintrp.By и постройте получившуюся плотность магнитного потока.

    BintrpX = reshape(Bintrp.Bx,size(X));
    BintrpY = reshape(Bintrp.By,size(Y));
    figure
    quiver(X,Y,BintrpX,BintrpY,'Color','red')

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

    В качестве альтернативы можно задать сетку при помощи матрицы точек запроса.

    querypoints = [X(:),Y(:)]';
    Bintrp = interpolateMagneticFlux(R,querypoints);

    Входные параметры

    свернуть все

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

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

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

    interpolateMagneticFlux преобразует точки запроса в вектор-столбцы xq(:) и yq(:). Это возвращает плотность магнитного потока как вектор-столбец, одного размера. Чтобы гарантировать, что размерности возвращенного решения сопоставимы с размерностями точек исходного запроса, использовать reshape. Например, используйте BintrpX = reshape(Bintrp.Bx,size(xq)).

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

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

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

    interpolateMagneticFlux преобразует точки запроса в вектор-столбцы xq(:) и yq(:). Это возвращает плотность магнитного потока как вектор-столбец, одного размера. Чтобы гарантировать, что размерности возвращенного решения сопоставимы с размерностями точек исходного запроса, использовать reshape. Например, используйте BintrpY = reshape(Bintrp.By,size(yq)).

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

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

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

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

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

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

    свернуть все

    Плотность магнитного потока в точках запроса, возвращенных как FEStruct объект со свойствами, представляющими пространственные компоненты плотности магнитного потока в точках запроса. Для точек запроса, которые находятся вне геометрии, Bintrp.Bx(i) и Bintrp.By(i) NaN. Свойства FEStruct объект только для чтения.

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