大家看看我的这个求悬臂梁 振型函数 的程序 错在那里？

tiange

大家看看我的这个求悬臂梁 振型函数 的程序 错在那里？

若一个悬臂梁，长度L=30;
那么，按照《结构动力学》里所给的振型函数公式，编程序求前三阶振型图像，但是怎么就和书中的不一样呢？？

程序（设C1=10）：
>> L=30;  %悬臂梁的度
>> C1=10;    %假设 C1
>> a1=1.875/L;  
>> a2=4.694/L;
>> a3=7.855/L;
>> x=0:1:500;
>> y1=C1*(cos(a1*x)-cosh(a1*x)+(cos(a1*L)+cos(a1*L))*(sin(a1*L)-sin(a1*x))/(sin(a1*L)+sinh(a1*L))); %第一阶振型函数
>> y2=C1*(cos(a2*x)-cosh(a2*x)+(cos(a2*L)+cos(a2*L))*(sin(a2*L)-sin(a2*x))/(sin(a2*L)+sinh(a2*L)));%第三阶振型函数
>> y3=C1*(cos(a3*x)-cosh(a3*x)+(cos(a3*L)+cos(a3*L))*(sin(a3*L)-sin(a3*x))/(sin(a3*L)+sinh(a3*L)));%第三阶振型函数
>> subplot(3,1,1)
>> plot(x,y1)%绘第一阶振型函数曲线
>> subplot(3,1,2)
>> plot(x,y2)%绘第三阶振型函数曲线
>> subplot(3,1,3)
>> plot(x,y3)%绘第三阶振型函数曲线
>>

附件给出了前三阶的振型曲线

21172485

你求的特征值不对，不同的振型函数应该对应不同的特征值。

toes

能不能把书中公式贴出来，怀疑振型函数的正确性，懒得找书了，呵呵。

[ 本帖最后由 toes 于 2006-10-7 23:20 编辑 ]

tiange

但是，我根据书中（《结构动力学》盛宏玉 合肥工业大学出版社 2005）所给出的公式，觉得自己的公式没有错误啊。请大家继续给予帮助
谢谢

附：书中的公式

tiange

toes

你的程序中x的赋值有问题，振型函数也不对，细心一点，呵呵。
我改的程序 test.m 和运行结果如下，请楼主对比一下自己的程序，有错包涵，另外，请自己或版主删除重复帖：

1. 
   
2. function test
   
3. L=30;  %悬臂梁的度
   
4. C1=10;    %假设 C1
   
5. a1=1.875/L;
   
6. a2=4.694/L;
   
7. a3=7.855/L;
   
8. x=0:1:30;% 因为L为30，所以x最大只能为30，否则没有意义。
   
9. 
   
10. y1=C1*( cos(a1*x)-cosh(a1*x) + (cos(a1*L)+cosh(a1*L))*(sinh(a1*x)-sin(a1*x))/(sin(a1*L)+sinh(a1*L)) );
    
11. %第一阶振型函数
    
12. y2=C1*( cos(a2*x)-cosh(a2*x) + (cos(a2*L)+cosh(a2*L))*(sinh(a2*x)-sin(a2*x))/(sin(a2*L)+sinh(a2*L)) );
    
13. %第三阶振型函数
    
14. y3=C1*( cos(a3*x)-cosh(a3*x) + (cos(a3*L)+cosh(a3*L))*(sinh(a3*x)-sin(a3*x))/(sin(a3*L)+sinh(a3*L)) );
    
15. %第三阶振型函数
    
16. 
    
17. subplot(3,1,1)
    
18. plot(x,y1)%绘第一阶振型函数曲线
    
19. subplot(3,1,2)
    
20. plot(x,y2)%绘第三阶振型函数曲线
    
21. subplot(3,1,3)
    
22. plot(x,y3)%绘第三阶振型函数曲线
    

复制代码

[ 本帖最后由 toes 于 2006-10-8 18:42 编辑 ]