弹性板受激振动响应MATLAB程序

通过MATLAB编程实现四边简支平板分别在受到点激励、线激励和面激励下位移响应。画出了位移幅值在垂直于板平面的三维图。

%%给定时间下的平板振动位移响应
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