Buonasera a tutti, premetto che sono nuovo nel forum. Ho cercato altre discussioni che mi permettessero di risolvere il mio problema ma non ho trovato nulla.
Sto preparando la mia tesi triennale in ingegneria elettrica sul calcolo del punto di lavoro di un circuito magnetico non lineare, pertanto devo trovare analiticamente il punto di intersezione tra un'equazione non lineare (il ciclo di isteresi di un materiale) ed un'equazione lineare (retta di carico del circuito magnetico). Il problema è che non riesco trovare un modo per farlo, in quanto la classica fsolve non mi è di aiuto.
Allego di seguito lo script che ho fin qui realizzato:
clear all
clc
%INIZIALIZZAZIONE DEI DATI IN INGRESSO-------------------------------------
fprintf('\n\nInserimento dei parametri del modello di J-A.\n');
Ms=input('Ms= ');
a=input('a= ');
alpha=input('alpha= ');
k=input('k= ');
c=input('c= ');
fprintf('\n\nInserimento delle dimensioni del circuito magnetico.\n')
Am=input('Am= ');
Ag=input('Ag= ');
lm=input('lm= ');
lg=input('lg= ');
%--------------------------------------------------------------------------
mu0=4.*pi.*1e-7; % Permeabilità magnetica nel vuoto
fprintf('Calcolo dei parametri: \na=%f(A/m), k=%f(A/m), c=%f, Ms=%e(A/m), alpha=%e \n\n',a,k,c,Ms,alpha);
% Valori di campo magnetico H - vettore colonna
H=[0:50000:1500000 1500000:-50000:-1500000 -1500000:50000:1500000]';
M0 = 0; % provino non magnetizzato
[Hmodel,Bmodel] = My_JAmodel(a,k,c,Ms,alpha,H,M0);
%Ciclo di isteresi ottenuto dal modello di J-A sul piano BH (blu) e retta
%di carico (magenta)
figure
hold on
plot(Hmodel./1000, Bmodel,'b','Linewidth',0.4);
xlabel('H [kA/m]');
ylabel('B [T]');
grid;
Hmag=linspace(-1500,0,100);
Bmag=-mu0.*(Ag./Am).*(lm./lg).*Hmag.*1000;
plot(Hmag, Bmag,'m','LineWidth',0.4);
hold off
---------------------------------------fine
"My_JAmodel" è la function per il calcolo dell'equazione non lineare.
Mentre "Bmag=-mu0.*(Ag./Am).*(lm./lg).*Hmag.*1000" è l'equazione lineare.
L'altro problema è che Bmag-Hmag, ed Bmodel-Hmodel, sono "indipendenti" in quanto una dipende dalle dimensioni del circuito magnetico (retta di carico) ed una dipende dal solo materiale (ciclo di isteresi).
Spero si sia capito quello di cui sto parlando. Grazie mille in anticipo!
Allegati:
29168_d999572c0f037edb2fd08b9e781f98a5.png