利用matlab绘制骨架曲线与滞回曲线并求解累积滞回耗能刚度退化和强度退化.zip

clear;clc format long H1=-5e5:100:(5e5);H1=[H1 -H1]; H_Oe=H1; H=H1*79.5775;%转换磁场单位为国际单位制 L=length(H); Ms=1.74e5;%饱和磁化强度 Hk=14;%学号后两位 Ku1=Hk*Ms/2; u0=4*pi*1e-7; phi1=(0:15:90)/180*pi;%磁场的方向 lphi=length(phi1); for k1=1:lphi phi=phi1(k1); x0_Res(1)=phi+pi;%设置初始磁场方向为正 M_Res(1)=1;%设置初始磁矩沿着磁场方向，磁化强度为1 for i=2:L i; H0=H(i); x1=-2*pi:pi/1000:2*pi; Lx1=length(x1); E=Ku1*sin(x1).^2-u0*Ms*H0.*cos(phi-x1); %SW模型下的能量公式E=Ku*cos(x1).^2-u0*Ms*H0.*cos(phi-x1); k=1; x0_solve=0; for j=2:(Lx1-1) if E(j-1)>E(j) & E(j+1)>E(j) x0_solve(k)=x1(j); k=k+1;%取到使得能量最小的x1值（θ） end end xabs=abs(x0_Res(i-1)-x0_solve);%磁场方向和磁化强度夹角 x0position=find(xabs==min(xabs));%夹角最小值 x0_Res(i)=x0_solve(x0position(1)); M1_Res(i)=1*cos(phi-x0_Res(i));% 沿磁场方向的磁化强度 end subplot(3,3,k1) plot(H_Oe(2:L),M1_Res(2:L)) axis([-5e5 5e5 -1 1]) end