%%給定時間下的平板振動位移響應
clc
clear
close?all
E=210E9;
ro=7800;
v=0.3;
dr=0.001;%%阻尼系數
h=0.025;
D=E*h^3/12/（1-v^2）;
a=10;%%板長
b=10;%%板寬
p=a/2;%%激勵點位置
q=b/2;
P=10;
%t=17.5E-3;%%激勵時間
%{
figure（1）
omiga=59*2*pi;
%%t=pi/2/omiga;
[xy]=meshgrid（0:0.01:1）;
w=0;
j=0;
for?nx=1:1:7
for?ny=1:1:7
????omigamn=sqrt（D/ro/h）*（（nx*pi/a）^2+（ny*pi/b）^2）;
???????f=omigamn/2/pi;
????if?f>1350
????????break;
????end
????f;nx;ny;
????j=j+1;
??%s=（4*P/ro/h/a/b）*（1/（omigamn^2-omiga^2））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）;
s=（4*P/ro/h/a/b）*（1/（omigamn^2-omiga^2））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）.*（sin（omiga*t）-sin（omigamn*t）*omiga/omigamn）;
s=（4*P*omiga/ro/h/a/b）*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）.*（dr*h*ro*exp（-（dr*t）/（2*h*ro））*（cosh（（t*（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2））/（h*ro））?-?（h*ro*sinh（（t*（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2））/（h*ro））*（dr/（2*h*ro）?-?（dr^2*h*ro?+?h^3*omiga^2*ro^3?-?h^3*omigamn^2*ro^3）/（dr*h^2*ro^2）））/（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2）））/（dr^2*omiga^2?+?h^2*omiga^4*ro^2?-?2*h^2*omiga^2*omigamn^2*ro^2?+?h^2*omigamn^4*ro^2）?-?（h^2*omiga^2*ro^2*sin（omiga*t）?-?h^2*omigamn^2*ro^2*sin（omiga*t）?+?dr*h*omiga*ro*cos（omiga*t））/（omiga*（dr^2*omiga^2?+?h^2*omiga^4*ro^2?-?2*h^2*omiga^2*omigamn^2*ro^2?+?h^2*omigamn^4*ro^2））;
w=w+s;
end
end
f=omigamn/2/pi;
surfc（xyw）
%}
figure（2）
P=10;
omiga=231.3*2*pi;
%%t=pi/2/omiga;
[xy]=meshgrid（0:0.1:10）;
w=0;
%%w2=0;
j=0;
for?nx=1:1:10
for?ny=1:1:10
????omigamn=sqrt（D/ro/h）*（（nx*pi/a）^2+（ny*pi/b）^2）;
?????f=omigamn/2/pi;
???
if?f>350
????????break;
end
????j=j+1;
?????f;nx;ny;
?????%}
s=（4*P/ro/h/a/b）*（1/（omigamn^2*（1+i*dr）-omiga^2））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）;
s=-（4*P/ro/h/a/b）*（1/（omigamn^2*（1+i*dr）-omiga^2））*（a/nx/pi）*（cos（nx*pi*5.5/a）-cos（nx*pi*4.5/a））*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）;
%s=（4*P/ro/h/a/b）*（1/（omigamn^2*（1+i*dr）-omiga^2））*（a*b/nx/ny/pi^2）*（cos（nx*pi*5.5/a）-cos（nx*pi*4.5/a））*（cos（ny*pi*5.5/b）-cos（ny*pi*4.5/b））*sin（nx*pi*x./a）.*sin（ny*pi*y./b）;
%%s=（4*P*omiga/ro/h/a/b）*（1/（omigamn^2-omiga^2））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）.*（sin（omiga*t）/omiga-sin（omigamn*t）/omigamn）;
%s=（4*P*omiga/ro/h/（a*b））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*sin（ny*pi*y./b）.*（dr*h*ro*exp（-（dr*t）/（2*h*ro））*（cosh（（t*（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2））/（h*ro））?-?（h*ro*sinh（（t*（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2））/（h*ro））*（dr/（2*h*ro）?-?（dr^2*h*ro?+?h^3*omiga^2*ro^3?-?h^3*omigamn^2*ro^3）/（dr*h^2*ro^2）））/（dr^2/4?-?h^2*omigamn^2*ro^2）^（1/2）））/（dr^2*omiga^2?+?h^2*omiga^4*ro^2?-?2*h^2*omiga^2*omigamn^2*ro^2?+?h^2*omigamn^4*ro^2）?-?（h^2*omiga^2*ro^2*sin（omiga*t）?-?h^2*omigamn^2*ro^2*sin（omiga*t）?+?dr*h*omiga*ro*cos（omiga*t））/（omiga*（dr^2*omiga^2?+?h^2*omiga^4*ro^2?-?2*h^2*omiga^2*omigamn^2*ro^2?+?h^2*omigamn^4*ro^2））;
%%u=（4*P*omiga/ro/h/a/b）*（1/（omigamn^2-omiga^2））*sin（nx*pi*p/a）*sin（ny*pi*q/b）*sin（nx*pi*x./a）.*s

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4371??2020-03-23?07:57??彈性板受激振動響應MATLAB程序.m