tiange 发表于 2006-10-7 09:27

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

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

若一个悬臂梁，长度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 发表于 2006-10-7 21:28

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

toes 发表于 2006-10-7 22:43

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

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

tiange 发表于 2006-10-8 08:54

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

附：书中的公式

tiange 发表于 2006-10-8 08:55

书中给的公式

toes 发表于 2006-10-8 18:30

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


function test
L=30;%悬臂梁的度
C1=10;    %假设 C1
a1=1.875/L;
a2=4.694/L;
a3=7.855/L;
x=0:1:30;% 因为L为30，所以x最大只能为30，否则没有意义。

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)) );
%第一阶振型函数
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)) );
%第三阶振型函数
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)) );
%第三阶振型函数

subplot(3,1,1)
plot(x,y1)%绘第一阶振型函数曲线
subplot(3,1,2)
plot(x,y2)%绘第三阶振型函数曲线
subplot(3,1,3)
plot(x,y3)%绘第三阶振型函数曲线


[ 本帖最后由 toes 于 2006-10-8 18:42 编辑 ]
页: [1]
查看完整版本: 大家看看我的这个求悬臂梁 振型函数 的程序 错在那里？