syms x y;
c=151;%('流向板长度');
l=275;%('板宽');
t=1;%('y向涡格数');
s=1;%('流向涡格数');
oo=4;%('x阶数');
pp=4;%('y阶数');
h=0.381;%('厚度');
r=2.84e-9;%('密度');
v=0.33;%('泊松比');
e=70600;%('杨氏模量');
for m=1:1:s;
for n=1:1:t;
%f1(m)=(sqrt(2))*sin((m-2/3)*pi*x/c+3*pi/4)+exp(-(m-2/3)*pi*x/c)+(power(-1,m+1))*exp(-(m-2/3)*pi*(1-x/c));
f1(m)=(sqrt(2))*sin((m+1)*pi*x/c+3*pi/4)+exp(-(m+1)*pi*x/c)+(power(-1,m+1))*exp(-(m+1)*pi*(1-x/c));
f2(n)=(sqrt(2))*sin((n-0.5)*pi*y/l-pi/4)+exp(-(n-0.5)*pi*y/l)+(power(-1,n+1))*exp(-(n-0.5)*pi*(1-y/l));
%f1(m)=-sin(((2*m+1)*pi/2)*x)+cos(((2*m+1)*pi/2)*x)+log(-((2*m+1)*pi/2)*x)-((-1)^(m))*log(-((2*m+1)*pi/2)*(1-x));
%自由
%f2(n)=sin(((2*n-1)*pi/2)*y)-cos(((2*n-1)*pi/2)*y)+exp(-((2*n-1)*pi/2)*y)+((-1)^(n+1))*exp(-((2*n-1)*pi/2)*(1-y));
%固支
%f1(m)=sin(m*pi*x/c);;%sin(m*pi*x/c)-sin((m+2)*pi*x/c)*(m^3)/(m+2)^3;
%自由i
%f2(n)=sin(n*pi*y/l);;%sin(n*pi*y/l)-sin((n+2)*pi*y/l)*(n^3)/((n+2)^3);%sin(n*pi*y/l);
%固支j
end;
end;
for mm=1:1:oo;
for nn=1:1:pp;
f3(mm)=(sqrt(2))*sin((mm+1)*pi*x/c+3*pi/4)+exp(-(mm+1)*pi*x/c)+(power(-1,mm+1))*exp(-(mm+1)*pi*(1-x/c));
f4(nn)=(sqrt(2))*sin((nn-0.5)*pi*y/l-pi/4)+exp(-(nn-0.5)*pi*y/l)+(power(-1,nn+1))*exp(-(nn-0.5)*pi*(1-y/l));
%f1(m)=-sin(((2*m+1)*pi/2)*x)+cos(((2*m+1)*pi/2)*x)+log(-((2*m+1)*pi/2)*x)-((-1)^(m))*log(-((2*m+1)*pi/2)*(1-x));
%自由
%f2(n)=sin(((2*n-1)*pi/2)*y)-cos(((2*n-1)*pi/2)*y)+exp(-((2*n-1)*pi/2)*y)+((-1)^(n+1))*exp(-((2*n-1)*pi/2)*(1-y));
%固支
%f3(mm)=sin(mm*pi*x/c);%sin(m*pi*x/c)-sin((m+2)*pi*x/c)*(m^3)/(m+2)^3;
%m
%f4(nn)=sin(nn*pi*y/l);;%sin(n*pi*y/l)-sin((n+2)*pi*y/l)*(n^3)/((n+2)^3);%sin(n*pi*y/l);
%n
end;
end;
f55=f2.'*f1;%ij
f5=reshape(f55,1,[]);
f6=f4.'*f3;%mn
fz6=reshape(f6,1,[]);
f7=f5.'*fz6;
m2=int(int(f7,y,0,l),x,0,c);%m矩阵 ?
m=m2*r*h;
f11=diff(f1,x,1);
f12=diff(f1,x,2);
f21=diff(f2,y,1);
f22=diff(f2,y,2);
f31=diff(f3,x,1);
f32=diff(f3,x,2);
f41=diff(f4,y,1);
f42=diff(f4,y,2);
ff55=f2.'*f12;
ff5=reshape(ff55,1,[]);
ff6=f4.'*f32;
ff66=reshape(ff6,1,[]);
ff7=ff5.'*ff66;%m2ni2j
fk55=f22.'*f1;
fk5=reshape(fk55,1,[]);
fk6=f42.'*f3;
fk66=reshape(fk6,1,[]);
fk7=fk5.'*fk66;%mn2ij2
fp55=f2.'*f12;
fp5=reshape(fp55,1,[]);
fp6=f42.'*f3;
fp66=reshape(fp6,1,[]);
fp7=fp5.'*fp66;%mn2i2j
fq55=f22.'*f1;
fq5=reshape(fq55,1,[]);
fq6=f4.'*f32;
fq66=reshape(fq6,1,[]);
fq7=fq5.'*fq66;%m2nij2
fw55=f21.'*f11;
fw5=reshape(fw55,1,[]);
fw6=f41.'*f31;
fw66=reshape(fq6,1,[]);
fw7=fw5.'*fw66;%m1n1i1j1
k1=ff7+fk7+v*(fp7+fq7)+2*(1-v)*fw7;
k2=int(int(k1,y,0,l),x,0,c);
k=(e*h*h*h/12/(1-v*v))*k2;
评论0