梁弯曲分析,自由振动分析的c程序

《梁弯曲分析,自由振动分析的c程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1、梁弯曲分析，自由振动分析的c程序程序：★include#includefloatzgxj(floatb,floath)/*轴惯性矩计算子程序*/{floatI;I=b*h*h*h/12;return(l);}voidmain()/*主函数*//**/{inti,j;floata[4][4];floatL,b,h,jz,1;/*设置轴的长度、周截面的高度h、宽带b、轴惯性矩、梁单元长度*/intN;/*设置节点数量N*/floatE;/*刚度系数模量*/printfC请输入刚度系数:

2、〃)；scanf("%f",&E);printf(/z");printfC请输入杆的长度:/z)；scanf("%f〃，&L);printfCW);printfC请输入杆的宽带和高度(数据之间用空格隔开)：〃)；scanff%f%fz/,&b,&h);printf("");printfC请需要输入划分的节点数量(建议不要太大<20)：〃)；scanf&N);printf(〃〃)；jz=zgxj(b,h)；/*计算轴惯性矩的值*/1二L/(N+1);/*求梁单元长度*/a[0][0]=12*E*jz/(1*1

3、*1);/*计算刚度矩阵*/a[0][l]=6*E*jz/(l*l);a[0][2]=(-12)*E*jz/(l*l*l);a[0][3]=6*E*jz/(l*l);a[l][0]=6*E*jz/(l*l);a[l][l]=4*E*jz/l;a[l][2]=(-6)*E*jz/(l*l);a[l][3]=2*E*jz/l;a[2][0]=(-12)*E*jz/(1*1*1)；a⑵[l]=(-6)*E*jz/(l*l);a[2][2]=12*E*jz/(1*1*1);a[2][3]=(-6)*E*jz/(l*l);a[3

4、][0]=6*E*jz/(l*l);a[3][l]=2*E*jz/(l*l);a[3][2]=(-6)*E*jz/(l*l);a[3][3]=4*E*jz/l;printfC梁单元刚度矩阵K：『);for(i=0;i<=3;i++)/*输出刚度矩阵*/{for(j=0;j<=3;j++){printfC%,If",a[i][j]);}printf("");#include#includefloatzgxj(floatb,floath)/*轴惯性矩计算子程序*/{floatI;I二

5、b*h*h*h/12;return(T);}(2)求振动频率的程序/*模特分析子程序，用floatzizhcn(floatE,float1,floatin,floatL,inti)士输出前各阶门振频率*/{floatf;floata[4]={3.516,22.03,61.7,120.9};f=a[i]*sqrt((E*l/m)/L/L);/*计算前四阶H振频率*/return(f);/*返回单前频率*/voidmain()/*主函数*/floatL,b,h,iz;inti,j;floatE,m;/*设置轴的长度、周截面

6、的高度h、宽带b、轴惯性矩*//*设宜节点数量N*/floata[4]={3.516,22.03,61.7,120.9},f[4];/*定义影响系数和自振频率的字符*/printfC请输入刚度系数:3;scanf("%f：&E);printf(〃〃)；printfC请输入杆的长度：”);scanf("%f〃，&L);printf("");printfCig输入杆的宽带和高度(数据Z间用空格隔开)：〃)；scanf("%f%f",&b,&h);printf("");printfC请输入线密度(单位:kg/立方

7、米)：〃)；scanf("%f",&m);printfCW);iz二zgxj(b,h);for(i=0;i<=3;i++)/*循环计算各阶频率*/{f[i]=zizhen(E,iz,m,L,i);j二i+1；printf(/z第%d阶的频率为%frad/s",j,f[i]);}运算结果如图1.1所示taC-J

8、ijF:swlderuanjijnCYuYanbinwwtemp-exe请输入刚度系数，-2请输入杆的长度：1请输入杆的宽带和高度《数据之间用空格隔开”0-12M.1为为为为SKO十8的的的的1234

9、第第第第请输入线密度〈单位沖3立方米＞,78600.000137rad/sQ.0Q0B61rad/30.002411rAd/s0,094724rad/5Pressanykeytocontinue图1.1