吉林大学材料力学课程设计
?n???n?n?n??n?22?7.4?n?2 所以满足疲劳条件
在x?2a右侧Mx??230.54N.m ?max?Mx1?4.5Mpa ?min?0 ?m??a??max?2.25Mpa?n????1?68.4 kW?p2????a????m ?nn?n????n22?7.4?n?2 所以满足疲劳条件
??n? (8)对于8点(x?5a)处
d?54mm k??1.81 ???0.81 k??1.61 ???0.76 ??2.4 ???0.1 此处由弯矩图可知弯矩为0,所以只需校核扭转切应力。Mx??230.54N.m ?Mxmax?W?7.4Mpa ??1min?0 m??a??max?3.7Mpa
p2 ?n????1k??42.6?n?2 所以满足疲劳条件 ????a????m
十、程序设计
1.程序流程图
- 13 -
吉林大学材料力学课程设计
2.程序
#include
#define sigma_1 300e6 #define tao_1 150e6 #define t 15 void main() { int m,i;
float d[t],My[t],Mz[t],Mx[t],k_sigma[t],k_tao[t],ipxl_sigma[t],ipxl_tao[t],beita[t],pusai_tao[t]; float M[t],W[t],Wp[t],sigma_max[t],tao_max[t],n_sigma[t],n_tao[t],n_sigma_tao[t]; printf(\共有几组数据:m=\scanf(\for(i=0;i printf(\请输入第%d组相关数据:\\n\ printf(\圆轴直径d:\printf(\弯矩My:\printf(\弯矩Mz:\printf(\扭矩Mx:\ printf(\正应力的有效应力集中系数k_sigma:\printf(\切应力的有效应力集中系数k_tao:\printf(\正应力的尺寸系数ipxl_sigma:\printf(\切应力的尺寸系数ipxl_tao:\printf(\表面质量系数beita:\printf(\敏感系数pusai_tao:\} d[i]='\\0';My[i]='\\0';Mz[i]='\\0';Mx[i]='\\0';k_sigma[i]='\\0';k_tao[i]='\\0';ipxl_sigma[i]='\\0';ipxl_tao[i]='\\0';beita[i]='\\0';pusai_tao[i]='\\0'; for(i=0;i float wanju(float My,float Mz); float kwjmxs(float d); float knjmxs(float d); float max_zhengyl(float M,float W); float max_qieyl(float Mx,float Wp); float zhengylaqys(float k_sigma,float ipxl_sigma,float beita,float sigma_max); float qieylaqys(float k_tao,float ipxl_tao,float beita,float pusai_tao,float tao_max); float wnzhaqys(float n_sigma,float n_tao); M[i]=wanju(My[i],Mz[i]); W[i]=kwjmxs(d[i]); Wp[i]=knjmxs(d[i]); - 14 - 吉林大学材料力学课程设计 sigma_max[i]=max_zhengyl(M[i],W[i]); tao_max[i]=max_qieyl(Mx[i],Wp[i]); n_sigma[i]=zhengylaqys(k_sigma[i],ipxl_sigma[i],beita[i],sigma_max[i]); n_tao[i]=qieylaqys(k_tao[i],ipxl_tao[i],beita[i],pusai_tao[i],tao_max[i]); n_sigma_tao[i]=wnzhaqys(n_sigma[i],n_tao[i]); if(My[i]==0&&Mz[i]==0) n_sigma_tao[i]=n_tao[i]; if(Mx[i]==0) n_sigma_tao[i]=n_sigma[i]; } M[i]='\\0';W[i]='\\0';Wp[i]='\\0';sigma_max[i]='\\0';tao_max[i]='\\0';n_sigma[i]='\\0';n_tao[i]='\\0';n_sigma_tao[i]='\\0'; i=0; printf(\校核结论如下:\\n\\n\while(n_sigma_tao[i]!='\\0') { printf(\在第%d个截面处\\n\if(My[i]==0&&Mz[i]==0) printf(\ n_sigma_tao[%d]=%g\\n\else printf(\ n_sigma_tao[%d]=%g\\n\if(n_sigma_tao[i]>n) printf(\在第%d个截面处,n_sigma_tao=%f>n=2,满足疲劳强度要求\\n\else printf(\在第%d个截面处,n_sigma_tao=%f printf(\if(i==m) printf(\所校核的截面均满足疲劳强度要求\} /*求弯矩的函数,M=sqrt(My*My+Mz*Mz)*/ float wanju(float My,float Mz) { return(sqrt(My*My+Mz*Mz)); } /*求抗弯截面系数的函数,W=1.0/32*PI*pow(d,3)*/ float kwjmxs(float d) { return(1.0/32*PI*pow(d,3)); } /*求抗扭截面系数的函数,Wp=1.0/16*PI*pow(d,3)*/ - 15 - 吉林大学材料力学课程设计 float knjmxs(float d) { return(1.0/16*PI*pow(d,3)); } /*求最大正应力的函数,sigma_max=M/W*/ float max_zhengyl(float M,float W) { return(M/W); } /*求最大切应力的函数,tao_max=Mx/Wp*/ float max_qieyl(float Mx,float Wp) { return(Mx/Wp); } /*求正应力作用下构件的安全工作因数,n_sigma=sigma_1/(k_sigma*sigma_max/(ipxl_sigma*beita))*/ float zhengylaqys(float k_sigma,float ipxl_sigma,float beita,float sigma_max) { float n_sigma; n_sigma=sigma_1/(k_sigma*sigma_max/(ipxl_sigma*beita)); return(n_sigma); } /*求切应力作用下构件的安全工作因数,n_tao=tao_1/(k_tao*tao_a/(ipxl_tao*beita)+pusai_tao*tao_m)*/ float qieylaqys(float k_tao,float ipxl_tao,float beita,float pusai_tao,float tao_max) { float tao_a,tao_m,n_tao; tao_a=tao_max/2.0; tao_m=tao_max/2.0; n_tao=tao_1/(k_tao*tao_a/(ipxl_tao*beita)+pusai_tao*tao_m); return(n_tao); } /*求弯扭组合作用下构件的安全因数,n_sigma_tao=n_sigma*n_tao/sqrt(n_sigma*n_sigma+n_tao*n_tao)*/ float wnzhaqys(float n_sigma,float n_tao) { float n_sigma_tao; n_sigma_tao=n_sigma*n_tao/sqrt(n_sigma*n_sigma+n_tao*n_tao); return(n_sigma); } - 16 -
相关推荐: