24小时热门版块排行榜

>论坛更新日志 (3944)
>文献求助 (430)
>虫友互识 (344)
>导师招生 (268)
>招聘信息布告栏 (146)
>考博 (145)
>硕博家园 (139)
>休闲灌水 (91)
>论文道贺祈福 (90)
>博后之家 (86)
>论文投稿 (74)
>考研 (67)
>教师之家 (57)
>基金申请 (54)
>公派出国 (45)
>找工作 (44)

返回列表

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

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