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

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

若一个悬臂梁，长度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)%绘第三阶振型函数曲线
>>


附件给出了各阶的振型曲线

[ 本帖最后由 lxq 于 2006-10-7 21:09 编辑 ]

feifeifool 发表于 2006-10-7 16:23

肯定是振型函数的方程写的不对！自己查查看。

fanghuikeer 发表于 2006-10-7 20:37

振型函数不对,找本结构力学书好好看看,画图的程序没有错误...

tiange 发表于 2006-10-8 09:02

我的振型函数好像没有问题啊

我是按照《结构动力学》（盛宏玉，合肥工业大学出版社 2005）中的公式做的
请大家继续给予指点
谢谢

附 书中的公式

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

http://forum.vibunion.com/forum/viewthread.php?tid=27402&page=1&extra=page%3D1#pid149457
页: [1]
查看完整版本: 大家看看我的这个求悬臂梁 振型函数 的程序 错在那里？