📄 xiangyang26.cpp
字号:
if(loop1>=0&&loop1<=NNN-1)
{ sgmx=sgmmax*(NNN-loop1-1)/NNN*(NNN-loop1-1)/NNN;
sgmmx=sgmx*142122.3034;
}
else if(loop1>=Nx-NNN&&loop1<Nx)
{ sgmx=sgmmax*(NNN+loop1-Nx)/NNN*(NNN+loop1-Nx)/NNN;
sgmmx=sgmx*142122.3034;
}
else
{sgmx=0;sgmmx=0;}
if(loop2>=0&&loop2<=NNN-1)
{ sgmy=sgmmax*(NNN-loop2-1)/NNN*(NNN-loop2-1)/NNN;
sgmmy=sgmy*142122.3034;
}
else if(loop2>=Ny-NNN&&loop2<Ny)
{ sgmy=sgmmax*(NNN+loop2-Ny)/NNN*(NNN+loop2-Ny)/NNN;
sgmmy=sgmy*142122.3034;
}
else
{sgmy=0;sgmmy=0;}
if(loop3>=0&&loop3<=NNN-1)
{ sgmz=sgmmax*(NNN-loop3-1)/NNN*(NNN-loop3-1)/NNN;
sgmmz=sgmz*142122.3034;
}
else if(loop3>=Nz-NNN&&loop3<Nz)
{ sgmz=sgmmax*(NNN+loop3-Nz)/NNN*(NNN+loop3-Nz)/NNN;
sgmmz=sgmz*142122.3034;
}
else
{sgmz=0;sgmmz=0;}
if(loop1>=Nx/2-DD2/2-1&&loop1<Nx/2-1&&loop2==Ny/2+movedis&&loop3==Nz/2)
{Exy[loop1][loop2][loop3]=0;Exz[loop1][loop2][loop3]=0;}
else if(loop1>Nx/2&&loop1<=Nx/2+DD2/2&&loop2==Ny/2+movedis&&loop3==Nz/2)
{Exy[loop1][loop2][loop3]=0;Exz[loop1][loop2][loop3]=0;}
else
{
Exy[loop1][loop2][loop3]=((2*e-sgmy*dt)/(2*e+sgmy*dt))*Exy[loop1][loop2][loop3]+(2*dt)/(2*e+sgmy*dt)*((Hzx[loop1][loop2][loop3]+Hzy[loop1][loop2][loop3]-Hzx[loop1][loop2-1][loop3]-Hzy[loop1][loop2-1][loop3])/dy);
Exz[loop1][loop2][loop3]=((2*e-sgmz*dt)/(2*e+sgmz*dt))*Exz[loop1][loop2][loop3]-(2*dt)/(2*e+sgmz*dt)*((Hyx[loop1][loop2][loop3]+Hyz[loop1][loop2][loop3]-Hyz[loop1][loop2][loop3-1]-Hyx[loop1][loop2][loop3-1])/dz);
}
}
if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
{ Exy[loop1][loop2][loop3]+= Exz[loop1][loop2][loop3];
Exz[loop1][loop2][loop3]=0;
}
e=ee;
}
}
}
/***********************************************************************************************************************************************************************************************************************************/
/*y向电场*/
for(loop1=1;loop1<Nx;loop1++)
{for(loop2=0;loop2<Ny;loop2++)
{for(loop3=1;loop3<Nz;loop3++)
/* if(loop1==N/2&&loop2==N/2&&loop3==N/2)
{ Eyz[loop1][loop2][loop3]=((2*e-sgmz*dt)/(2*e+sgmz*dt))*Eyz[loop1][loop2][loop3]+(2*dt)/(2*e+sgmz*dt)*((Hxy[loop1][loop2][loop3]+Hxz[loop1][loop2][loop3]-Hxy[loop1][loop2][loop3-1]-Hxz[loop1][loop2][loop3-1])/dz)-sin(2*Pi*f*dt*(M-0.5));
Eyx[loop1][loop2][loop3]=((2*e-sgmx*dt)/(2*e+sgmx*dt))*Eyx[loop1][loop2][loop3]-(2*dt)/(2*e+sgmx*dt)*((Hzx[loop1][loop2][loop3]+Hzy[loop1][loop2][loop3]-Hzx[loop1-1][loop2][loop3]-Hzy[loop1-1][loop2][loop3])/dx);
}
else
*/
{
if(loop1>=0&&loop1<=NNN-1)
{ sgmx=sgmmax*(NNN-loop1-1)/NNN*(NNN-loop1-1)/NNN;
sgmmx=sgmx*142122.3034;
}
else if(loop1>=Nx-NNN&&loop1<Nx)
{ sgmx=sgmmax*(NNN+loop1-Nx)/NNN*(NNN+loop1-Nx)/NNN;
sgmmx=sgmx*142122.3034;
}
else
{sgmx=0;sgmmx=0;}
if(loop2>=0&&loop2<=NNN-1)
{ sgmy=sgmmax*(NNN-loop2-1)/NNN*(NNN-loop2-1)/NNN;
sgmmy=sgmy*142122.3034;
}
else if(loop2>=Ny-NNN&&loop2<Ny)
{ sgmy=sgmmax*(NNN+loop2-Ny)/NNN*(NNN+loop2-Ny)/NNN;
sgmmy=sgmy*142122.3034;
}
else
{sgmy=0;sgmmy=0;}
if(loop3>=0&&loop3<=NNN-1)
{ sgmz=sgmmax*(NNN-loop3-1)/NNN*(NNN-loop3-1)/NNN;
sgmmz=sgmz*142122.3034;
}
else if(loop3>=Nz-NNN&&loop3<Nz)
{ sgmz=sgmmax*(NNN+loop3-Nz)/NNN*(NNN+loop3-Nz)/NNN;
sgmmz=sgmz*142122.3034;
}
else
{sgmz=0;sgmmz=0;}
Eyz[loop1][loop2][loop3]=((2*e-sgmz*dt)/(2*e+sgmz*dt))*Eyz[loop1][loop2][loop3]+(2*dt)/(2*e+sgmz*dt)*((Hxy[loop1][loop2][loop3]+Hxz[loop1][loop2][loop3]-Hxy[loop1][loop2][loop3-1]-Hxz[loop1][loop2][loop3-1])/dz);
Eyx[loop1][loop2][loop3]=((2*e-sgmx*dt)/(2*e+sgmx*dt))*Eyx[loop1][loop2][loop3]-(2*dt)/(2*e+sgmx*dt)*((Hzx[loop1][loop2][loop3]+Hzy[loop1][loop2][loop3]-Hzx[loop1-1][loop2][loop3]-Hzy[loop1-1][loop2][loop3])/dx);
if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
{Eyz[loop1][loop2][loop3]+= Eyx[loop1][loop2][loop3];
Eyx[loop1][loop2][loop3]=0;}
e=ee;
}
}
}
/****************************************************************************************************************************************************************************************************************************************/
/*z向电场*/
for(loop1=0;loop1<Nx+1;loop1++)
{for(loop2=0;loop2<Ny+1;loop2++)
{for(loop3=0;loop3<Nz;loop3++)
{if(loop1==Nx/2&&loop2==Ny/2&&(loop3==Nz/2||loop3==Nz/2-1)) //前一个5是预留网格,后一个5是介质厚度
{
Ezx[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezx[loop1][loop2][loop3]+(2*dt)/(2*e*er+sig*dt)*((Hyz[loop1][loop2][loop3]+Hyx[loop1][loop2][loop3]-Hyz[loop1-1][loop2][loop3]-Hyx[loop1-1][loop2][loop3])/dx)+sin(2*Pi*f*(M*dt-0.5*dt));
// Ezx[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezx[loop1][loop2][loop3]+(2*dt)/(2*e*er+sig*dt)*((Hyz[loop1][loop2][loop3]+Hyx[loop1][loop2][loop3]-Hyz[loop1-1][loop2][loop3]-Hyx[loop1-1][loop2][loop3])/dx)-(sin(2*Pi*(M-0.5)*f*dt)+sin(2*Pi*(M-0.5)*(f+20000000)*dt)+sin(2*Pi*(M-0.5)*(f+40000000)*dt)+sin(2*Pi*(M-0.5)*(f+60000000)*dt)+sin(2*Pi*(M-0.5)*(f+80000000)*dt)+sin(2*Pi*(M-0.5)*(f+100000000)*dt)+sin(2*Pi*(M-0.5)*(f+120000000)*dt)+sin(2*Pi*(M-0.5)*(f+140000000)*dt)+sin(2*Pi*(M-0.5)*(f+160000000)*dt)+sin(2*Pi*(M-0.5)*(f+180000000)*dt)+sin(2*Pi*(M-0.5)*(f+200000000)*dt)+sin(2*Pi*(M-0.5)*(f+220000000)*dt)+sin(2*Pi*(M-0.5)*(f+240000000)*dt)+sin(2*Pi*(M-0.5)*(f+260000000)*dt)+sin(2*Pi*(M-0.5)*(f+280000000)*dt)+sin(2*Pi*(M-0.5)*(f+300000000)*dt));
Ezy[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezy[loop1][loop2][loop3]-(2*dt)/(2*e*er+sig*dt)*((Hxy[loop1][loop2][loop3]+Hxz[loop1][loop2][loop3]-Hxy[loop1][loop2-1][loop3]-Hxz[loop1][loop2-1][loop3])/dy);
}
/* else if(loop1==Nx/2&&loop2==Ny/2&&loop3==Nz/2-D/2-DD/4) //前一个5是预留网格,后一个5是介质厚度
{
Ezx[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezx[loop1][loop2][loop3]+(2*dt)/(2*e*er+sig*dt)*((Hyz[loop1][loop2][loop3]+Hyx[loop1][loop2][loop3]-Hyz[loop1-1][loop2][loop3]-Hyx[loop1-1][loop2][loop3])/dx)+sin(2*Pi*f*(M*dt-0.5*dt));
// Ezx[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezx[loop1][loop2][loop3]+(2*dt)/(2*e*er+sig*dt)*((Hyz[loop1][loop2][loop3]+Hyx[loop1][loop2][loop3]-Hyz[loop1-1][loop2][loop3]-Hyx[loop1-1][loop2][loop3])/dx)-(sin(2*Pi*(M-0.5)*f*dt)+sin(2*Pi*(M-0.5)*(f+20000000)*dt)+sin(2*Pi*(M-0.5)*(f+40000000)*dt)+sin(2*Pi*(M-0.5)*(f+60000000)*dt)+sin(2*Pi*(M-0.5)*(f+80000000)*dt)+sin(2*Pi*(M-0.5)*(f+100000000)*dt)+sin(2*Pi*(M-0.5)*(f+120000000)*dt)+sin(2*Pi*(M-0.5)*(f+140000000)*dt)+sin(2*Pi*(M-0.5)*(f+160000000)*dt)+sin(2*Pi*(M-0.5)*(f+180000000)*dt)+sin(2*Pi*(M-0.5)*(f+200000000)*dt)+sin(2*Pi*(M-0.5)*(f+220000000)*dt)+sin(2*Pi*(M-0.5)*(f+240000000)*dt)+sin(2*Pi*(M-0.5)*(f+260000000)*dt)+sin(2*Pi*(M-0.5)*(f+280000000)*dt)+sin(2*Pi*(M-0.5)*(f+300000000)*dt));
Ezy[loop1][loop2][loop3]=((2*e*er-sig*dt)/(2*e*er+sig*dt))*Ezy[loop1][loop2][loop3]-(2*dt)/(2*e*er+sig*dt)*((Hxy[loop1][loop2][loop3]+Hxz[loop1][loop2][loop3]-Hxy[loop1][loop2-1][loop3]-Hxz[loop1][loop2-1][loop3])/dy);
}*/
else
{
if(loop1>=0&&loop1<=NNN-1)
{ sgmx=sgmmax*(NNN-loop1-1)/NNN*(NNN-loop1-1)/NNN;
sgmmx=sgmx*142122.3034;
}
else if(loop1>=Nx-NNN&&loop1<Nx)
{ sgmx=sgmmax*(NNN+loop1-Nx)/NNN*(NNN+loop1-Nx)/NNN;
sgmmx=sgmx*142122.3034;
}
else
{sgmx=0;sgmmx=0;}
if(loop2>=0&&loop2<=NNN-1)
{ sgmy=sgmmax*(NNN-loop2-1)/NNN*(NNN-loop2-1)/NNN;
sgmmy=sgmy*142122.3034;
}
else if(loop2>=Ny-NNN&&loop2<Ny)
{ sgmy=sgmmax*(NNN+loop2-Ny)/NNN*(NNN+loop2-Ny)/NNN;
sgmmy=sgmy*142122.3034;
}
else
{sgmy=0;sgmmy=0;}
if(loop3>=0&&loop3<=NNN-1)
{ sgmz=sgmmax*(NNN-loop3-1)/NNN*(NNN-loop3-1)/NNN;
sgmmz=sgmz*142122.3034;
}
else if(loop3>=Nz-NNN&&loop3<Nz)
{ sgmz=sgmmax*(NNN+loop3-Nz)/NNN*(NNN+loop3-Nz)/NNN;
sgmmz=sgmz*142122.3034;
}
else
{sgmz=0;sgmmz=0;}
if(loop1==Nx/2&&loop2==Ny/2&&loop3>=Nz/2-DD1/2-1&&loop3<Nz/2-1)
{Ezx[loop1][loop2][loop3]=0;Ezy[loop1][loop2][loop3]=0;}
else if(loop1==Nx/2&&loop2==Ny/2&&loop3<=Nz/2+DD1/2&&loop3>Nz/2)
{Ezx[loop1][loop2][loop3]=0;Ezy[loop1][loop2][loop3]=0;}
else
{
Ezx[loop1][loop2][loop3]=((2*e-sgmx*dt)/(2*e+sgmx*dt))*Ezx[loop1][loop2][loop3]+(2*dt)/(2*e+sgmx*dt)*((Hyz[loop1][loop2][loop3]+Hyx[loop1][loop2][loop3]-Hyz[loop1-1][loop2][loop3]-Hyx[loop1-1][loop2][loop3])/dx);
Ezy[loop1][loop2][loop3]=((2*e-sgmy*dt)/(2*e+sgmy*dt))*Ezy[loop1][loop2][loop3]-(2*dt)/(2*e+sgmy*dt)*((Hxy[loop1][loop2][loop3]+Hxz[loop1][loop2][loop3]-Hxy[loop1][loop2-1][loop3]-Hxz[loop1][loop2-1][loop3])/dy);
}
}
if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
{Ezx[loop1][loop2][loop3]+= Ezy[loop1][loop2][loop3];
Ezy[loop1][loop2][loop3]=0;}
e=ee;
}
}
}
return(1);
}
/***********************************************************************************************************************************************************************************************************************************/
/*x向磁场*/
int chuanbo2(int M)
{int loop1,loop2,loop3;
for(loop1=0;loop1<Nx+1;loop1++)
{for(loop2=0;loop2<Ny;loop2++)
{for(loop3=0;loop3<Nz;loop3++)
{ if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
{sgmx=0;sgmy=0;sgmz=0;
sgmmx=0;sgmmy=0;sgmmz=0;
}
else
{
if(loop1>=0&&loop1<=NNN-1)
{ sgmx=sgmmax*(NNN-loop1-1)/NNN*(NNN-loop1-1)/NNN;
sgmmx=sgmx*142122.3034;
}
else if(loop1>=Nx-NNN&&loop1<Nx)
{ sgmx=sgmmax*(NNN+loop1-Nx)/NNN*(NNN+loop1-Nx)/NNN;
sgmmx=sgmx*142122.3034;
}
else
{sgmx=0;sgmmx=0;}
if(loop2>=0&&loop2<=NNN-1)
{ sgmy=sgmmax*(NNN-loop2-1)/NNN*(NNN-loop2-1)/NNN;
sgmmy=sgmy*142122.3034;
}
else if(loop2>=Ny-NNN&&loop2<Ny)
{ sgmy=sgmmax*(NNN+loop2-Ny)/NNN*(NNN+loop2-Ny)/NNN;
sgmmy=sgmy*142122.3034;
}
else
{sgmy=0;sgmmy=0;}
if(loop3>=0&&loop3<=NNN-1)
{ sgmz=sgmmax*(NNN-loop3-1)/NNN*(NNN-loop3-1)/NNN;
sgmmz=sgmz*142122.3034;
}
else if(loop3>=Nz-NNN&&loop3<Nz)
{ sgmz=sgmmax*(NNN+loop3-Nz)/NNN*(NNN+loop3-Nz)/NNN;
sgmmz=sgmz*142122.3034;
}
else
{sgmz=0;sgmmz=0;}
}
if(loop1==Nx/2&&loop2==Ny/2&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2-1)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2)))
{
Hxy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hxy[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmy*dt)*(2.0/log(dy/r0))*((Ezx[loop1][loop2+1][loop3]+Ezy[loop1][loop2+1][loop3])/dy);
Hxz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hxz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*((Eyz[loop1][loop2][loop3+1]+Eyx[loop1][loop2][loop3+1]-Eyz[loop1][loop2][loop3]-Eyx[loop1][loop2][loop3])/dz);
}
else if(loop1==Nx/2&&loop2==Ny/2-1&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2-1)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2)))
{
Hxy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hxy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*(2.0/log(dy/r0))*((Ezx[loop1][loop2][loop3]+Ezy[loop1][loop2][loop3])/dy);
Hxz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hxz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*((Eyz[loop1][loop2][loop3+1]+Eyx[loop1][loop2][loop3+1]-Eyz[loop1][loop2][loop3]-Eyx[loop1][loop2][loop3])/dz);
}
else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2+1))&&loop2==Ny/2+movedis&&loop3==Nz/2)
{
Hxy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hxy[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmy*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Ezx[loop1][loop2+1][loop3]*dz+Ezy[loop1][loop2+1][loop3]*dz-Ezx[loop1][loop2][loop3]*dz/2*log(dz/r1)-Ezy[loop1][loop2][loop3]*dz/2*log(dz/r1)));
Hxz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hxz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Eyz[loop1][loop2][loop3+1]*dy+Eyx[loop1][loop2][loop3+1]*dy-Eyz[loop1][loop2][loop3]*dy/2*log(dy/r1)-Eyx[loop1][loop2][loop3]*dy/2*log(dy/r1)));
}
else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2+1))&&loop2==Ny/2+movedis&&loop3==Nz/2-1)
{
Hxy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hxy[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmy*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Ezx[loop1][loop2+1][loop3]*dz+Ezy[loop1][loop2+1][loop3]*dz-Ezx[loop1][loop2][loop3]*dz/2*log(dz/r1)-Ezy[loop1][loop2][loop3]*dz/2*log(dz/r1)));
Hxz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hxz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Eyz[loop1][loop2][loop3+1]*dy/2*log(dy/r1)+Eyx[loop1][loop2][loop3+1]*dy/2*log(dy/r1)-Eyz[loop1][loop2][loop3]*dy-Eyx[loop1][loop2][loop3]*dy));
}
else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2+1))&&loop2==Ny/2-1+movedis&&loop3==Nz/2)
{
Hxy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hxy[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmy*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Ezx[loop1][loop2+1][loop3]*dz/2*log(dz/r1)+Ezy[loop1][loop2+1][loop3]*dz/2*log(dz/r1)-Ezx[loop1][loop2][loop3]*dz-Ezy[loop1][loop2][loop3]*dz));
Hxz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hxz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*1.0/(dy*dz-Pi*r1*r1/4)*((Eyz[loop1][loop2][loop3+1]*dy+Eyx[loop1][loop2][loop3+1]*dy-Eyz[loop1][loop2][loop3]*dy/2*log(dy/r1)-Eyx[loop1][loop2][loop3]*dy/2*log(dy/r1)));
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -