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

1、材料力学课程设计设计计算说明书设计题目：7.4车床主轴设计学号:41111317姓名：张振铭指导教师：魏媛老师目录一、课程设计目的 03二、课程设计任务和要求 03三、课程设计题目 04四、课程设计计算过程1. 对主轴静定情况校核 05A. 根据第三强度理论校核 07B. 根据刚度进行校核 07C. 疲劳强度校核 142. 对主轴超静定情况校核 14A. 根据第三强度理论校核 17B. 根据刚度进行校核 17C. 疲劳强度校核 23五、循环计算程序六、课程设计总结 2428一、课程设计目的材料力学课程设计的目的是在于系统的学习材料力学之后，能结合工程中的实际问题，运用材料力学设计的基本原理和计

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

3、选题力求结合专业实际，因而课程设计可以把材料力学的知识和专业4. 综合运用以前所学的各门课程知识（高等数学、工程图学、理论力学、算 法语言、计算机等），是相关学科知识有机的联系起来。5. 初步了解和掌握工程实践中的设计思想和设计方法。6. 为以后课程的学习打下基础。二、课程设计任务和要求参加设计者要系统复习材料力学课程的全部基本理论和方法， 独立分析、判 断设计题目的已知条件和所求问题， 画出受力分析计算简图和内力图，列出理论 依据并导出计算公式，独立编制计算程序，通过计算机给出计算结果，并完成设 计计算说明书。三、课程设计题目设计题目：车床主轴设计某车床主轴尺寸及受力情况如图1所示。在A、B

4、、C三个支座的中间支座B处，轴承与轴承座之间有间隙，正常工作时，B处轴承不起支撑作用，此时轴 处于A、C两支座下的静定状态。当B截面处弯曲变形大于间隙:时，轴处于A、 B、C三支座下的静不定状态。轴截面 E处装有斜齿轮，其法向压力角为：，螺 旋角为1，工作处的切削力有Fx、Fy、Fz （在进行强度、刚度计算时，可以不计轴向力Fx的影响，而以弯曲、扭转变形为主）。轴的材料为优质碳素结构钢（45 钢），表面磨削加工，氮化处理。其他已知数据见表1。1. 试按静定梁（A、C支撑）的强度、刚度条件设计等截面空心圆轴外径 D（d D 值可见数据表2），并计算这时轴上B截面处的实际位移。2. 在安装齿轮的E

5、截面处有一铣刀加工的键槽，试校核此截面处的疲劳强度。规定的安全系数n=3 （二=420MPa，=240MPa ）。3. 对静不定情况（A、B、C支撑），同时根据强度、刚度条件设计外径 D,并用疲劳强度理论校核。设计数据：表1:可m6/MPafD】/mfE】/m日c/rad2010_40.50150_43.3043.5X00.0028注意：设计中不考虑轴的旋转静定要求和热变形的影响，并且将各轴承视为刚体，且不产生刚体位移，不 考虑制造工艺和尺寸链等因素。表2:（设计计算数据表U 7）Ijm-m打/ma/mb mR m70.180.580.149()n/(r min )P/

6、kwd/DFHy/NFhz/N7454005.50.6545002200图一:a)L6-1 -h丄h1111b)四、课程设计计算过程1.对主轴静定情况校核由公式可知 Me =9549pkw =9549 5.5 =131.30 N m n r/ min400片二阻=131.299=937.85NR 0. 14由斜齿轮受力分析得:Fr = = 937.850.364 =346.58NCOS P则有：FEy 二 Ft sin，- Fr cos =418.10N0. 985Fez 二 Ft cos t Fr sin v =908.24NMy 二 Fzb 二卩出匕=2200 0.18=396 N mM

7、Dz = Fyb = FHy b=4500 0.18=810 N m由图1受力分析求支座反力FAy、FAz、FCy、Fez：二.M Cz (F )= FAyh I2 FEya M Dz FHy b =0M Az(F)二 FeyIlI2FEy IlI2 -a M DzFHyIlI2I3=01M Cy ( F )二 FAzh 12FEzaM DyFHzI -0:一MAy(F=FCzI12FEz I12 _a _ M Dy_FHz11I 2I3=0解上面的方程，则有：FAy=2129.98N, Faz =-1276.39N , Fey =-7048.08N , Fcz=2568.15N根据已知分别作

8、出Y、Z方向的剪力图与弯矩图，如下图所示：2130.0r2548 L，jl4500+，XMz1575.01192.8810.0-Fn2200. OpI x-368.-1276.4-714.5.-770.0.*3960由剪力图及弯矩图可知C点为危险点且:Me 二 J24027432 =1753.15 N mMe =131.30 N mA. 根据第三强度理论校核（忽略剪力）Me Mecr =t =W 2W-34WD）jMc2 + Me2 r i c r3（W代入数据解得：u 一5.2610JmB. 由刚度对轴进行校核：利用图乘法丄业对各点进行刚度校核：z EI1）根据D点刚度计算轴径，在D点分别沿

9、y、z轴加一单位力,有弯矩图如下：=333.98G| 0.134=.0892 二 1192. 80. 15=178.921 = 10. 561192. 81 2120. 151575. 0 - 1192.8 =18.67巫=20. 17 - 0. 1340. 134 =0.1583810. 000. 17=137.7Ml41-0. 1721575. 00 - 810. 00 =65.025MC5-20. 17=0.1133Dy580.7=5MC5 二 ei=Mci1 U1 M)1 亠2肌2 亠:3% 亠:4MUi 4 EI EI1 x 0. 56 x 714. 8=200 1442= 20.

10、134 =0.089c1 3M3714. 8 況 0. 15=107.22M2 Jc2 20. 1340. 17 =0.152I0.1570.0714.1423 。.134 0.134=0.158=396. 00 x 0. 17=67.320. 17=0.08510. 172Dz1-百1MC1： =2% v3MC3 : u4MC4.l 5MC5=44. 17ElJDy2-f DzElf J - 3. 310-mqE =210 10 Pa ID4(1 一 : 4) =0.04D46492. 0_ 210_10_0-Q4_3?3_10D2 _ 7. 610 -m2)根据E点刚度计算轴径，在E点分别

11、沿y、z轴加一单位力，有弯矩图如121 0. 56 1192. 8=333.98Mc10. 12=0.082312 = 1192.8 0. 15=178.92MC20. 12 =0.06113=20. 151575. 00 - 1192. 8 =28.665 W = 30. 12=0.0423fEy二晋诂2底旋二誉(770. 00 396. 00 )=31.79皿= 汉 0. 17=0.113My232 = 714.8 0. 15=107.22MC- 0.12=0.061 x 0. 15 汇(770.0 714.8=4.14M - 10. 12=0.043Dz_y 孙 MCi_i 4 El右

12、M 2%3MC3 =19. 7ET2Ey 2Ez- f E = 3. 510-4mE=210 109Pa I D4(1 - : 4) =0.04D464D4 -210 109a0043 3.5 10 -41 小3）根据C点刚度计算轴径，在C点处加一单位力偶，有弯矩图如下：XMz单位力偶：.i i0. 561192. 8=333.982 = 1192.80.15=178.92X：0. 151575. 00 - 1192. 8 =28.665 M d = M c2 二 M c3 = 1%y3I EIi MCi1541 565二百=1町亠： 2%亠：3皿3EP64My单位为偶X0. 56714. 8

13、=200.1442 = 714.80. 15=107.220. 15770. 0 - 714. 8 =4.14 M c1 = M c2 = M c3 = 1iMi右 M2玩3皿 311.504El鳶cz624. 76 空耳二 0. 0028radElE =210 109Pa ID4（1 - ： 4） =0.04D451D4 _ 7. 1810m综上所述：D = maxD-i， D2， D3， D4 J - 7. 6010-m624. 76210_10_Q-Q4_0. 0028当D=7.6010-m时，计算B点的实际位移：（应用图乘法）My0,14.10 18 2眄=厂0.曲再9234.506M

14、c-3X14=0-0930 182 二睛嗦80 56 - 0.18 =145921% = 20. 140. 04 =0.091 （ 0 18 叭二訂（。56 - 0. 18）沖192. 8 一 睛狀 1192. 8 |=153.7861MC3O.14。4Og。07314 = 1192. 80. 15=178.92MC40. 04=0.02吧冷.15575 - 119228.665M = 10. 04 =0.0133z亠： ：2皿2亠： ：3皿3亠： ：4皿4亠： ：5皿5=i Ai MCiElEl1MC131.50EI1 0 18 2= 2X0.16 X 714.8 =20-68叫飞714=.

15、932 二 0.8714. 80. 56 - 0.18 =87.310. 561MU = 0. 140. 04 =0.091(018、叭=訂(0. 56-0.18 H ；714.8 -睛14.8严.161。14 - 0.04O.。407314 = 714.80. 15=107.22MC40. 04=0.024 25 = 10. 15770 - 714. 8 =4.142MC = 10. 04 =0.0133r丄珂MCi1 矿丄 k丄 k丄 匸厂丄 二18. 71fBz 二-皿1223皿3445MC5i 珀 ElElEI964E =210 10 Pa ID4(1 - : 4) = 0.04D4B

16、y fBz36. 6436.64El210 109 0. 047.60 10-= 1.3110-4mC.疲劳强度校核: 若不计键槽对抗弯截面系数的影响，则危险截面处抗弯截面系数:32(.5410计由弯矩M不变可知该循环为对称循环，则有:二J minM 119282 714.82 P 39.28MPa3. 5410一5W = 2/y = 2W = 7.0810出-max131. 30Pa 二 1.85MPa7.0810查表确定铣加工的键槽危险截面处疲劳强度的影响系数:K；.60K 88 ；：一 = 0.75 ； =0.731.8J max420MPa肓% 39.28MPa二 9.02240MPa

17、=90. 67-max1.85MPa0.731.8n.C=8. 983故E处满足疲劳强度要求。由、.=0.510m ： f2.对超静定情况进行校核b =1.31 10，故此轴为超静定，且为一次静不定。由变形协调条件可知：fFB二fB -。分别沿y、Z轴加一单位力并作FBy、Fbz及单位力的弯矩图有:MyXFBy1 1百 0.18 0.14FBy2 130.14匸0.53屯0. 140. 00464ElFBy31. 50并且已知:ByEl代入上式有：FBy31.50 _ El = 5018. 61N0. 00464-6F巳2EI0. 00464EIFbz丄0. 180. 14FBz-0. 14-

18、0. 532z320. 141318. 71Bz并且已知：El代入上式有：FBz二 巴匸1旦 二1750. 69N0. 00464从而求A C点的支反力有：二.M CyF- F Ay111 2FBy I 2FEya _M Dz FHy l0二M AyF= FcyhI2FBy11FEy 11- a M Dz FHyl1 l2 l3 -0M CzF= FAzl12FBzl 2FEza M Dy FHzl = 0二 M Az F 二Fczh 丄FBzl1FEzl12 - a-MDy - F 吃 l 123- 0.FAy 二-1616. 31N Faz 二-2583. 24 N怎 -8320. 40N

19、 FCz = 2124.31N做剪力图FQy、FQz如下所示：-尸3820,40*-3402.3+，I1616.31450075.692200832朋2583,24做弯矩图My、Mz如下图所示:My X4百4.9霁78L35P770439反5由上图有：Me = M：MCZ 二 1753. 15N mc点为危险点M NA.第三强度理论校核有:=1 M；Mf叮且034 .EDI )代入数据解得：U _ 5. 2510 mB.由刚度对轴进行校核:n B M 利用图乘法归对各点进行刚度校核：y EI1）根据D点的刚度对主轴进行校核，分别沿 y、z轴加一单位力得到如下图所示弯矩图：Mn

20、6412A: 10.18290.94=268斫鹤 0.043=09-2290. 943402. 3290. 94=12.44MC2130. 064 一 0. 0430. 043 =0.050. 53 -3402.3290 940.15 -1001. 93=147.5280. 134-0. 0640. 064 =0.1111001.93 0. 15=150.29MC41 0. 13420. 17 =0.1520. 151575. 0 - 1001.93 =42.980. 17 - 0. 1340. 134 =0.1588100. 17=137.7= *0. 17=0.0850. 17=0.113

21、66.4410. 171575. 00 - 810. 00 =65.025 ML_ 1 ML三2 MC2L_ 3肌3 .三 4 仏DyElElElMy464.98,781.35770.0396.0-1 0.18 464.98=41.85 Mi0.043=0.0292 32 = 464. 980. 53 - 0. 15 =176.6910. 0430. 134 =0.08913 = 20. 53 - 0.15781.35 - 464. 98 =60.11023 o.134-。0430.043=004770. 000. 15=115.50MI41-0. 1340. 17 =0.1521X：3396

22、0. 17=67.32兀二-0.17=0.085DZ0. 17770. 00-396. 00 =31.79 皿7-0. 17=0.1133血1肌1 +可2% +可3皿3 +可4肌450. 19ElEl二-,f DyfD： fD = 3. 310*mElD10. 15781.35 - 770.00 =0.850. 17 - 0. 1340. 134=0.146E =210 109Pa ID4(14) = 0.04D464D4-83. 262101090. 043. 310一4D2-7. 4010，m2）根据E点的刚度对主轴进行校核，分别沿 y、z轴加一单位力得到如下图 所示弯矩图：1 = -0.

23、 18290.94=26.18瓦=20.04=0.027231290 94-= 2 x 340x 290.94=12.44 %=3沁.06 - 0.04)+ 0.04=0.0471(290 94轨飞10.53-0.15-侖1001.93 = 147.532。.06心.1014 = 1001.930. 15=150.29 皿40. 12=0.06二？ 0. 12=0.083i MciElElM 2%3ML34ML428.50ElMy464.98i781.35770.0站氐8冷 O15 1575.00 一 1001.93 =42.9812眄=一 x 0. 18 x 464. 98 =41.85M

24、= - x 0 04 =0.0272 32 二 464. 980. 53 - 0. 15 =176.69肌2 = 0. 040. 12 =0.08130. 53 - 0.15781.35 - 464. 98 =60.112 MU 工一 0. 12 - 0. 040. 04 =0.093313 = 770. 00 0. 15=115.50MC40. 12 =0.06124 0. 15781.35 - 770.00 =0.85 MU0. 12=0.0823r 二 四 MCi1 匚38.02fEz 八 一二-1皿12皿23%4皿45皿5 =-i _4 EIEIEIf E 二 f Ey f Ez 二

25、E|- fE】=3 5 10_m0.04D4E =210 109Pa I =d4CI4)=64D3 - 6. 341047. 52210 1090. 043. 5 1043）根据C点刚度计算轴径，在C点处加一单位力偶，有弯矩图如下：Mn290.94*Mz单位力偶XX：0. 18290.94=26.181=12.44聾 290. 940.53- 0.15 -叱3402.31001. 93=147.531001.93 汉 0. 15 =150.2910. 1521575. 00 - 1001. 93 =42.98M c1 = M c2 - M c3：cyJ i MCi1 .379 42百 1皿12

26、%3皿34皿4=5%MyJX464.9878135770.039氐6ElMz单位力偶X1 0. 18464. 98=41.8523402. 32 = 464. 980. 53 - 0. 15 =176.6913= 20.53 - 0.15781.35 - 464.98 =60.1114 =770.000. 15=115.5050. 15781.35 - 770.00 =0.852M C1 二 M C2 二 M c3二 M c4 二 M c5 = 1TCzJ i MCi2-1一 EI1 M?1亠12Mc2亠i 3M?3亠i 4M34亠：i 5MC5-395. 00El547. 71 空比=0.

27、0028radEl9E =210 10 Pa I444D4（1 - ： 4） =0.04D4 64D4 -547. 712101090. 040. 0028D4 _ 6. 9510综上所述：D = max, D2，D3，D二 7. 40 10，mC. 疲劳强度校核: 若不计键槽对抗弯截面系数的影响，则危险截面处抗弯截面系数:3（1 - ： 4）二 3. 2810f32由弯矩M不变可知该循环为对称循环，则有:二 maxminM1575.002770.002 pa = 53.45MPa_53. 28106.5610maxMXWP131 30Pa = 2.00MPa6. 5610查表确定铣加工的键槽

28、危险截面处疲劳强度的影响系数:K.：；=1.60=1.88 ；.：； = 0.75 % =0.73 1 =1.8则:420MPa- maxl6053.45MPa0.751.8-6. 63n.240MPa:-max1 882. 00MPa0.731.8二 83. 87n.c2 2八n.-6.613故满足强度条件。五、循环计算程序#in clude#in clude#defi ne pi 3.141592654#defi ne ip 0.017453292doubleL1=0.15,L2=0.53, L3=0.14,a=0.12,b=0.16,A0=45, n=400,P=5.5,i=0.65,R

29、=0.12,Fy=4000,F z=4000,Fby=0,Fbz=0,Fay,Faz,Fcy,Fdz,Fdy,Fcz,Fey,Fez,Me,Mby,Mbz,Mdy,Mdz,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;Me=9

30、549*P/n;Ft=Me/R;Fy=Ft*ta n(An *ip)/(cos(B n*ip);Fey=Ft*si n(A0*ip)-Fy*cos(A0*ip);Fez=Ft*cos(A0*ip)+Fy*si n(A0*ip);Md=-Me;Fdy=Fy;Fdz=-Fz;Mdy=Fdy*b;Mdz=Fdz*b;void waili() Fay=(Fdy*L3+Mdy-Fey*a)/(L1+L2);Fcy=(-Fdy*(L1+L2+L3)-Mdy-Fey*(L1+L2-a)/(L1+L2);Faz=(-Fdz*L3-Mdz-Fez*a)/(L1+L2);Fcz=(Fdz*(L1+L2+L3)+M

31、dz-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 qia ngdu()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/

32、2)*(Xmax/2)+Mmax*Mmax);Xmax 仁-Xmax/2+sqrt(Xmax/2)*(Xmax/2)+Mmax*Mmax); ddd=sqrt(Xmax1-Xmax2)*(Xmax1-Xmax2)+(Xmax2-Xmax3)*(Xmax2-Xmax3)+(Xmax3-Xmax1) *(Xmax3-Xmax1)/2)/xxx;ddd=pow(ddd,0.333333333333333333333);ddd=ddd/100;ddd=(float)ddd;Void 1n aodu()long double fffd1,fffdy,fffd z;fffdy=L1*Mby*L1*aMc/

33、(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+L2)/(0.294*eee);fffdy+=L3*aMc*(Mdy-2*Mcy)/(0.294*eee); fffdz=L1*Mbz*L1*aMc/(L1+L2)/0.147

34、/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*eee);fffdz+=a*(Mez*(2*aMc+(L1+L2-a)*aMc/(L1+L2)+aMc)+Mez*(2*aMc+(L1+L2-a)*aMc/(L 1+L2)/(0.294*eee);fffdz+=L3*aMc*(Mdz+2*Mcz)/(0.294*eee);fffd仁 pi*sqrt(fffdy*fffdy+fffdz*fffdz)/www/32/f

35、ffd;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;fff

36、ez+=(L2-a)*(Mbz*(2*aMe*L1/(L1+L2-a)+aMe)+Mez*(2*aMe*L1/(L1+L2-a)/(0.294*eee);fffez+=a*aMe*(Mcz-2*Mez)/(0.294*eee);fffe仁 pi*sqrt(fffey*fffey+fffez*fffez)/www/32/eee;ddd=pow(fffe1,0.25);ddd=ddd/1000;ddd=(float)ddd;voidzhuanjiaoc()long double aaac1,aaacy,aaacz;aaacy=L1*Mby*L1*aaac/(L1+L2)/0.147/eee;aaa

37、cy+=(L2-a)*(Mby*(2*aaaac*L1/(L1+L2)+aaaac*(L1+L2-a)/(L1+L2)+Mey*(2*aaaa c*(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*(

38、L1+L2-a)/(L1+L2)+Mez*(2*aaaa c*(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;voidfa nlib()long double ff

39、fb1,fffby,fffb z;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

40、*(Mcz+2*Mez)/(0.294*eee);fby0=fffby;fbz0=fffb z;fffb仁 pi*sqrt(fffby*fffby+fffbz*fffbz)*pi/www/big/big/big/big/640e11;fffbb=(float)fffb1;prin tf(B点的实际位移：S=%fmn,fffbb);/*void waili()Fby=fffb*SD*SD*SD*SD*(1-i*i*i*i)/FFF*10e12;Fbz=Fby*fbzO/fbyO;*/voidxun hua n()double AAAA, big1,big2,big3,big4;dowaili();qia ngdu();big1=ddd;big=big1;n aodue();big2=ddd;if(bigbig1)big=big2;n aodue();big3=ddd;if(bigbig3)big=big3;zhuanjiaoc();big4=(float)ddd;if(big=0.0001);prin tf(在超静定情况下的直径：D=%fmn,big);voidmai n()float big1=0,big2,big3,big4;printf(输入原始数据:n);sca nf(%f,%f,%f,%f,%f,%f,%f,%f