梁彎曲分析,自由振動分析的c程序

《梁彎曲分析,自由振動分析的c程序》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、1、梁彎曲分析，自由振動分析的c程序程序：★include#includefloatzgxj(floatb,floath)/*軸慣性矩計算子程序*/{floatI;I=b*h*h*h/12;return(l);}voidmain()/*主函數(shù)*//**/{inti,j;floata[4][4];floatL,b,h,jz,1;/*設(shè)置軸的長度、周截面的高度h、寬帶b、軸慣性矩、梁單元長度*/intN;/*設(shè)置節(jié)點數(shù)量N*/floatE;/*剛度系數(shù)模量*/printfC請輸入剛度系數(shù):

2、〃)；scanf("%f",&E);printf(/z");printfC請輸入桿的長度:/z)；scanf("%f〃，&L);printfCW);printfC請輸入桿的寬帶和高度(數(shù)據(jù)之間用空格隔開)：〃)；scanff%f%fz/,&b,&h);printf("");printfC請需要輸入劃分的節(jié)點數(shù)量(建議不要太大<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()/*主函數(shù)*/floatL,b,h,iz;inti,j;floatE,m;/*設(shè)置軸的長度、周截面

6、的高度h、寬帶b、軸慣性矩*//*設(shè)宜節(jié)點數(shù)量N*/floata[4]={3.516,22.03,61.7,120.9},f[4];/*定義影響系數(shù)和自振頻率的字符*/printfC請輸入剛度系數(shù):3;scanf("%f：&E);printf(〃〃)；printfC請輸入桿的長度：”);scanf("%f〃，&L);printf("");printfCig輸入桿的寬帶和高度(數(shù)據(jù)Z間用空格隔開)：〃)；scanf("%f%f",&b,&h);printf("");printfC請輸入線密度(單位:kg/立方

7、米)：〃)；scanf("%f",&m);printfCW);iz二zgxj(b,h);for(i=0;i<=3;i++)/*循環(huán)計算各階頻率*/{f[i]=zizhen(E,iz,m,L,i);j二i+1；printf(/z第%d階的頻率為%frad/s",j,f[i]);}運算結(jié)果如圖1.1所示taC-J

8、ijF:swlderuanjijnCYuYanbinwwtemp-exe請輸入剛度系數(shù)，-2請輸入桿的長度：1請輸入桿的寬帶和高度《數(shù)據(jù)之間用空格隔開”0-12M.1為為為為SKO十8的的的的1234

9、第第第第請輸入線密度〈單位沖3立方米＞,78600.000137rad/sQ.0Q0B61rad/30.002411rAd/s0,094724rad/5Pressanykeytocontinue圖1.1