原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
1、梁弯曲分析，自由振动分析的C程序
程序：
#includemath.h
#includestdio.h
floatzgxj(floatb,floath) /*轴惯性矩计算子程序*/
{
floatI;
I=b*h*h*h/12;
return(I);
}
voidmain()/*主函数*//**/
{
inti,j;
floata[4][4];
floatL,b,h,jz,l; /*设置轴的长度、周截面的高度h、宽带b、轴惯性矩、梁单元长度*/
intN; /*设置节点数量N*/
floatE;/*刚度系数模量*/
printf(请输入刚度系数：);
scanf(%f,E);
printf(\n);
printf(请输入杆的长度：);
scanf(%f,L);
printf(\n);
printf(请输入杆的宽带和高度(数据之间用空格隔开)：);
scanf(%f%f,b,h);
printf(\n);
printf(请需要输入划分的节点数量（建议不要太大20）：);
scanf(%d,N);
printf(\n);
jz=zgxj(b,h);/*计算轴惯性矩的值*/
l=L/(N+1);/*求梁单元长度*/
a[0][0]=12*E*jz/(l*l*l);/*计算刚度矩阵*/
a[0][1]=6*E*jz/(l*l);
a[0][2]=(-12)*E*jz/(l*l*l);
a[0][3]=6*E*jz/(l*l);
a[1][0]=6*E*jz/(l*l);
a[1][1]=4*E*jz/l;
a[1][2]=(-6)*E*jz/(l*l);
a[1][3]=2*E*jz/l;
a[2][0]=(-12)*E*jz/(l*l*l);
a[2][1]=(-6)*E*jz/(l*l);
a[2][2]=12*E*jz/(l*l*l);
a[2][3]=(-6)*E*jz/(l*l);
a[3][0]=6*E*jz/(l*l);
a[3][1]=2*E*jz/(l*l);
a[3][2]=(-6)*E*jz/(l*l);
a[3][3]=4*E*jz/l;
printf(梁单元刚度矩阵K：\n);
for(i=0;i=3;i++)/*输出刚度矩阵*/
{
for(j=0;j=3;j++)
{
printf(%.1f,a[i][j]);
}
printf(\n);
}
}
#includestdio.h
#includemath.h
floatzgxj(floatb,floath) /*轴惯性矩计算子程序*/
{
floatI;
I=b*h*h*h/12;
return(I);
}
（2）求振动频率的程序
floatzizhen(floatE,floatI,floatm,floatL,inti) /*模特分析子程序,用于输出前各阶自振频率*/
{
floatf;
floata[4]={3.516,22.03,61.7,120.9};
f=a[i]*sqrt((E*I/m)/L/L);/*计算前四阶自振频率*/
return(f);/*返回单前频率*/
}
voidmain() /*主函数*/
{
floatL,b,h,iz; /*设置轴的长度、周截面的高度h、宽带b、轴惯性矩*/
inti,j; /*设置节点数量N*/
floatE,m;
floata[4]={3.516,22.03,61.7,120.9},f[4];/*定义影响系数和自振频率的字符*/
printf(请输入刚度系数：);
scanf(%f,E);
printf(\n);
printf(请输入杆的长度：);
scanf(%f,L);
printf(\n);
printf(请输入杆的宽带和高度(数据之间用空格隔开)：);
scanf(%f%f,b,h);
printf(\n);
printf(请输入线密度(单位:kg/立方米)：);
scanf(%f,m);
printf(\n);
iz=zgxj(b,h);
for(i=0;i=3;i++)/*循环计算各阶频率*/
{
f[i]=zizhen(E,iz,m,L,i);
j=i+1;
printf(第%d阶的频率为%frad/s\n,j,f[i]);
}
}
运算结果如图1.1所示
图1.1
文档评论（0）