材料力学课程设计-车床主轴设计

材料力学课程设计
设计题目：车床主轴设计
数据序号: Ⅰ12 班级： 10级工机2班
学号: 20106994
姓名：冯华钊
目录
一、材料力学课程设计的目的
二、材料力学课程设计的任务和要求
三、设计题目
四、对主轴静定情况校核
1.根据第三强度理论校核
2.根据刚度进行校核
3.疲劳强度校核
五、对主轴超静定情况校核
1.根据第三强度理论校核
2.根据刚度进行校核
3.疲劳强度校核
六、循环计算程序
七、课程设计总结
一、设计目的
材料力学课程设计的目的是在于系统的学习材料力学之后，能结合工程中的实际问题，运用材料力学设计的基本原理和计算方法，独立计算工程中的典型零部件，已达到综合运用材料力学的知识解决工程实际问题的能力。

同时，可以使我们将材料力学的理论和现代的计算方法及手段融为一体。

即从整体上掌握了基本理论和现代的计算方法，又提高了分析问题，解决问题的能力；即把以前学到的知识综合的运用，又为以后的学习打下了基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。

1.使我们的材料力学知识系统化，完整化。

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

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

4.综合运用以前所学的各门课程知识，是相关学科知识有机的联系起来。

5.初步了解和掌握工程实践中的设计思想和方法，为以后打下基础。

二、设计的任务和要求
1.画出受力分析计算简图和内力图
2.列出理论依据和导出的计算公式
3.独立编制计算机程序，通过计算机给出计算结果
4.完成设计说明书。

三、设计题目
车床主轴设计---
某车床主轴尺寸及受力情况如图1所示。

在A、B、C三个支座的中间支座B处，轴承与轴承座之间有间隙δ，正常工作时，B处轴承不起支撑作用，此时轴处于A、C两支座下的静定状态。

当B截面处弯曲变形大于间隙δ时，轴处于A、B、C三支座下的静不定状态。

轴截面E处装有斜齿轮，其法向压力角为α，螺旋角为β，工作处的切削力有Fx、Fy、Fz（在进行强度、刚度计算时，可以不计轴向力Fx的影响，而以弯曲、扭转变形为主）。

轴的材料为优质碳素结构钢（45钢），表面磨削加工，氮化处理。

其他已知数据见表1。

1、试按静定梁（A、C支撑）的强度、刚度条件设计等截面空心圆轴外
径D(d/D值可见数据表2)，并计算这时轴上B截面处的实际位移。

2、在安装齿轮的E截面处有一铣刀加工的键槽，试校核此截面处的疲
劳强度。

规定的安全系数n=3（1-σ=420a MP ，1-τ=240a MP ）。

3、 对静不定情况（A 、B 、C 支撑），同时根据强度、刚度条件设计外径
D ，并用疲劳强度理论校核。

表1：
)(︒α )(︒β m /δ
MPa /][σ m f D /][ m f E /][ rad c /][θ
20
10
0.5410-⨯ 150 3.3410-⨯ 3.5410-⨯ 0.0028
注意：设计中不考虑轴的旋转静定要求和热变形的影响，并且将各轴承视为刚体，且不产生刚体位移，不考虑制造工艺和尺寸链等因素。

表2：（设计计算数据表Ⅰ12）
1
l /m
2l /m
3l
/m
A /m
B /m
R /m
θ (
)︒ n/
(r/min )
P
/kw
D
d
y H F /N
Z H F /N
12 0.1
6 0.48
0.15
0.12
0.16
0.12
45 500 5.4 0.7
4000 2500
图一：
一、 对主轴静定情况校核
由公式可知Me=9549⨯
min
/}{}{r kw n p = 9549 5.4
500⨯=103.13N m ∙
∴F t =
103.13
0.12
=859.41N 由斜齿轮受力分析得： F r =
t tan cos F αβ=859.410.364
0.985
⨯=317.59N 则有：F y E =F t sin θ-F r cos θ=383.12N F Z E =F t cos θ+F r sin θ=832.26N F y b= F Hy b=4000⨯0.16=640N ∙m
F Z b= F Hz b=2500⨯0.16=400N ∙m
由图1受力分析求支座反力F Ay 、F Az 、F Cy 、F Cz ：
)(F M Cz
∑= F Ay (L 1+L 2)+F Ey a-640- F y L 3=0
∴ F Ay =2091.30N
∑)(F M Az
= F Cy (L 1+L 2)+ F Ey (L 1+L 2-a)+640+ F y ( L 1+L 2+L 3)=0
∴ F Cy =-6241.32N
∑)(F M Cy
= F Az (L 1+L 2)+ F Z E a+400+ F Z L 3=0
∴F Az =-1157.53N
∑)(F M
Ay
= F Cz (L 1+L 2)+ F Ez (L 1+L 2-a)-400- F z ( L 1+L 2+L 3)=0
∴ F Cz =2681.14N
根据已知分别作出Y 、Z 方向的剪力图与弯矩图，如下图所示：
由剪力图及弯矩图可知c 点为危险点且： Mc=227041408+=1574.19N ∙m Me=81.17N ∙m 1.根据第三强度理论校核： W
Me r 223
Mc +=
σ][σ≤ 且 )1(3243
απ-=D W 代入数据解得： D 1≥39.5210-⨯m 2.由刚度对轴进行校核： 利用图乘法∑
=-
-=∆n
i ci
i EI
M 1ϖ对各点进行刚度校核：
1)根据D 点刚度计算轴径，在D 点分别沿y 、z 轴加一单位力有扭矩图如下图
=Dy f 149.078.25714.021128.03222.115055.021
1
⨯⨯⨯+⨯⨯⨯⨯⎢⎣⎡EI
()3216.070416.02122.115014.016.0128.021
⨯⨯⨯⨯+⨯⨯++
EI 89
.672116.016.0704=
⎥⎦⎤⨯⨯⨯+
()]EI
EI f Dz 99.3616.03235216.02108.035216.069.03214.055.064.63670414.02169.062.016.064.63614.0128.03264.63655.02
11-=⨯⨯⨯⨯+⨯⨯+⨯+⨯
-⨯⨯+⨯⨯⨯+⨯⨯⨯⨯⎢⎣⎡-
=
=+2
2Dz Dy D f f f EI
63.78m f D 4103.3][-⨯=≤ E=210Pa 910⨯ I=44403.0)1(64
D D =-απ
4
9410
3.303.01021063.78-⨯⨯⨯⨯≥∴D 2
21084.7-⨯≥∴D m 2) 根据E 点刚度计算轴径，在E 点分别沿y 、Z 轴加一单位力有扭矩图如下图
()EI
EI f Ey 63
.41]22.1150140814.021140.03122.115014.0140.021[1=
-⨯⨯⨯⨯+⨯⨯⨯=
()EI
EI f Ez 80.22]14.03164.63670414.02114.064.63614.021
64.6363214.055.021[1-=⨯⨯-⨯⨯+⨯⨯⨯+⨯⨯⨯⨯-
= ][46
.472
2E Ez Ey E f EI
f f f ≤=
+= 即：
44
9105.303.01021046
.47-⨯≤⨯⨯D
解得：D321081.6-⨯≥m
3）根据C 点刚度计算直径，在C 点处加一单位力偶得如下图所示弯矩图：
EI
EI cy 39
.495]114.0)22.11501408(21155.022.115021[1=
⨯⨯++⨯⨯⨯=
θ
EI
EI Cz 92
.268]114.0)70464.636(21164.63655.021[1-
=⨯⨯+⨯+⨯⨯⨯-
=θ 0028
.0][67.5632
2=≤=+=C cz cy C EI
θθθθ 即：
0028.003.01021067
.5634
9≤⨯⨯D
解得：≥4D 21052.7-⨯m
综上所述：D=max[D 1、D 2、D 3、D 4]=7.84210-⨯m 当D= 6.97210-⨯m 时，计算B 点的实际位移：（应用图乘法）
]128.052
.007.02122.115014.052.038
.031
14.0128.038.022.115019.055.038.052
.033.0128.038.022.115055.017.055.017.017.03222.1150128.021[1⨯⨯⨯⨯+⨯+⨯⨯⨯⨯⨯+⨯
⨯⨯⨯+⨯⨯⨯⨯⨯=EI f By =EI
06
.25
()128
.052
.014.03121
64.63670414.0128
.052
.007.064.63614.0128.052.038.03114.064.63655.038.019.052
.019.014.0128.064.63655.017.038.055.017.0128.03264.63617.021[1⨯⨯⨯⨯-⨯+⨯⨯⨯+⨯⨯+⨯⨯⨯++⨯⨯⨯⨯+⨯⨯⨯⨯⨯-
=EI f Bz =-EI 50.17 m EI f f f Bz By B 44
29221028.1)
1084.7(03.01021057
.3057.30--⨯=⨯⨯⨯⨯==
+=
3.疲劳强度校核:
若不计键槽对抗弯截面系数的影响，则危险截面处抗弯截面系数：
3643
1032.32)1(32
--⨯=-=
m D απϖ
由弯矩M 不变可知该循环为对称循环，则有： MPa Pa W M 68.4010
32.3264.63622.115062
2min
max =⨯+==-=-σσ
MPa Pa W M P X 26.110
65.6417
.816
max =⨯==
-τ 查表确定铣加工的键槽危险截面处疲劳强度的影响系数：
60.1=σK 88.1=τK 75.0=σε 73.0=τε
8.1=β
则：71.868.408
.175.060
.1420max
1
=⨯⨯=
=
-MPa
MPa
K n σβ
εσσσ
σ
13.13326.18
.173.088
.1240max 1
=⨯⨯=
=
-MPa
MPa
K n τβ
ετττ
τ
369.82
2
>=+=τ
στστσn n n n n 故E 处满足疲劳强度要求。

二、 对超静定情况进行校核
由m f m B 441028.1105.0--⨯=<⨯=δ，故此轴为超静定，且为一次静不定。

由变形协调条件可知： δ-=B F f f B 。

分别沿y 、z 轴加一单位力并作By F 、Bz F 、及单位力的弯矩图有：
By By By By FBy F EI
F F EI f f 00377.0]128.032128.052.021128.032128.017.021[1=⨯⨯⨯⨯+⨯⨯⨯⨯=
-=δ 又EI f By
06
.25=
; 代入上式有：N EI F By 15.349500377
.006.25=-=δ
Bz Bz Bz Bz FBz F EI
F F EI f f 00377.0]128.032128.052.021128.032128.017.021[1=⨯⨯⨯⨯+⨯⨯⨯⨯=-=δ
又EI
f Bz 50.17=;代入上式有： N EI F Bz 18.148700377
.050.17=-=δ 从而求A 、C 点的支反力有：
1223()()6400Cy Ey Ay By y M F F a F L L F L F L ∑=+++--=
N F Ay 73.542-=∴
12112123()()()640()0
Ay Cy By Ey y M F F L L F L F L L a F L L L =++++-++++=∑ N F Cy 44.7102-=∴
1223()()4000Cz Az Bz Ez Z M F F L L F L F a F L =+++++=
N F Az 30.2278-=∴
12112123()()()400()0Az Cz Bz Ez Z M F F L L F L F L L a F L L L =++++---++=∑
N F Cz 73.2314=∴
做剪力图Qy F 、Qz F 如下所示：
做弯矩图y M 、Mz 如下图所示：
由上图有：m
N M M Mc m N M M Me cz cy ez ey ∙=+=∙=+=19.157433.123822
22
故C 点为危险点
1）.第三强度理论校核有： ][1
223σϖσ≤+=e c r M M 且 )1(3243απ
ϖ-=D
代入数据解得：m D 211084.5-⨯≥
2）.由刚度对轴进行校核： 利用图乘法∑=--=∆n i ci i EI M 1ϖ对各点进行刚度校核：
1. 根据D 点的刚度对主轴进行校核，分别沿Y 、Z 轴加一单位力得到如
下图所示弯矩图：
()()EI
EI f Dy 84.5316.032704140816.02116.02
170416.016.069.052.014.03266.1029140814.02116.069.062.066.102914.03
266.102935.02
155.018.0128.026.9203.021128.055.017.03226.9217.021[1=⨯⨯-⨯⨯+⨯⨯⨯+⨯+⨯⨯-⨯⨯+⨯⨯⨯+⨯⨯⨯+⨯⨯⨯⨯-⨯⨯⨯⨯⨯-=
()]EI
EI f Dz 28.4016.03216.03522116.02116.035216.069
.055.014.03214.094.6877042116.069.062.014.094.68738.0)31.38794.687(2
169.036.016.031.38738.016.069.017.03231.38717.021[1-=⨯⨯⨯⨯+⨯⨯⨯+⨯+⨯⨯⨯-⨯+⨯⨯⨯+⨯⨯-+⨯⨯⨯+⨯⨯⨯⨯⨯-= 则：m f EI
f f f D Dz Dy D 422103.3][24.67-⨯=≤=+= 解得：m D 211054.7-⨯≥ 2. 根据E 点的刚度对轴校核：有静定情况可知23D D ≤
3. 根据C 点的刚度对轴校核：在C 点分别加一绕Y 、Z 方向的单位偶有
扭矩图如下：
166.102935.02
1126.9203.021126.9217.021[1⨯⨯⨯+⨯⨯⨯-⨯⨯⨯-=
EI Cy θ EI 60.341]114.0)140866.1029(21-=⨯⨯++
]114.0)70494.687(2138.0)94.68731.387(21131.38717.021[1⨯⨯+⨯+⨯+⨯+⨯⨯⨯-
=EI Cz θ
EI
17.335-= 则：0028.0][57.47822=≤=
+=C Cz Cy C EI θθθθ 解得：m D 241022.7-⨯≥ 综上所述：m D D D D D 243211054.7],,,m ax [-⨯==
3）疲劳强度校核：
查机械手册得到：60.1=σK 88.1=τK 75.0=σε 73.0=τε 8.1=β则： 3643
1075.28)1(32m D W -⨯=-=απ；
MPa W M 07.4010
75.2833.12386min max =⨯==-=-σσ MPa W M P X 41.1max ==
τ 84.807.408.175.060.1420max 1
=⨯⨯==-σβ
εσσσσK n ； 97.11841.18
.173.088.1240max 1
=⨯⨯==-τβετττσK n ； 382.822>=+=τστ
σn n n n n ；故满足强度条件。

三、 循环计算程序
#include<stdio.h>
#include<math.h>
#define pi 3.141592654
#define ip 0.017453292
floatL1,L2,L3,a,b,A0,n,P,i,Fy,Fz,Fby=0,Fbz=0,Fcy,Fcz,Fey,Fez,Mby,Mbz, Mcy,Mcz,Mey,Mez,,Md,Mc,www,xxx=150,eee=0.21,ddd,fby0,fbz0,fby1,fbz1, Fffe=0.00035,fffd=0.00033,aaac=0.0028,aMb,aMc,aMe,aaaac,fffb=0.00005, Fffbb,big,SM,SN,B1,SD,FFF;
void zaihe()
{long double Ft,Fy,An=20.0,Bn=10.0;
Ft=30*P*1000/(n*pi*R);
Fy=Ft*tan(An*ip)/(cos(Bn*ip);
Fey=Ft*sin(A0*ip)-Fy*cos(A0*ip;
Fez=Ft*cos(A0*ip)+Fy*sin(A0*ip);
Me=Ft*R;Md=-Me;Fdy=Fy;Fdz=-Fz;
Mdy=Fdy*b;Mdz=Fdz*b;}
void waili()
{Fay=(Fdz*L3+Mdz-Fez*a-Fbz*L2)/(L1+L2);Fcy=(-Fdy*(L1+L2+L3)-Mdy-Fey(L 1+L2-a)/(L1+L2);Faz=(Fdy*L3+Mdy-Fey*a-Fby*L2)/(L1+L2);Fcz=(-Fdz*(L1+L 2+L3)-Mdz-Fez(L1+L2-a) /(L1+L2);
Mby=Fay*L1;Mbz=Faz*L1;
Mey=Fay*(L1+L2-a)+Fby*(L2-a);Mez=Faz*(L1+L2-a)+Fbz*(L2-a);
Mcy=Fay*(L1+L2)+Fby*L2+Fey*a;Mcz=Faz*(L1+L2)+Fbz*L2+Fez*a;} void qiangdu()
{long double Xmax,Mmax,Xmax1,Xmax3,Xmax2=0.0;
Mc=sqrt(Mcy*Mcy+Mcz*Mcz);Xmax=Mc/www;Mmax=Me/(2*www);
Xmax3=-Xmax/2-sqrt((Xmax/2)*(Xmax/2)+Mmax*Mmax);
Xmax1=-Xmax/2+sqrt((Xmax/2)*(Xmax/2)+Mmax*Mmax);
ddd=sqrt(((Xmax1-Xmax2)*(Xmax1-Xmax2)+(Xmax2-Xmax3)*(Xmax2-Xmax3)+(Xm ax3-Xmax1)*(Xmax3-Xmax1))/2)/xxx;ddd=pow(ddd,0.333333333333333333333) ;ddd=ddd/100;ddd=(float)ddd;}
voidnaodu(){longdoublefffd1,fffdy,fffdz;fffdy=L1*Mby*L1*aMc/(L1+L2)/ 0.147/eee;fffdy+=(L2-a)*(Mby*(2*aMc*L1/(L1+L2)+aMc*(L1+L2-a)/(L1+L2)) +Mey*(2*aMc*(L1+L2-a)/(L1+L2)+aMc*L1/(L1+L2)))/(0.294*eee);fffdy+=a*( Mey*(2*aMc+(L1+L2-a)*aMc/(L1+L2)+aMc)+Mcy*(2*aMc+(L1+L2-a)*aMc/( L1+L 2)))/(0.294*eee);fffdy+=L3*Mac*(Mdy-2*Mcy)/(0.294*eee);fffdz=L1*Mbz*L 1*aMc/(L1+L2)/0.147/eee;fffdz+=(L2-a)*(Mbz*(2*aMc*L1/(L1+L2)+aMc*(L1+ L2-a)/(L1+L2))+Mez*(2*aMc*(L1+L2-a)/(L1+L2)+aMc*L1/(L1+L2)))/(0.294*e ee);
fffdz+=a*(Mez*(2*aMc+(L1+L2-a)*aMc/(L1+L2)+aMc)+Mez*(2*aMc+(L1+L2-a)* aMc/(L1+L2)))/(0.294*eee);
fffdz+=L3*Mac*(Mdz+2*Mcz)/(0.294*eee);
fffd1=pi*sqrt(fffdy*fffdy+fffdz*fffdz)/www/32/fffd;
ddd=pow(fffd1,0.25);ddd=ddd/1000;ddd=(float)ddd;}
void naodue()
{long double fffe1,fffey,fffez;
fffey=L1*Mby*L1*aMe/(L1+L2-a)/0.147/eee;
fffey+=(L2-a)*(Mby*(2*aMe*L1/(L1+L2-a)+aMe)+Mey*(2*aMe*aMe*L1/(L1+L2-a)))/ 0.294/eee; fffey+=a*aMe*(Mcy+2*Mey)/(0.294*eee);
fffez=L1*Mbz*L1*aMe/(L1+L2-a)/0.147/eee;
fffdz+=(L2-a)*(Mbz*(2*aMe*L1/(L1+L2-a)+aMe)+Mez*(2*aMe*L1/(L1+L2-a))) /(0.294*eee);
fffdz+=a*aMe*(Mcz-2*Mez)/(0.294*eee);
fffe1=pi*sqrt(fffey*fffey+fffez*fffez)/www/32/eee;
ddd=pow(fffe1,0.25);ddd=ddd/1000;ddd=(float)ddd;
void zhuanjiaoc()
{long double aaac1,aaacy,aaacz;
aaacy=L1*Mby*L1*aaac/(L1+L2)/0.147/eee;
aaacy+=(L2-a)*(Mby*(2*aaaac*L1/(L1+L2)+aaaac*(L1+L2-a)/(L1+L2))+Mey*( 2*aaaac*(L1+L2-a)/(L1+L2)+aaaac*L1/(L1+L2)))/(0.294*eee);
aaacy+=a*(Mey*(2*(L1+L2-a)*aaaac/(L1+L2)+aaaac)+Mcy*(2*aaaac+(L1+L2-a )/(L1+L2) ))/(0.294*eee);
aaacz=L1*Mbz*L1*aaac/(L1+L2)/0.147/eee;
aaacz+=(L2-a)*(Mbz*(2*aaaac*L1/(L1+L2)+aaaac*(L1+L2-a)/(L1+L2))+Mez*( 2*aaaac*(L1+L2-a)/(L1+L2)+aaaac*L1/(L1+L2)))/(0.294*eee);
aaacz+=a*(Mez*(2*(L1+L2-a)*aaaac/(L1+L2)+aaaac)+Mcz*(2*aaaac+(L1+L2-a )/( L1+L2) ))/(0.294*eee);
aaac1=pi*sqrt(aaacy*aaacy+aaacz*aaacz)/www/32/aaac;
ddd=pow(aaac1,0.25);ddd=ddd/1000;ddd=(float)ddd;}
void fanlib()
{long double fffb1,fffby,fffbz;
fffby=L1*Mby*aMb/0.147/eee;
fffby+=(L1+L2-a)*(Mby*(2*aMb+(L2-a)/L2*aMb)+Mey*(2*(L2-a)/L2*aMb+aMb) )/0.294/eee;
fffby+=a*(L2-a)*aMb*(Mcz+2*Mey)/(0.294*eee);
fffbz=L1*Mbz*aMb/0.147/eee;
fffbz+=(L1+L2-a)*(Mbz*(2*aMb+(L2-a)/L2*aMb)+Mez*(2*(L2-a)/L2*aMb+aMb) )/ (0.294*eee);
fffbz+=a*(L2-a)/L2*aMb*(Mcz+2*Mez)/(0.294*eee);
fby0=fffby;fbz0=fffbz;
fffb1=pi*sqrt(fffby*fffby+fffbz*fffbz)*pi/www/big/big/big/big/640e11; fffbb=(float)fffb1;
printf("B点的实际位移：S=%fm\n",fffbb)
void waili()
{Fby=fffb*SD*SD*SD*SD*(1-i*i*i*i)/FFF*10e12;
Fbz=Fby*fbz0/fby0;}
void xunhuan()
{double AAAA;fkoat big1,big2,big3,big4;
do
waili();
qiangdu();
big1=ddd;big=big1;
naodud();
big2=ddd;if(big<big1) big=big2;
naodue();
big3=ddd;if(big<big3) big=big3;
zhuanjiaoc();
big4=(float)ddd;if(big<big4) big=big4;
SM=SN;SN=BIG;SD=(SN+SM)/2;AAAA=fabs(SD-big);
fanlib();
while(AAAA>=0.0001);
printf("在超静定情况下的直径 :D=%fm\n",big);}
void main()
{float big1,big2,big3,big4;
printf(" 输入原始数据:\n");
scanf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f",&L1,&L2,&L3,&a,&b,&R, &A0,&n,&P,&i,&Fy,&Fz);
printf(" 所输入数据为:\nL1=%fm,L2=%fm,L3=%fm,a=%fm,b=%fm,R=%fm,A0=%f\n",L1,L2,L3,a,b,R,A 0);
printf("n=%fr/min,P=%fKW,i=%f,Fy=%fN,Fz=%fN\n",n,P,i,Fy,Fz);
aMc=1*L3;
aMe=(L1+L2-a)*a/(L1+L2)*1;
aaaac=1;
aMb=1*L1*12/(L1+L2);
www=pi*(1-i*i*i)/32;
zaihe();
waili();
printf("辅助计算外力（静定）:\nFay=%fN,Faz=%fN,Fby=%fN,Fbz=%fN,Fcy=%fN,Fcz=%fN,",Fay,Faz,Fby, Fbz,Fcy,Fcz);
printf("Fdy=%fN,Fdz=%fN,Fey=%fN,Fez=%fN,Mby=%fN.m,Mbz=%fN.m,Mcy=% fN.m,Mcz=%fN.m, ",Fdy,Fdz,Fey,Fez,Mby,Mbz,Mcy,Mcz);
printf("Mdy=%fN.m,Mdz=%fN.m,Mey=%fN.m,Mez=%fN.m,Me=%fN.m,Md=%fN.m \n",Mdy,Mdz,Mey,Mez,Me,Md);
qiangdu();
printf("D1=%fm\n",big1);big=big1;
naodud();
printf("D2=%fm\n",ddd);big2=(float)ddd;
if(big<big2) big=big2;
naodue();
printf("D3=%fm\n",ddd);big3=(float)ddd;
if(big<big3) big=big3;
zhuanjiaoc();
printf("D4=%fm\n",ddd);big4=(float)ddd;
if(big<big4) big=big4;
printf("D=%fm\n",big);
fanlib();
if(fffbb<=fffb) printf("可判定此结构为静定结构\n");
else
{printf(" 可判定此结构为超静定结构\n");
SM=0;SN=big;SD=(SN+SM)/2;
B1=1+fbz0*fbz0/(fby0*fby0);
B1=sqrt(B1);
FFF=(L1+L2)*aMb*aMb/(eee*0.147);
fanlib();
xunhuan();
printf(" \n");}
}
四、课程设计总结
本次课程设计涉及到了很多以前学习过的知识，包括材料力学、理论力学、AutoCAD、C语言编程等，透过本次课程设计，使我能更熟练的运用所学内容解决实际问题。

这次课程设计是我以前学习的一个阶段性总结，从中我看到了自己知识方面的不足以及学科综合的重要性。

此后，我会更加扎实的学习所学课程并广泛涉猎其他学科，在提高专业知识水平的基础上进一步提高自己的综合素质。