' 程序名称：圆柱直齿轮设计
' Author:ncf
' Date:2016-10-29
' Email:1476879092@qq.com
' log:
' 2016-10-29:Complete
Public Check3_flag As BooleanPrivate Sub Check1_Click()If Check2.Value = 1 Then Check2.Value = 0
End SubPrivate Sub Check2_Click()If Check1.Value = 1 Then Check1.Value = 0
End SubPrivate Sub Check3_Click()If Check3.Value = 1 ThenCheck3_flag = True    '开式为TrueCheck3.Caption = "开式传动"ElseCheck3_flag = False    '开式为TrueCheck3.Caption = "闭式传动"End If
End SubPrivate Sub Command1_Click()
'数据的输入p_0 = Val(Text1.Text)    '输入的功率N_1 = Val(Text5.Text)    '小齿轮转速S_a = Val(Text21.Text)    '工作机载荷性质代码P_w = Val(Text9.Text)    '原动机代码W_N = Val(Text10.Text)    '工作两班H_0 = Val(Text22.Text)    '一班的时间D_A = Val(Text11.Text)    '工作时间j = Val(Text23.Text)    'j 为每转一圈齿轮啮合次数MA1 = Val(Text6.Text)    '小齿轮材料热处理代码MA2 = Val(Text18.Text)    '大齿轮材料热处理代码HB_x = Val(Text7.Text)    '小齿轮硬度HB_d = Val(Text19.Text)    '大齿轮硬度K_ht = Val(Text8.Text)     '接触试选载荷系数K_ft = Val(Text20.Text)     '弯曲试选载荷系数FD = Val(Text15.Text)     '齿宽系数RF = 1    '设计准则的选用标志,注：这里我们直接同时采用两个标准设计，不采用一个设计，一个校核pai = 3.1415926536    '圆周率piha = Val(Text2.Text)     '齿顶高系数alpha = Val(Text14.Text) / 180 * pai   '分度圆压力角JZ = Val(Text3.Text)     '精度等级CJ = Val(Text13.Text)     '支撑结构选择'确定计算方式flag_js = 0If Check1.Value = 1 Then flag_js = 1If Check2.Value = 1 Then flag_js = 2If flag_js <> 0 Then'检测空缺数据flag_space = 1If flag_js = 1 Then    '齿数已知ZV1 = Val(Text4.Text)    '小齿轮齿数ZV2 = Val(Text16.Text)    '大齿轮齿数Text12.Text = 0If ZV1 = 0 Or ZV2 = 0 Or p_0 = 0 Or N_1 = 0 Or S_a = 0 Or P_w = 0 Or W_N = 0 Or H_0 = 0 Or D_A = 0 Or j = 0 Or MA1 = 0 Or MA2 = 0 Or HB_x = 0 Or HB_d = 0 Or K_ht = 0 Or K_ft = 0 Or FD = 0 Or ha = 0 Or alpha = 0 Or JZ = 0 Or CJ = 0 ThenMsgBox ("请确保设计参数填写正确")Elseu = ZV2 / ZV1Text12.Text = uText17.Text = Val(Text5.Text) / u    '大齿轮转速flag_space = 0    '数据无空缺End IfElse    '传动比已知If Check3_flag = False ThenIf HB_x > 350 Then ZV1 = 18 Else ZV1 = 24    '闭式硬齿面默认18齿，闭式软齿24齿ElseZV1 = 18    '开式18齿End Ifu = Val(Text12.Text)    '齿数比If u = 0 Or p_0 = 0 Or N_1 = 0 Or S_a = 0 Or P_w = 0 Or W_N = 0 Or H_0 = 0 Or D_A = 0 Or j = 0 Or MA1 = 0 Or MA2 = 0 Or HB_x = 0 Or HB_d = 0 Or K_ht = 0 Or K_ft = 0 Or FD = 0 Or ha = 0 Or alpha = 0 Or JZ = 0 Or CJ = 0 ThenMsgBox ("请确保设计参数填写正确")ElseZV2 = Round(ZV1 * u)    '大齿轮齿数Text4.Text = ZV1  '小齿轮齿数Text16.Text = ZV2  '小齿轮齿数Text17.Text = Val(Text5.Text) / u    '大齿轮转速flag_space = 0    '数据无空缺End IfEnd If'计算If flag_space = 0 Then    '检测空缺数据'通用数据的计算'转矩的计算T_1 = 9.55 * (10 ^ 6) * p_0 / N_1'许用应力的计算'应力循环次数的计算N1 = 60 * j * N_1 * (W_N * H_0 * D_A)N2 = N1 / u'寿命系数Call Zichengxu1(MA1, N1, KB1)Call Zichengxu1(MA2, N2, KB2)Call Zichengxu3(MA1, N1, KH1)Call Zichengxu3(MA2, N2, KH2)'许用应力Call Zichengxu2(MA1, HB_x, KB1, CB1)Call Zichengxu2(MA2, HB_d, KB2, CB2)Call Zichengxu4(MA1, HB_x, KH1, CH1)Call Zichengxu4(MA2, HB_d, KH2, CH2)'计算重合度Call arccoa(ZV1 * Cos(alpha) / (ZV1 + 2 * ha), alpha1)Call arccoa(ZV2 * Cos(alpha) / (ZV2 + 2 * ha), alpha2)e_alpha = (ZV1 * (Tan(alpha1) - Tan(alpha)) + ZV2 * (Tan(alpha2) - Tan(alpha))) / 2 / pai'使用系数KA的查询Call KA_find(S_a, P_w, KA)'按接触强度设计'计算重合度系数Z_e = Sqr((4 - e_alpha) / 3)'计算ze zhU1 = 0.3    '泊松比U2 = 0.3E1 = 20.6 * 10000    '锻钢的弹性模量E2 = 20.6 * 10000Call Zichengxu7(alpha, U1, U2, E1, E2, ZH, ZE)'带入计算时数据的选用If CH1 > CH2 Then CH = CH2 Else CH = CH1    '选接触疲劳许用应力小的'试算分度圆直径d_1t = ((2 * K_ht * T_1 / FD) * ((u + 1) / u) * ((ZH * ZE * Z_e) / CH) ^ 2) ^ (1 / 3)'计算系数KVCall KV_js(ZV1, d_1t, N_1, JZ, KV)'调整小齿轮分度圆直径Call Zichengxu8(HB_x, HB_d, CJ, FD, KHB, KFB)KH = KV * KA * KHB    '忽略齿间载荷分配系数K_alphaD1 = d_1t * ((KH / K_ht) ^ (1 / 3))m_H = D1 / ZV1'按齿根弯曲强度设计'计算重合度系数Y_e = 0.25 + 0.75 / e_alpha'修正系数的计算Call Zichengxu5(ZV1, YFA1)Call Zichengxu5(ZV2, YFA2)Call Zichengxu6(ZV1, YSA1)Call Zichengxu6(ZV2, YSA2)YK1 = YFA1 * YSA1 / CB1YK2 = YFA2 * YSA2 / CB2'带入计算时数据的选用If YK2 > YK1 Then YK = YK2 Else YK = YK1    '在弯曲疲劳强度的设计中也是选小的'试算模数m_1t = (2 * K_ft * T_1 * Y_e / (FD * ZV1 ^ 2) * YK) ^ (1 / 3)d_1t = m_1t * ZV1'计算系数KVCall KV_js(ZV1, d_1t_F, N_1, JZ, KV)'调整小齿轮模数Call Zichengxu8(HB_x, HB_d, CJ, FD, KHB, KFB)KF = KV * KA * KFBm_F = m_1t * ((KF / K_ft) ^ (1 / 3))'两种设计结果的处理'已知传动比If flag_js = 2 Then    '传动比已知If m_H > m_F ThenCall m_find(m_F, m_b)    '标准化模数Z_1 = Round(D1 / m_b)    '修改小齿轮齿数Z_2 = Round(Z_1 * u)    '大齿轮齿数d_1 = Z_1 * m_b    '小齿轮直径d_2 = Z_2 * m_b    '大齿轮直径A = (d_1 + d_2) / 2    '中心距B_d = FD * d_1    '大齿轮齿宽B_x = B_d + 7.5    '小齿轮齿宽ElseZ_1 = ZV1 '不修改齿数Z_2 = ZV2Call m_find(m_F, m_b)    '标准化模数d_1 = Z_1 * m_b    '小齿轮直径d_2 = Z_2 * m_b    '大齿轮直径A = (d_1 + d_2) / 2    '中心距B_d = FD * d_1    '大齿轮齿宽B_x = B_d + 7.5    '小齿轮齿宽End IfElseIf flag_js = 1 Then    '齿数已知If m_H > m_F Thenm_ = m_HElsem_ = m_FEnd IfZ_1 = ZV1 '不修改齿数Z_2 = ZV2Call m_find(m_, m_b)    '标准化模数d_1 = Z_1 * m_b    '小齿轮直径d_2 = Z_2 * m_b    '大齿轮直径A = (d_1 + d_2) / 2    '中心距B_d = FD * d_1    '大齿轮齿宽B_x = B_d + 7.5    '小齿轮齿宽End If'显示输出Text24.Text = m_bText25.Text = AText26.Text = Z_1Text27.Text = Z_2Text28.Text = d_1Text29.Text = d_2Text30.Text = B_xText31.Text = B_dEnd IfElseMsgBox ("请选择：计算方式")End If
End Sub'自定义子程序
'KA的获取
Sub KA_find(S_a, P_w, KA)S_a = S_a - 1P_w = P_w - 1Dim KA_array1()KA_array1 = Array(1#, 1.1, 1.25, 1.5)Dim KA_array2KA_array2 = Array(1.25, 1.35, 1.5, 1.75)Dim KA_array3KA_array3 = Array(1.5, 1.6, 1.75, 2)Dim KA_array4KA_array4 = Array(1.75, 1.85, 2, 2.25)If S_a = 0 ThenKA = KA_array1(P_w)ElseIf S_a = 1 ThenKA = KA_array2(P_w)ElseIf S_a = 2 ThenKA = KA_array3(P_w)ElseIf S_a = 3 ThenKA = KA_array4(P_w)ElseEnd If
End Sub
'm的查询
Sub m_find(m, m_b)Dim m_bzm_bz = Array(1, 1.25, 1.5, 2, 2.5, 3, 4, 5, 6, 8, 10)num_m = UBound(m_bz) - LBound(m_bz)For i = 0 To num_mIf m < m_bz(i) Thenm_b = m_bz(i)Exit ForEnd IfNext i
End Sub
'反余弦函数
Sub arccoa(X1, alpha1)alpha1 = Atn(-X1 / Sqr(-X1 * X1 + 1)) + 2 * Atn(1)
End Sub
'计算KV
Sub KV_js(Z1, D1, N0, JZ, KV)V = Z1 * 3.14159 * D1 * NO / 6000000!KV = (0.0104 + 0.003 * (JZ - 6)) * JZ * V + 1
End Sub'老师提供的子程序
Sub Zichengxu1(MA, N, KB)If MA > 27 Then GoTo 1035If N <= 10000! Then GoTo 1025If N > 10000! And N < 3000000! Then GoTo 1030If N > 3000000! Then KB = 1!: GoTo 1125
1025:KB = 2.5: GoTo 1125
1030:KB = (3! * 10 ^ 6 / N) ^ (0.115): GoTo 1125
1035:If MA > 36 Or MA < 27 Then GoTo 1065If N <= 1000! Then GoTo 1055If N > 1000! And N < 3000000! Then GoTo 1060If N > 3000000! Then KB = 1!: GoTo 1125
1055:KB = 2.5: GoTo 1125
1060:KB = (3! * 10 ^ 6 / N) ^ (0.115): GoTo 1125
1065:If MA > 46 Or MA < 36 Then GoTo 1095If N <= 1000! Then GoTo 1085If N > 1000! And N < 3000000! Then GoTo 1090If N > 3000000! Then KB = 1!: GoTo 1125
1085:KB = 1.6: GoTo 1125
1090:KB = (3! * 10 ^ 6 / N) ^ (0.05900001): GoTo 1125
1095:If MA < 46 Then GoTo 1125If N < 1000! Then GoTo 1115If N > 1000! And N < 3000000! Then GoTo 1120If N > 3000000! Then KB = 1!: GoTo 1125
1115:KB = 1.1: GoTo 1125
1120:KB = (3! * 10 ^ 6 / N) ^ (0.012)
1125:
End Sub
Sub Zichengxu2(MA, HB, KB, CB)If MA > 15 Then GoTo 1160On MA - 10 GoTo 1140, 1145, 1150, 1155
1140:CB = 0.66 * (HB - 140) + 120: GoTo 1230
1145:CB = 0.6850001 * (HB - 140) + 350: GoTo 1230
1150:CB = 0.2 * (HB - 110) + 300: GoTo 1230
1155:CB = 0.6 * (HB - 110) + 340: GoTo 1230
1160:If MA > 27 Then GoTo 1190On MA - 22 GoTo 1170, 1175, 1180, 1185
1170:CB = 0.37 * (HB - 115) + 340: GoTo 1230
1175:CB = 0.75 * (HB - 115) + 380: GoTo 1230
1180:CB = 0.8 * (HB - 190) + 520: GoTo 1230
1185:CB = 0.47 * (HB - 190) + 460: GoTo 1230
1190:If MA > 36 Then GoTo 1210On MA - 33 GoTo 1200, 1205
1200:CB = 10 * (HB - 48) + 640: GoTo 1230
1205:CB = 3 * (HB - 56) + 840: GoTo 1230
1210:If MA = 44 Or MA = 45 Then GoTo 1220If MA = 54 Or MA = 55 Then CB = 2! * (HB - 54) + 640: GoTo 1230
1220:CB = 2! * (HB - 54) + 740
1230:CB = KB * CB / 1.375
End Sub
Sub Zichengxu3(MA, N, KH)If MA > 37 And FI = 1 Then GoTo 1545If N <= 600000! Then GoTo 1530If N > 600000! And N < 10000000# Then GoTo 1535If N > 10000000# And N < 1000000000# Then GoTo 1540If N > 1000000000# Then KH = 1!: GoTo 1640
1530:KH = 1.6: GoTo 1640
1535:KH = (3! * 10 ^ 8 / N) ^ (0.0756): GoTo 1640
1540:KH = (10 ^ 9 / N) ^ (0.05700001): GoTo 1640
1545:If MA > 37 Or MA < 22 Then GoTo 1575If N < 100000! Then GoTo 1565If N > 100000! And N < 50000000# Then GoTo 1570If N > 50000000# Then KH = 1!: GoTo 1640
1565:KH = 1.6: GoTo 1640
1570:KH = (5 * 10 ^ 7 / N) ^ (0.0756): GoTo 1640
1575:If MA > 47 Or MA < 42 Then GoTo 1605If N <= 100000! Then GoTo 1595If N > 100000! And N < 2000000! Then GoTo 1600If N > 2000000! Then KH = 1!: GoTo 1640
1595:KH = 1.3: GoTo 1640
1600:KH = (2 * 10 ^ 6 / N) ^ (0.0875): GoTo 1640
1605:If MA < 52 Then GoTo 1640If N <= 100000! Then GoTo 1625If N > 100000! And N < 2000000! Then GoTo 1630If N >= 2000000! Then KH = 1!: GoTo 1640
1625:KH = 1.1: GoTo 1640
1630:KH = (2 * 10 ^ 6 / N) ^ (0.0318)
1640:
End Sub
Sub Zichengxu4(MA, HB, KH, CH)If MA > 15 Then GoTo 1680On MA - 10 GoTo 1663, 1665, 1670, 1675
1663:CH = 1.154 * (HB - 140) + 300: GoTo 1800
1665:CH = 1.406 * (HB - 140) + 400: GoTo 1800
1670:CH = 190 + HB: GoTo 1800
1675:CH = 260 + HB: GoTo 1800
1680:If MA > 27 Then GoTo 1710
1685 On MA - 22 GoTo 1690, 1695, 1700, 1705
1690:CH = 290 + HB: GoTo 1800
1695:CH = 350 + HB: GoTo 1800
1700:CH = 1.465 * (HB - 190) + 620: GoTo 1800
1705:CH = 1.4 * (HB - 190) + 550: GoTo 1800
1710:If MA > 36 Then GoTo 1745If MA = 34 Then GoTo 1730If MA = 35 Then GoTo 1735If MA = 35 And HB > 58 Then GoTo 1740
1730:CH = 25 * (HB - 48) + 1120: GoTo 1800
1735:CH = 60 * (HB - 56) + 1350: GoTo 1800
1740:CH = 1470: GoTo 1800
1745:If MA > 50 Then GoTo 1775If MA = 44 Or MA = 45 Then CH = 11.1 * (HB - 36) + 900If MA = 44 And HB > 45 Then GoTo 1765If MA = 45 And HB > 45 Then GoTo 1770GoTo 1800
1765:CH = 1000: GoTo 1800
1770:CH = 1000: GoTo 1800
1775:If MA = 54 Or MA = 55 Then CH = 18 * (HB - 54) + 1225
1780 If MA = 54 And HB > 59 Then GoTo 1790
1786 If MA = 55 And HB > 59 Then GoTo 1795
1787 GoTo 1800
1790:CH = 1315: GoTo 1800
1795:CH = 1315
1800:CH = KH * CH
End Sub
Sub Zichengxu5(ZV, YFA)YFA = 2.984508 - 0.014134 * ZV + 0.000051 * ZV ^ 2
End Sub
Sub Zichengxu6(ZV, YSA)YSA = 1.472047 + 0.00497 * ZV - 0.000016 * ZV ^ 2
End Sub
Sub Zichengxu7(AL, U1, U2, E1, E2, ZH, ZE)ZH = Sqr(2 / (Sin(AL) * Cos(AL)))ZE = 0.564 * Sqr(1 / ((1 - U1 ^ 2) / E1 + (1 - U2 ^ 2) / E2))
End Sub
'修改程序8
Sub Zichengxu8(H1, H2, CJ, FD, KHB, KFB)If H1 > 350 And H2 > 350 Then GoTo 2450On CJ GoTo 2415, 2420, 2425, 2430, 2435, 2440, 2445
2415:KHB = 0.24 * FD + 0.4456 * FD ^ 2 + 1.0054: GoTo 2490
2420:KHB = 0.29184 * FD + 0.06218 * FD ^ 2 + 0.99569: GoTo 2490
2425:KHB = 0.10482 * FD + 0.04448 * FD ^ 2 + 1.0036: GoTo 2490
2430:KHB = 0.0395 * FD + 0.04968 * FD ^ 2 + 1.013: GoTo 2490
2435:KHB = 0.02576 * FD + 0.04125 * FD ^ 2 + 1.011: GoTo 2490
2440:KHB = -0.007 * FD + 0.03977 * FD ^ 2 + 1.157: GoTo 2490
2445:KHB = -0.0177 * FD + 0.0337 * FD ^ 3 + 1.00977: GoTo 2490
2450 On CJ GoTo 2455, 2460, 2465, 2470, 2475, 2480, 2485
2455:KHB = 1.05 * FD + 1: GoTo 2490
2460:KHB = 0.41026 * FD + 0.51246 * FD ^ 2 + 1.0033: GoTo 2490
2465:KHB = 0.29503 * FD + 0.07772 * FD ^ 2 + 1.00587: GoTo 2490
2470:KHB = 0.18042 * FD + 0.08549999 * FD ^ 2 + 1.00498: GoTo 2490
2475:KHB = 0.0907 * FD + 0.09494 * FD ^ 2 + 1.00352: GoTo 2490
2480:KHB = 0.01265 * FD + 0.08847999 * FD ^ 2 + 1.0095: GoTo 2490
2485:KHB = -0.009289999 * FD + 0.06270001 * FD ^ 2 + 1.00862
2490:KFB = 1.5 * KHB - 0.5
End Sub

百度云链接：