www.pudn.com > pinsanquxian.rar > lamepA.asv, change:2004-12-16,size:1175b


%求解Lame波反对称部分的基本模态Ao 
format long 
u=0.3; %泊松比u,越大对应 
E=1.88e10; %弹性模量E,越大对应速度越大 
p=1800; 
cd=((1-u)*E/(1+u)/(1-2*u)/p)^0.5;%cd为纵波波速,cs为横波波速 
cs=(E/2/(1+u)/p)^0.5;  
cp=[]; 
err=[]; 
for kh=0:0.1:2   %k为波数,h为总的板厚 
figure(fix(10*kh+1)); 
syms c;   %定义速度c为符号函数    
y1=4*tanh(kh/2*(1-(c/cs)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2); 
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cd)^2)^(1/2)); 
ezplot(y1,[0,2000]); %作图求解 
hold on 
grid on 
ezplot(y2,[0,2000]); 
[xx,yy]=ginput(1); 
[x(1),fval,exitflag]=fzero(inline(y1-y2),xx(1),[]); 
c=x(1); 
y1=4*tanh(kh/2*(1-(c/cs)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2); 
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cd)^2)^(1/2)); 
m=y1-y2; 
cp=[cp;c]; 
err=[err;m]; 
end 
%计算Ao群速度 
kh=[0:0.1:2];  
for i=2:1:21   
    cg(i)=cp(i)+kh(i)*(cp(i)-cp(i-1))/(kh(i)-kh(i-1)); %波的相速度cp与群速度cg的转换公式 
end 
cg(1)=0; 
%作图 
fh=[]; 
for kh=0:0.1:2 
    f=kh*cp(fix(10*kh+1))/2/pi; %将kd与cp、cg的关系转换为fh与cp、cg的关系 
    fh=[fh;f]; 
end     
figure(22)     %fh与cp的图 
plot(fh/2,cp); 
grid on;hold on 
figure(23)     %fh与cg的图 
plot(fh/2,cg); 
grid on;hold on