材料力学课程设计168

材料力学课程设计

设计题目： 传动轴的静强度、变形

及疲劳强度的计算

学院班级： 汽车学院发动机12班 姓名学号： 王富强（168） 42101203

指导教师: 李峰

目 录

1．材料力学课程设计的目的…………………………..................3

2．材料力学课程设计的任务和要求……………………………..3

3．材料力学课程设计的题目……………………………………..3

3－1 设计题目 …………………………………………………3 3—2传动轴的零件图 ..............................................……………4 3—3设计数据 ........................................................……………4

4.材料力学课程设计的具体设计方案…………………………….4

4—1绘出传动轴的受力简图………………………………………4 4—2 作扭矩图及弯矩图 …………………………………………6 4—3根据强度条件设计等直轴的直径 ……………………………8 4—4计算齿轮处轴的挠度（均按直径Φ1的等直杆计算）…………11 4—5对阶梯传动轴进行疲劳强度计算（若不满足，采取改进措施 使其满足疲劳强度） ……………………………………………13 4—6对所取数据的理论根据作必要的说明 ………………………20 4—7设计所使用软件名称……………………… ……………….20

附录：本题所编写的C程序…………………………………….20

参考文献…………………………………………………………..30

材料力学课程设计的体会和收获………………………………..30

2012 年 9 月 21 日

1．材料力学课程设计的目的

本课程设计的目的是在于系统学完材料力学之后，能结合工程中的实际问题，运用材料力学的基本理论和计算方法，独立地计算工程中的典型零部件，以达到综合运用材料力学的知识解决工程实际问题之目的。同时，可以使学生将材料力学的理论和现代计算方法及手段融为一体。既从整体上掌握了基本理论和现代的计算方法，又提高了分析问题，解决问题的能力；既把以前所学的知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等）综合运用，又为后继课程（机械设计、专业课等）打下基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。具体的有以下六项：

1．使学生的材料力学知识系统化、完整化；

2．在系统全面复习的基础上，运用材料力学知识解决工程中的实际问题；

3．由于选题力求结合专业实际，因而课程设计可以把材料力学知识和专业需要结 合起来；

4．综合运用了以前所学的个门课程的知识（高数、制图、理力、算法语言、计算机等等）使相关学科的知识有机地联系起来；

5．初步了解和掌握工程实践中的设计思想和设计方法； 6．为后继课程的教学打下基础。

2．材料力学课程设计的任务和要求

要求参加设计者，要系统地复习材料力学的全部基本理论和方法，独立分析、判断、设计题目的已知条件和所求问题。画出受力分析计算简图和内力图，列出理论依据和导出计算公式，独立编制计算程序，通过计算机给出计算结果，并完成设计计算说明书。

3．材料力学课程设计的题目

传动轴的强度、变形及疲劳强度计算

3－1 设计题目

传动轴的材料为优质碳素结构钢（牌号45），许用应力[σ]=80MPa，经高频淬火处理，其σb=650MPa，σ-1=300MPa,τ-1=155MPa,磨削轴的表面，键槽均为端铣加工，阶梯轴过渡圆弧r均为2，疲劳安全系数n=2，要求：

1）绘出传动轴的受力简图; 2）作扭矩图及弯矩图；

3）根据强度条件设计等直轴的直径； 4）计算齿轮处轴的挠度；（按直径Φ1的等直杆计算） 5）对阶梯传动轴进行疲劳强度计算；（若不满足，采取改进措施使其满足疲劳强度）； 6）对所取数据的理论根据作必要的说明。

说明 ：

a) 坐标的选取均按下图所示； b) 齿轮上的力F与节圆相切；

c) 数据表中P为直径D的皮带轮传递的功率，P1为直径为D1的皮带轮传递的功率。

3—2传动轴的零件图

Φ1 为静强度条件所确定的轴径，尺寸最后一位数准确到mm，并取偶数。 设

?1?2?3???1.1 ?2?3?4图号如下：

3—3设计数据

本次课程设计采用第16组数据，如下：

P/kW 33.1

P1/kW 16.1 n/(r/min) 1500 D/mm 900 D1/mm 300 D2/mm 220 G2/N 960 G1/N 250 a/mm 300 a(o) 25 4.材料力学课程设计的具体设计方案

4—1绘出传动轴的受力简图

分析传动轴的零件图（下图）和受力图（右图），P为直径D的皮带轮传递的功率，所

P=210.715Nm , P1为直径为D1的皮带轮传递的功nP1率，所以直径D2的皮带轮传递的力矩M1=9549=102.493Nm。

n以直径D的皮带轮传递的力矩M=9549

在传动轴旋转方向上由力矩守衡可得平衡方程 F×D2/2+(2F1-F1)×D1/2+(F2-2F2)×D2/2=0 其中

M=(2F2-F2)D/2 , M1=(2F1-F1)D1/2

故可解得F=2(M-M1)/D2=983.836N

传动轴的受力图：

传动轴的零件图：

现绘出传动轴的受力简图（如下图所示）：

4—2 作扭矩图及弯矩图

由传动轴的受力简图可求支反力得

Fy1=(4Fcos?+2G1+6F1+G2)/5=1825.286N

Fz1=(4Fsin?+3F2)/5=613.583N Fy2=(Fcos?+3G1+9F1+4G2)/5=2326.243N Fz2=(4Fsin?+12F2)/5=1456.441N

并作出传动轴各截面的内力图： 沿y轴方向的剪力图：

沿z轴方向的剪力图：

扭矩图：

沿y轴方向的弯矩图：

沿z轴方向的弯矩图：

4—3根据强度条件设计等直轴的直径

I．由于传动轴的材料为优质碳素结构钢（牌号45），因此需要选用第三强度理论进行强度计算。

根据第三强度理论?r3=?2?4?2?12W Mx?My2?Mz2???? 3其中 W???132

由扭矩图与弯矩图可确定危险截面在D截面右侧与E截面左侧。

在D截面右侧MDy?1387.829Nm，MDz?352.646Nm，MDx?210.715Nm，则有

M22D,max?MDy?MDz?M22Dx?(1387.829Nm)2?(352.646Nm)?(210.715Nm)2?1447.352Nm

在E截面左侧MEy?697.873Nm，MEz?436.932Nm，MEx?210.715Nm，则有

ME,max?MEy?MEz?MEx?(697.873Nm)2?(436.932Nm)2?(210.715Nm)2?1248.903Nm222MD,max?ME,max，所以等直轴只需要满足D截面右侧即可。因此

?D,max?132MD,max??1447.352Nm?????80MPa 3W??1解得?1?56.906mm，取?1?57mm。由

?1?2?3???1.1得 ?2?3?4?2?51.733mm，取?2?52mm；

?3?47.030mm，取?3?48mm； ?4?42.755mm，取?4?43mm；

II．校核?2是否满足静强度条件。 此时需对U截面左侧进行校核。其中W?2???2332；

在U截面左侧MUy?1177.767Nm，MUz?310.504Nm，MUx?108.222Nm，则有 MU,max?MUy?MUz?MUx?(1177.767Nm)2?(310.504Nm)2?(108.222Nm)2?1222.808Nm222因此?U,max?132MU,max??1222.808Nm?????80MPa 3W?2??2解得?2?53.797mm?52mm，所以?2不满足静强度条件。 取?2?54mm，由

?1?2?3???1.1得 ?2?3?4?1?59.176mm，取?1?60mm；

?3?48.906mm，取?3?49mm； ?4?44.460mm，取?4?45mm.

III．校核?3是否满足静强度条件。

此时需对Q截面左侧，V截面右侧和E截面左侧进行校核。很明显MV,max?MQ,max，

其中W?3???3332。

在V截面左侧MVy?1042.851Nm，MVz?394.789Nm，MVx?210.715Nm，则有 MV,max?MVy?MVz?MVx?(1042.815Nm)2?(394.789Nm)2?(210.715Nm)2?1134.812Nm222在E截面左侧MEy?697.873Nm，MEz?436.932Nm，MEx?210.715Nm，则有

ME,max?MEy?MEz?MEx?(697.873Nm)2?(436.932Nm)2?(210.715Nm)2?849.904Nm222MVy?MEy，因此?V,max?132MV,max??1134.812Nm?????80MPa 3W?3??3解得?3?52.474mm?49mm，所以?3不满足静强度条件。 取?3?53mm，由

?1?2?3???1.1得 ?2?3?4?1?63.494mm，取?1?64mm； ?2?57.721mm，取?2?58mm； ?4?47.704mm，取?4?48mm。

IV．校核

?4是否满足静强度条件。

由III.中数据对比可知，?4满足静强度条件。

综上所述，?1?64mm，?2?58mm，?3?53mm，?4?48mm。

4—4计算齿轮处轴的挠度（均按直径Φ1的等直杆计算）

图中直径为D2的轮为齿轮。

I．可以在该轮处（图中B点位置）沿y轴方向加一单位力F=1，并作出单位力作用下的弯矩图M图。 沿y轴方向的弯矩图：

M图：

其中E=200GPa(数据来源：《材料力学》(第二版)（机械工业出版社）P29页表2-2)，

Iy?Iz???1464

此时可以利用图形互乘法求齿轮处该轴沿y轴方向的挠度

fy?1123123[?MBya?M?MBy2a?M??(MDy?MBy)2a?M?MEya?M?EIz2342381511?(MDy?MEy)a?M??MEya?M]?1.144mm21226

II．再在该轮处沿z轴方向加一单位力F=1，并作出单位力作用下的弯矩图M图。 沿z轴方向的弯矩图：

M图：

此时可以利用图形互乘法求齿轮处该轴沿z轴方向的挠度

fz?11251111[?MBza?M?MBz3a?M??(MEz?MBz)3a?M??MEza?M] EIy2382226?0.438mmIII． f?

fy?fz?4.0122?1.8962mm?1.225mm

224—5对阶梯传动轴进行疲劳强度计算（若不满足，采取改进措施使其满足疲劳强度）

I．首先对传动轴键槽进行疲劳强度计算

因为该轴键槽为端铣加工，σb=650MPa，所以根据《材料力学》（第二版）（机械工业出版社）P369页图13-10a可查得K?=1.8，根据《材料力学》（第二版）（机械工业出版社）

P369页图13-10b可查得K?=1.48。

因为该轴经高频淬火处理，σb=650MPa，K?=1.8，所以根据《材料力学》（第二版）（机械工业出版社）P370页表13-4可查得?=2.4。

由于此传动轴工作在弯扭组合交变应力状态下，因此在进行疲劳强度计算时疲劳强度条件可写成n???n?n?n??n?222?n?2。

2?maxM?max?WMy?MzW，?maxMx??3??3，W?，WP?。 ?3216WP?min???max，故弯矩循环系数r=-1,循环特征为对称循环；

?min??max，故扭矩循环系数r=0,循环特征为静应力。

所以n????1K?，n????1K?。

????max????a????m其中?a?0，?m??max，?????1??02?02。

参照《材料力学》（第二版）（机械工业出版社）P373页表13-5可选取???0.10。

在D截面右侧处：

，根据《材料力学》（第二版）?1?64mm，传动轴的材料为优质碳素结构钢（牌号45）（机械工业出版社）P369页表13-2可查得???0.78，???0.74。

MDy?1387.829Nm，MDz?352.646Nm，MDx?210.715Nm

?D,max?13222MDy?MDz?1387.829Nm2?352.646Nm2?56.981MPa 3W1??1?max?MDx16?210.715Nm?4.193MPa 3WP??1则n????1K??300MPa1.856.981MPa0.78?2.4?5.475

????maxn????1K??155MPa1.48x0?0.10?4.197MPa0.74?2.45.475?369.706?369.706

???n???

?a????mn?n??n??n?225.475?369.70622?5.475?n?2，安全。

在B截面右侧和E截面左侧处，?3?53mm，传动轴的材料为优质碳素结构钢（牌号45），根据《材料力学》（第二版）（机械工业出版社）P369页表13-2可查得???0.81，

???0.76。

在B截面右侧处：

MBy?547.580Nm，MBz?184.075Nm，MBx?108.222Nm

?B,max?1W3MBy?MBz?2232??33547.580Nm2?184.075Nm2?40.725MPa

?B,max?则n??MBx16?108.222Nm?3.815MPa 3WP3??3??1K??300MPa1.840.725MPa0.81?2.4?155MPa?7.956

???n??K??max??1????a????mn?n?n??n?221.48x0?0.10x3.185MPa0.76?2.47.956?406.331?406.331

n???

?7.956?406.33122?7.954?n?2，安全。

在E截面左侧处：

MEy?697.873Nm，MEz?436.932Nm，MEx?210.715Nm

?E,max?1W3MEy?MEz?2232??33697.873Nm2?436.932Nm2?58.044MPa

?E,max?MEx16?210.715Nm?7.427MPa WP3??33

则n????1K??300MPa1.858.044MPa0.81?2.4?155MPa?5.582

???n??K??max??1????a????mn?n?n??n?22?208.690

1.48x0?0.10x7.472MPa0.76?2.45.582?208.690?5.580?n?2，安全。

n???

?5.582?208.69022II．对传动轴阶梯轴进行疲劳强度计算

由于σb=650MPa，?1?72mm，?2?64mm，?3?58mm，?4?54mm,

?1?2?3（第二版）（机械工业???1.1，阶梯轴过渡圆弧r均为2mm，根据《材料力学》

?2?3?4出版社）P369页图13-9a，图13-9c，图13-9d，图13-9e可查得： 在P截面处

r?4r??2mm?0.0380，3?1.1，所以K?=1.73，K?=1.40；

54mm?4?2mm?0.0346，2?1.1，所以K?=1.76，K?=1.45；

58mm?3?2mm?0.0314，1?1.1，所以K?=1.80，K?=1.48；

64mm?2在Q截面处

?3r?在U截面处

?2r?在V截面处

?3r??2mm?0.0346，1?1.21，所以K?=2.25，K?=1.70；

58mm?3?2mm?0.0380，3?1.1，所以K?=1.73，K?=1.40；

54mm?4在W截面处

?4?在P截面处：

?4?48mm，传动轴的材料为优质碳素结构钢（牌号45），根据《材料力学》（第二版）

（机械工业出版社）P369页表13-2可查???0.84，???0.78。

MPy?273.790Nm，MPz?92.038Nm，MPx?0Nm

?P,max?1W4MPy?MPz?2232??43273.790Nm2?92.038Nm2?27.103MPa

?P,max?MPx16?0Nm?0MPa 3WP4??4则n????1K?????n??K??max300MPa?12.899

1.7327.103MPa0.84?2.4?155MPa???

1.40x0?0.10x00.78?2.412.899?(??)12.899?(??)22??1????a????mn?n?n??n?22n???

??12.899?n?2，安全。

在Q截面处：

，根据《材料力学》（第二版）?3?53mm，传动轴的材料为优质碳素结构钢（牌号45）（机械工业出版社）P369页表13-2可查得???0.81，???0.76。

MQy?757.642Nm，MQz?226.218Nm，MQx?108.222Nm

?Q,max?1W3MQy?MQz?2232??33757.642Nm2?226.218Nm2?55.741MPa

?Q,max?则n??MQxWP3K??16??33108.222Nm?3.815MPa 300MPa?5.945

1.7655.741MPa0.81?2.4??406.331

1.45x0?0.10x3.815MPa0.76?2.45.945?406.3315.945?406.33122??1????max?n????1K?155MPa???n????a????mn?n??n??n?22?5.944?n?2，安全。

在U截面处：

?2?58mm，传动轴的材料为优质碳素结构钢（牌号45），根据《材料力学》（第二版）

（机械工业出版社）P369页表13-2可查得???0.81，???0.76。

MUy?1177.767Nm，MUz?310.504Nm，MUx?108.222Nm

?U,max?1W2MUy?MUz?2232??231177.767Nm2?310.504Nm2?114.286MPa

?U,max?MUx16?108.222Nm?5.077MPa 3WP2??2则n????1K?????n??K??max300MPa?2.835

1.80114.286MPa0.81?2.4?155MPa1.48x0?0.10x5.0770.76?2.42.835?305.2822.835?305.28222??1????a????mn?n?n??n?22?305.282

n???

??2.835?n?2，安全。

在V截面处：

，根据《材料力学》（第二版）?3?53mm，传动轴的材料为优质碳素结构钢（牌号45）（机械工业出版社）P369页表13-2可查得???0.81，???0.76。

MVy?1042.851Nm，MVz?394.789Nm，MVx?210.715Nm

?V,max?1W3MVy?MVz?2232??331042.851Nm2?394.789Nm2?78.609MPa

?V,max?MVx16?210.715Nm?7.427MPa WP3??33则n????1K?????n??K??max300MPa?3.297

2.2578.609MPa0.81?2.4?155MPa1.70x0?0.10x7.4270.76?2.43.297?208.6903.297?208.69022??1????a????mn?n?n??n?22?208.690

n?????3.296?n?2，安全。

在W截面处：

，根据《材料力学》（第二版）?4?48mm，传动轴的材料为优质碳素结构钢（牌号45）（机械工业出版社）P369页表13-2可查得???0.84，???0.78。

MWy?348.937Nm，MWz?218.466Nm，MWx?0Nm

?W,max?1W4MWy?MWz?2232??43348.937Nm2?218.466Nm2?38.629MPa?W,max?MWx16?0Nm?0MPa 3WP4??4则n????1K?????n??K??max300MPa?9.050

1.7338.629MPa0.84?2.4?155MPa???

1.400MPa0MPa?0.100.78?2.4229.050?(??)9.050?(??)22??1????a????mn?n?n??n?22n?????9.050?n?2，安全。

现将各校核截面参数整理后列表如下： 校 核 点 D B E P Q U V W 参 数 初始应力 集中系数 尺寸系数 表面质量 系数? 敏感 系数 直径 ?（mm） 72 58 58 54 58 64 58 54 K? 1.80 1.80 1.80 1.73 1.76 1.80 2.25 1.73 K? 1.48 1.48 1.48 1.40 1.45 1.48 1.70 1.40 ?? 0.78 0.81 0.81 0.84 0.81 0.81 0.81 0.84 ?? 0.74 0.76 0.76 0.78 0.76 0.76 0.76 0.78 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 ?? 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1

各校核截面计算结果如下：

校核点 D B E P Q U V W

?max（MPa） ?max（MPa）56.981 40.725 58.044 27.103 55.741 114.286 78.609 38.629 4.193 3.815 7.427 0 3.815 5.077 7.427 0 n? 5.475 7.956 5.528 12.899 5.945 2.835 3.297 9.050 n? 369.706 406.331 208.690 +? 406.331 305.282 208.690 +? n?? 5.475 7.954 5.580 12.899 5.994 2.835 3.296 9.050 综上所述，阶梯传动轴各个截面符合疲劳强度条件。

由于阶梯传动轴各个截面均符合疲劳强度条件，故本题不需要采取改进措施来改善

疲劳强度。

4—6对所取数据的理论根据作必要的说明

本次课程设计所取的数据均选取于参考文献（1）。

4—7设计所使用软件名称 制图：AUTO CAD 2004

数据计算：VISUAL C++ 6.0 编译程序 截图：红蜻蜓抓图精灵

附录：本题所编写的C程序

该程序的源程序如下所示。只要输入该题的任何一组数据，便可得到所求的答案。本次课程设计所得数据均来自该程序。

#include #define Pi 3.141593 float a,b,m;

float max(float a,float b) { if(a>b) m=a;

else m=b; return m; }

float Mw(float a,float b) {

return(sqrt(pow(a,2)+pow(b,2))); }

main() {

float d1,d2,d3,d4; float M,M1,F,F1,F2; float Fy1,Fz1,Fy2,Fz2;

float P,P1,n,D,D1,D2,G2,G1,a,Alpha;

float MxA,MxB,MxC,MxD,MxE,MxF,MxP,MxQ,MxU,MxV,MxW; float MyA,MyB,MyC,MyD,MyE,MyF,MyP,MyQ,MyU,MyV,MyW; float MzA,MzB,MzC,MzD,MzE,MzF,MzP,MzQ,MzU,MzV,MzW; float Pa=80*pow(10,6),E=200*pow(10,9);

float MmaxA,MmaxB,MmaxC,MmaxD,MmaxE,MmaxF,mmaxP,MmaxQ,MmaxU,MmaxV,MmaxW; float Iy,Iz,fy,fz,f,MM,T1,T2,T3,T4; double a1[10],b[10],c[10]; int num; char h;

double ko,kt,bat,d,Eyz,Ex,NO,i2; double omax,tmax,no,nt,w,n1;

double o1=300000000,t1=155000000,T=0.10; NO=0;

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

printf(\scanf(\

Alpha*=Pi/180; D/=1000; D1/=1000; D2/=1000; a/=1000;

M=9549*P/n; M1=9549*P1/n; F=2*(M-M1)/D2; F1=2*M1/D1; F2=2*M/D;

printf(\

Fy1=(4*F*cos(Alpha)+2*G1+6*F1+G2)/5; Fz1=(4*F*sin(Alpha)+3*F2)/5;

Fy2=(F*cos(Alpha)+3*G1+9*F1+4*G2)/5; Fz2=(4*F*sin(Alpha)+12*F2)/5;

printf(\

/*Draw picture*/ MxC=F*D2/2; MxE=F2*D/2;

MyB=Fy1*a;

MyD=(Fy1-F*cos(Alpha))*3*a+Fy1*a; MyE=Fy2*a;

MzB=Fz1*a; MzE=Fz2*a;

printf(\3fNm,MzE=%0.3fNm\\n\

/*Calculate diameter*/

/*Calculate according to the axis which has the same diameter*/ MyD=MyD;

MzD=MzB/3+MzE*2/3; MxD=MxE;

MmaxD=sqrt(pow(MyD,2)+pow(MzD,2)+pow(MxD,2)); MyE=MyE; MzE=MzE; MxE=MxE;

MmaxE=sqrt(pow(MyE,2)+pow(MzE,2)+pow(MxE,2));

printf(\

printf(\printf(\printf(\d1=pow(32*max(MmaxD,MmaxE)/Pi/Pa,1.0/3); d2=d1/1.1; d3=d2/1.1; d4=d3/1.1;

printf(\d4*1000);

/*Calculate according to d2*/ MyU=MyB/4+MyD*3/4; MzU=MzB/2+MzE/2; MxU=MxC;

MmaxU=sqrt(pow(MyU,2)+pow(MzU,2)+pow(MxU,2)); printf(\

printf(\printf(\if(pow(32*MmaxU/Pi/Pa,1.0/3)>d2) {

printf(\d2=pow(32*MmaxU/Pi/Pa,1.0/3); d1=d2*1.1; d3=d2/1.1; d4=d3/1.1;

printf(\0,d4*1000); }

else printf(\

/*Calculate according to d3*/ MyQ=MyB*3/4+MyD/4; MzQ=MzB*5/6+MzE/6; MxQ=MxC;

MmaxQ=sqrt(pow(MyQ,2)+pow(MzQ,2)+pow(MxQ,2)); MyV=MyD/2+MyE/2; MzV=MzB/6+MzE*5/6; MxV=MxE;

MmaxV=sqrt(pow(MyV,2)+pow(MzV,2)+pow(MxV,2)); printf(\

printf(\printf(\printf(\

printf(\

if(pow(32*max(MmaxQ,max(MmaxV,MmaxE))/Pi/Pa,1.0/3)>d3) {

printf(\

d3=pow(32*max(MmaxQ,max(MmaxV,MmaxE))/Pi/Pa,1.0/3); d2=d3*1.1; d1=d2*1.1; d4=d3/1.1;

printf(\0,d4*1000); }

else printf(\printf(\

d1=%0.3fmm,d2=%0.3fmm,d3=%0.3fmm,d4=%0.3fmm\\n\;

/*Calculate deflection*/ MM=a;

Iz=Iy=Pi*pow(d1,4)/64;

fy=(1.0/2.0*MyB*a*2.0/3.0*MM+MyB*2.0*a*3.0/4.0*MM+1.0/2.0*(MyD-MyB)*2.0*a*2.0/3.0*MM+MyE*a*3.0/8.0*MM+1.0/2.0*(MyD-MyE)*a*5.0/12.0*MM+1.0/2.0*MyE*a*1.0/6.0*MM)/E/Iz;

fz=(1.0/2.0*MzB*2.0/3.0*MM+MzB*3.0*a*5.0/8.0*MM+1.0/2.0*(MzE-MzB)*3.0*a*1.0/2.0*MM+1.0/2.0*MzE*a*1.0/6.0*MM)/E/Iy; f=sqrt(pow(fy,2)+pow(fz,2));

printf(\printf(\printf(\

T1=d1;T2=d2;T3=d3;T4=d4;

a1[0]=MyP=MyB/2; /*赋值语句*/ b[0]=MzP=MzB/2; c[0]=MxP=0; a1[1]=MyB; b[1]=MzB; c[1]=MxB=MxC; a1[2]=MyQ; b[2]=MzQ; c[2]=MxQ;

a1[3]=MyD; b[3]=MzD; c[3]=MxD; a1[4]=MyV; b[4]=MzV; c[4]=MxV; a1[5]=MyE; b[5]=MzE; c[5]=MxE;

a1[6]=MyW=MyE/2; b[6]=MzW=MzE/2; c[6]=MxW=0; a1[7]=MyU; b[7]=MzU; c[7]=MxU;

for(i2=0;i2<28;i2++) {

printf(\请你输入要测定疲劳强度的截面\\n\ scanf(\ printf(\截面\\n\if(h=='P'||h=='p') {

num=0;

ko=1.73,kt=1.40,bat=2.4,Eyz=0.84,Ex=0.74; /* d=T4; }

else if(h=='B'||h=='b') {

num=1;

ko=1.80,kt=1.48,bat=2.4,Eyz=0.81,Ex=0.76; d=T3; }

else if(h=='Q'||h=='q') {

num=2;

ko=1.76,kt=1.45,bat=2.4,Eyz=0.81,Ex=0.76; d=T3; }

else if(h=='D'||h=='d') {

num=3;

ko=1.80,kt=1.48,bat=2.4,Eyz=0.78,Ex=0.74;

赋值*/

d=T1; }

else if(h=='V'||h=='v') {

num=4;

ko=2.25,kt=1.70,bat=2.4,Eyz=0.81,Ex=0.76; d=T3; }

else if(h=='E'||h=='e') {

num=5;

ko=1.80,kt=1.48,bat=2.4,Eyz=0.81,Ex=0.76; d=T3; }

else if(h=='W'||h=='w') {

num=6;

ko=1.73,kt=1.40,bat=2.4,Eyz=0.84,Ex=0.78; d=T4; }

else if(h=='U'||h=='u') {

num=7;

ko=1.80,kt=1.48,bat=2.4,Eyz=0.81,Ex=0.76; d=T4; } else {

continue; }

printf(\

w=Pi*d*d*d; /*检测疲劳强度过程*/ omax=32*Mw(a1[num],b[num])/w; tmax=16*c[num]/w; if(tmax

no=o1/ko/omax*bat*Eyz; nt=t1/(T*tmax); n1=no*nt/Mw(no,nt); if(tmax!=NO) {

if(omax!=NO)

{

printf(\ if(n1>=2)

printf(\在%c截面处安全\\n\ else

printf(\在%c截面处不安全\\n\ }

else {

n1=nt;

printf(\no=无限大 \\n\ if(n1>=2)

printf(\在%c截面处安全\\n\ else

printf(\在%c截面处不安全\\n\ }

} else {

n1=no;

printf(\，tmax=%0.3f，no=%0.3f，n=%0.3f, nt=无限大\\n\

if(n1>=2)

printf(\在%c截面处安全\\n\ else

printf(\在%c截面处不安全\\n\ }}}

运行结果如下:

参考文献

（1）聂玉琴，孟广伟主编. 材料力学（第二版）. 北京：机械公业出版社，2009. （2）谭浩强主编. C程序设计（第四版）. 北京：清华大学出版社，2010.

材料力学课程设计的体会和收获：

通过这次课程设计，我深切体会到了材料力学在工程实际中的作用，并且加深了对基本概念的理解，巩固了所学的知识。

在课程设计中，遇到了很多困难，但是在克服困难后我体会到了成功的喜悦，在其他方面我也获得了很多，从整体上掌握了基本理论和现代的计算方法，提高了分析问题，解决问

题的能力，把以前所学的知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等）综合运用，初步掌握工程中的设计思想和设计方法，提高了实际的工作能力。 感谢各位老师在日常学习中给予我的帮助!

本文来源：https://www.bwwdw.com/article/fl88.html