24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lxshisan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 7
在线: 10.5小时
虫号: 2626796
注册: 2013-08-31
专业: 金属功能材料

[求助] 塑性变形对铁磁材料磁化过程的影响物理模型

最近看关于塑性变形对铁磁材料的磁化过程的影响物理模型，但是根据相关论文中推导的公式都无法得到其论文中数值模拟得到的图形，磁场强度M计算结果图形趋势差不多，但是数值大小差了10^4，最终导致有效磁通密度B成直线，与文献中的结果相差太大，还请有也做相关研究或者数值模拟方面的大神们指导一下。

CODE:

function [x,y]=lx15()ms=1.585e6;  %饱和磁化强度r1=7e-18;r11=-1e-25;r2=-3.3e-30;r22=2.1e-38;u0=4*pi*(1e-7);G1=0.81e-6;G2=3.8e-12;G3=G1;G4=G2;alphak=0.76;C11=2.863e11;C12=1.41e11;poisson=C12/(C11+C12);yield=C11-2*C12*poisson;G=yield/(2*(1+poisson));b=2.5e-10;d=20e-6;k0_u0=500;a0=600;zetad0=1e12;k0=k0_u0*4*pi*10^(-7);f=0.075;sigmay=235e6;epsilonp=0.0001;sigma=yield*epsilonp+sigmay;zetad=(f*yield*epsilonp/(alphak*G*b)+zetad0^0.5)^2;k=(G1+G2/d)*zetad^(1/2)*k0;a=(G3+G4/d)*zetad^(1/2)*a0;alpha=8.44e-6;%求解x12=linspace(0,7e3,100);   %H=0:7000A/mx13=linspace(7e3,-7e3,200);  %H=7000:-7000A/m;x14=linspace(-7e3,7e3,200);  %H=-7000:7000A/m;[x,y11]=ode45(@ode_fun11,x12,[0.0001,0]); %anhysteretic loop when H=0:7000A/m[x,y12]=ode45(@ode_fun12,x13,[y11(100,1),y11(100,2)]);  %hysteretic process when H=7000:-7000A/m[x,y13]=ode45(@ode_fun11,x14,[y12(200,1),y12(200,2)]);  %hysteretic process when H=-7000:7000A/mMan_coth1=ms*(coth((x13'+alpha*y12(:,1)+3*sigma/u0*((r1+r11*sigma)*y12(:,1)+2*(r2+r22*sigma)*y12(:,1).^3))/a)-a./(x13'+alpha*y12(:,1)+3*sigma/u0*((r1+r11*sigma)*y12(:,1)+2*(r2+r22*sigma)*y12(:,1).^3))); %Man公式计算结果Man_coth2=ms*(coth((x14'+alpha*y13(:,1)+3*sigma/u0*((r1+r11*sigma)*y13(:,1)+2*(r2+r22*sigma)*y13(:,1).^3))/a)-a./(x14'+alpha*y13(:,1)+3*sigma/u0*((r1+r11*sigma)*y13(:,1)+2*(r2+r22*sigma)*y13(:,1).^3)));  %Man公式计算结果figureplot(x13',Man_coth1,x14',Man_coth2)title('Man公式计算结果');figureplot(x13',y12(:,2),'r',x14',y13(:,2),'b');grid on;xlabel('H(A/m)');ylabel('M(A/m)');title('toatal magnetization M');figureplot(x13',y12(:,1),'r',x14',y13(:,1),'b');grid on;xlabel('H(A/m)');ylabel('Man(A/m)');title('anhystersis magnetization Man');He12=x13'+alpha*y12(:,2)+3*sigma/u0*((r1+r11*sigma)*y12(:,2)+2*(r2+r22*sigma)*y12(:,2).^2);He13=x14'+alpha*y13(:,2)+3*sigma/u0*((r1+r11*sigma)*y13(:,2)+2*(r2+r22*sigma)*y13(:,2).^2);figureplot(x13',He12,'r',x14',He13,'b');grid on;xlabel('H(A/m)');ylabel('He(A/m)');title('effective magnetic field He');figureBe1=u0*(x13'+y12(:,2)+3*sigma/u0*((r1+r11*sigma)*y12(:,2)+2*(r2+r22*sigma)*y12(:,2).^2));Be2=u0*(x14'+y13(:,2)+3*sigma/u0*((r1+r11*sigma)*y