Процесс распространения гармонических волн расширения и сдвига в окрестности кругового отверстия
Дипломная работа - Физика
Другие дипломы по предмету Физика
p>}I(int n, double x)
{res=0, add=100;m=0;
{=(pow(-1, m)*pow(x/2,(2*m)+n))/(G(m+1)*G(m+n+1));+=add;++;
}(fabs(add)>eps);res;
}H(int n, double x)
{res(0,0);= sqrt(2/(Pi*x))*(complex(0,(x-((double)(n+1))*(Pi/4))));res;
}D(double n, double r)
{complex((n*n + n - 0.5*beta*beta*r*r)*I(n, alpha*r)
(alpha*r)*I(n-1,alpha*r),0);
}E(double n, double r)
{complex(n*(n+1)*I(n, alpha*r) - n*alpha*r*I(n+1,alpha*r),0);
}F(double n, double r)
{complex(-(n*n+n-alpha*alpha*r*r+0.5*beta*beta*r*r)*I(n,alpha*r)-alpha*r*I(n-1,alpha*r),0);
}Dm(double n, double r)
{(n*n + n - 0.5*beta*beta*r*r)*H(n, alpha*r)
(alpha*r)*H(n-1,alpha*r);
}e(double n, double r)
{n*(n+1)*H(n, alpha*r) - n*alpha*r*H(n+1,alpha*r);
}Fm(double n, double r)
{-(n*n+n-alpha*alpha*r*r+0.5*beta*beta*r*r)*H(n,alpha*r)-alpha*r*H(n-1,alpha*r);
}K(double n, double r)
{-n*(n+1)*H(n,beta*r)+n*beta*r*H(n-1,beta*r);
}Km(double n, double r)
{-(n*n+n-0.5*beta*beta*r*r)*H(n, beta*r)+beta*r*H(n-1,beta*r);
}det2(complex a[4])
{(a[0]*a[3])-(a[1]*a[2]);
}det3(complex a[9])
{t[4];res(0,0);[0]=a[4];[1]=a[5];[2]=a[7];[3]=a[8];+=a[0]*det2(t);[0]=a[3];[1]=a[5];[2]=a[6];[3]=a[8];=a[1]*det2(t);[0]=a[3];[1]=a[4];[2]=a[6];[3]=a[7];+=a[2]*det2(t);res;
}det(double n, double r)
{t[4];[0]=D(n,r);[1]=K(n,r);[2]=e(n,r);[3]=Km(n,r);det2(t);
}detrr(double n, double r)
{t[9];[0]=Dm(n,r);[1]=D(n,r);[2]=K(n,r);[3]=Dm(n,r);[4]=D(n,r);[5]=K(n,r);[6]=E(n,r);[7]=e(n,r);[8]=Km(n,r);det2(t);
}detrt(double n, double r)
{t[9];[0]=E(n,r);[1]=e(n,r);[2]=Km(n,r);[3]=Dm(n,r);[4]=D(n,r);[5]=K(n,r);[6]=E(n,r);[7]=e(n,r);[8]=Km(n,r);det2(t);
}dettt(double n, double r)
{t[9];[0]=F(n,r);[1]=Fm(n,r);[2]=-K(n,r);[3]=D(n,r);[4]=D(n,r);[5]=K(n,r);[6]=E(n,r);[7]=e(n,r);[8]=Km(n,r);det2(t);
}en(int n)
{(n==0)1;2;
}drr(double r, double teta, double t)
{res=0, add=100, prev=0;n=0;
{=((2*mu*A)/(r*r))*en(n)*pow(complex(0,1),n)*(detrr(n,R)/det(n,R))*cos(teta*n)*exp(complex(0,(-omega)*t));+=add;++;
}(n<steps);res;
}drt(double r, double teta, double t)
{res=0, add=100, prev=0;n=0;
{=((2*mu*A)/(r*r))*en(n)*pow(complex(0,1),n)*(detrt(n,R)/det(n,R))*cos(teta*n)*exp(complex(0,(-omega)*t));+=add;++;
}(n<steps);res;
}dtt(double r, double teta, double t)
{res=0, add=100, prev=0;n=0;
{=((2*mu*A)/(r*r))*en(n)*pow(complex(0,1),n)*(dettt(n,R)/det(n,R))*cos(teta*n)*exp(complex(0,(-omega)*t));+=add;++;
}(n<steps);res;
}main()
{();h = 0.1, rho = 100000, c1 = 2, c2 = 1;=2;=1;= 0.1;= omega/c1;= omega/c2;= (omega * (h/(2*Pi)))*(omega * (h/(2*Pi)))*rho;(int r=1;r<3;r++)
{(int t=0;t<3;t++)
{("--------------- t=%d r=%d ---------------\n",t,r);(" d_r_r = .7f .7fi\n", real(drr(r,Pi,t)), imag(drr(r,Pi,t)));(" d_r_teta = .7f .7fi\n", real(drt(r,Pi,t)), imag(drt(r,Pi,t)));("d_teta_teta = .7f .7fi\n", real(dtt(r,Pi,t)), imag(dtt(r,Pi,t)));
}
}
}
Приложение 2
-------------- t=0 r=1 ---------------_r_r = -0.2862360 -2.8618209i_r_teta = 0.0084582 -2.0247672i_teta_teta = -7.7657672 -1.5775293i
-------------- t=1 r=1 ---------------_r_r = -0.5705114 -2.8189478i_r_teta = -0.1937235 -2.0154963i_teta_teta = -7.8844609 -0.7943651i
-------------- t=2 r=1 ---------------_r_r = -0.8490864 -2.7479087i_r_teta = -0.3939695 -1.9860871i_teta_teta = -7.9243756 -0.0032640i
-------------- t=0 r=2 ---------------_r_r = -0.0715590 -0.7154552i_r_teta = 0.0021146 -0.5061918i_teta_teta = -1.9414418 -0.3943823i
-------------- t=1 r=2 ---------------_r_r = -0.1426278 -0.7047370i_r_teta = -0.0484309 -0.5038741i_teta_teta = -1.9711152 -0.1985913i
-------------- t=2 r=2 ---------------_r_r = -0.2122716 -0.6869772i_r_teta = -0.0984924 -0.4965218i
d_teta_teta = -1.9810939 -0.0008160i