梁彎曲分析自由振動分析的C程序.doc

1、1、梁彎曲分析，自由振動分析的c程序程序：#include <math.h>#include <stdio.h>float zgxj(float b,float h)/*軸慣性矩計算子程序*/float i;i=b*h*h*h/12;return(i);void main() /*主函數(shù)*/*/ int i,j;float a44;float l,b,h,jz,l; /*設置軸的長度、周截面的高度h、寬帶b、軸慣性矩、梁單元長度*/int n;/*設置節(jié)點數(shù)量n*/float e; /*剛度系數(shù)模量*/printf("請輸入剛度系數(shù)：");scanf

2、("%f",&e);printf("n");printf("請輸入桿的長度：");scanf("%f",&l);printf("n");printf("請輸入桿的寬帶和高度(數(shù)據(jù)之間用空格隔開)：");scanf("%f %f",&b,&h);printf("n");printf("請需要輸入劃分的節(jié)點數(shù)量（建議不要太大<20）：");scanf("%d",&

3、amp;n);printf("n");jz=zgxj(b,h);/*計算軸慣性矩的值*/l=l/(n+1);/*求梁單元長度*/a00=12*e*jz/(l*l*l); /*計算剛度矩陣*/a01=6*e*jz/(l*l);a02=(-12)*e*jz/(l*l*l);a03=6*e*jz/(l*l);a10=6*e*jz/(l*l);a11=4*e*jz/l;a12=(-6)*e*jz/(l*l);a13=2*e*jz/l;a20=(-12)*e*jz/(l*l*l);a21=(-6)*e*jz/(l*l);a22=12*e*jz/(l*l*l);a23=(-6)*e*j

4、z/(l*l);a30=6*e*jz/(l*l);a31=2*e*jz/(l*l);a32=(-6)*e*jz/(l*l);a33=4*e*jz/l;printf("梁單元剛度矩陣k：n");for(i=0;i<=3;i+) /*輸出剛度矩陣*/for(j=0;j<=3;j+)printf(" %.1f ",aij);printf("n");#include<stdio.h>#include<math.h>float zgxj(float b,float h)/*軸慣性矩計算子程序*/float i

5、;i=b*h*h*h/12;return(i);（2）求振動頻率的程序float zizhen(float e,float i,float m,float l,int i) /*模特分析子程序,用于輸出前各階自振頻率*/float f;float a4=3.516,22.03,61.7,120.9;f=ai*sqrt(e*i/m)/l/l);/*計算前四階自振頻率*/return(f);/*返回單前頻率*/void main()/*主函數(shù)*/float l,b,h,iz; /*設置軸的長度、周截面的高度h、寬帶b、軸慣性矩*/int i,j;/*設置節(jié)點數(shù)量n*/float e,m;float

6、 a4=3.516,22.03,61.7,120.9,f4;/*定義影響系數(shù)和自振頻率的字符*/printf("請輸入剛度系數(shù)：");scanf("%f",&e);printf("n");printf("請輸入桿的長度：");scanf("%f",&l);printf("n");printf("請輸入桿的寬帶和高度(數(shù)據(jù)之間用空格隔開)：");scanf("%f %f",&b,&h);printf("n");printf("請輸入線密度(單位:kg/立方米)：");scanf("%f",&m);printf("n");iz=zgxj(b,h); for(i=0;i<=3;i+) /*