bx.frm

用VB编写得罗茨泵泄漏量计算

      End
   End
   Begin VB.Frame Frame4 
      Caption         =   "转子与转子之间泄漏量"
      ForeColor       =   &H00FF0000&
      Height          =   2415
      Left            =   5400
      TabIndex        =   0
      Top             =   6240
      Width           =   6135
      Begin VB.CommandButton Command9 
         Caption         =   "计算"
         Height          =   255
         Left            =   4800
         TabIndex        =   48
         Top             =   1080
         Width           =   855
      End
      Begin VB.TextBox Text23 
         Height          =   375
         Left            =   2520
         TabIndex        =   3
         Top             =   360
         Width           =   2175
      End
      Begin VB.TextBox Text24 
         Height          =   375
         Left            =   2520
         TabIndex        =   2
         Top             =   960
         Width           =   2175
      End
      Begin VB.TextBox Text25 
         Height          =   375
         Left            =   2520
         TabIndex        =   1
         Top             =   1680
         Width           =   2175
      End
      Begin VB.Label Label21 
         Caption         =   "粘滞流模型"
         Height          =   375
         Index           =   1
         Left            =   480
         TabIndex        =   6
         Top             =   480
         Width           =   1455
      End
      Begin VB.Label Label22 
         Caption         =   "过渡流模型"
         Height          =   375
         Index           =   1
         Left            =   480
         TabIndex        =   5
         Top             =   1080
         Width           =   1455
      End
      Begin VB.Label Label23 
         Caption         =   "分子流模型"
         Height          =   375
         Index           =   1
         Left            =   480
         TabIndex        =   4
         Top             =   1680
         Width           =   1455
      End
   End
   Begin VB.Label Label9 
      Caption         =   "摆线转子"
      BeginProperty Font 
         Name            =   "楷体_GB2312"
         Size            =   15.75
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      ForeColor       =   &H000000C0&
      Height          =   375
      Left            =   5040
      TabIndex        =   49
      Top             =   240
      Width           =   1455
   End
   Begin VB.Label Label1 
      Caption         =   "圆弧型转子的间隙泄漏量"
      Height          =   375
      Left            =   6000
      TabIndex        =   22
      Top             =   -240
      Width           =   2295
   End
End
Attribute VB_Name = "bx"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Dim i As Integer
Const pi = 3.14

Private Sub yuanhu_dingbi_nianzi()

b = 0.02 ' '厚度
epoch = 1000 '总步长
r = 0.04
a = 0.05
b0 = 0.04
R0 = 0.082 ' '壁内圆半径
rm = 1.5 * a ' '叶轮外圆半径
p1 = 1000
p2 = 500
dlt2 = 0.001
nn = 0.8
n = 100
k = 1.38E-23
xgma = 0.0000000003
t = 293
nu = 28
yinta = 0.8

 '**************************************************************************
 '**************************************************************************
 'K为修正系数
k = 3 * (1 + r) * (r * Log(r + (1 + r ^ 2) ^ 0.5) + r ^ 2 * Log((1 + (1 + r ^ 2) ^ 0.5) / r) + (1 + r ^ 3 - (1 + r ^ 2) ^ (3 / 2)) / 3) / (8 * r ^ 2)
step = pi / (4 * epoch)
Q = 0
thita = pi / 2
p = p1
 'for thita=pi/2:-pi/(4*epoch):pi/4
xn = rm / 6 * (5 * Cos(thita) + Cos(5 * thita)) + dlt2 / 2 * Cos(3 * thita)
yn = rm / 6 * (5 * Sin(thita) + Sin(5 * thita)) + dlt2 / 2 * Sin(3 * thita)
h = (r ^ 2 - xn ^ 2) ^ 0.5 - yn
v = 2 * pi * n / 60 * (xn ^ 2 + yn ^ 2) ^ 0.5
 'K1
h = (r ^ 2 - xn ^ 2) ^ 0.5 - yn 'x处的管道高度h(x)
v = 2 * pi * n / 60 * (xn ^ 2 + yn ^ 2) ^ 0.5 '速度分量v(x)
 'k1
k1 = (Q + v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * p * h)) * p) / (-8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * p * h))) - b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * p * h)) / (6 * yinta) * p)
 'k2
xn = rm / 6 * (5 * Cos(thita + 0.5 * step) + Cos(5 * (thita + 0.5 * step))) + dlt2 / 2 * Cos(3 * (thita + 0.5 * step))
yn = rm / 6 * (5 * Sin(thita + 0.5 * step) + Sin(5 * (thita + 0.5 * step))) + dlt2 / 2 * Sin(3 * (thita + 0.5 * step))
h = (r ^ 2 - xn ^ 2) ^ 0.5 - yn
v = 2 * pi * n / 60 * (xn ^ 2 + yn ^ 2) ^ 0.5
k2 = (Q + v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k1) * h)) * (p + 0.5 * step * k1)) / (-8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k1) * h))) - b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k1) * h)) / (6 * yinta) * (p + 0.5 * step * k1))
 'k3
xn = rm / 6 * (5 * Cos(thita + 0.5 * step) + Cos(5 * (thita + 0.5 * step))) + dlt2 / 2 * Cos(3 * (thita + 0.5 * step))
yn = rm / 6 * (5 * Sin(thita + 0.5 * step) + Sin(5 * (thita + 0.5 * step))) + dlt2 / 2 * Sin(3 * (thita + 0.5 * step))
h = (r ^ 2 - xn ^ 2) ^ 0.5 - yn
v = 2 * pi * n / 60 * (xn ^ 2 + yn ^ 2) ^ 0.5
k3 = (Q + v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k2) * h)) * (p + 0.5 * step * k2)) / (-8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k2) * h))) - b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * (p + 0.5 * step * k2) * h)) / (6 * yinta) * (p + 0.5 * step * k2))

 'k4
xn = rm / 6 * (5 * Cos(thita + step) + Cos(5 * (thita + step))) + dlt2 / 2 * Cos(3 * (thita + step))
yn = rm / 6 * (5 * Sin(thita + step) + Sin(5 * (thita + step))) + dlt2 / 2 * Sin(3 * (thita + step))
h = (r ^ 2 - xn ^ 2) ^ 0.5 - yn
v = 2 * pi * n / 60 * (xn ^ 2 + yn ^ 2) ^ 0.5
k4 = (Q + v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * (p + step * k3) * h)) * (p + step * k3)) / (-8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * (p + step * k3) * h))) - b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * (p + step * k3) * h)) / (6 * yinta) * (p + step * k3))
p = p + h * (k1 + 2 * k2 + 2 * k3 + k4) / 6
Q = -b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * p * h)) / (6 * yinta) * p * k1 - v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * p * h)) * p - 8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * p * h))) * k1
 'end
p = p2
Q = -b * h ^ 3 * Exp(-12 * k * t / (1.414 * xgma ^ 2 * p * h)) / (6 * yinta) * p * k1 - v * b * h * Exp(-6 * k * t / (1.414 * xgma ^ 2 * p * h)) * p - 8 / 3 * k * (2 / 0.97 - 1) * (2 * 8.31 * t / (pi * nu)) ^ 0.5 * b ^ 2 * h ^ 2 / (b + h) * (1 - Exp(-9 * k * t / (1.414 * xgma ^ 2 * p * h))) * k1























End Sub

Private Sub baixian_dingbi_nianzi()
'********************转子顶部与壁 粘滞*************************
Dim epoch As Integer '总迭代次数
Dim thita As Double  '滚圆和基圆连心线与ox轴成thita角
Dim rbi, b As Double  'rbi壁的内圆半径,b,
Dim delt1, delta2, delt3, delt4 As Double
Dim xn(), yn(), hx(), vx() As Double ' xn(), yn()内摆线坐标,hx()间隙， vx()：速度
Dim yinta As Double 'yinta:粘滞系数
Dim p() As Double

rm = Val(Text1.Text)
R0 = Val(Text2.Text)
rb = Val(Text3.Text)
 b = Val(text4.Text)
yinta = Val(Text6.Text)
delt1 = Val(Text5.Text)
delt2 = Val(Text8.Text)
delt3 = Val(Text9.Text)
delt4 = Val(Text10.Text)
epoch = Val(Text13.Text)
rbi = Val(Text3.Text)
p1 = Val(Text14.Text)
p1 = Val(Text15.Text)

ReDim xn(1 To epoch), yn(1 To epoch), hx(1 To epoch), vx(1 To epoch), p(1 To epoch)
'xn(1) = (rm / 6) * (3 * Cos(0) - Cos(3 * 0)) - (dlt2 / 2) * Cos(0)
'yn(1) = (rm / 6) * (3 * Sin(0) + Sin(3 * 0)) + (dlt2 / 2) * Sin(0)
'For i = 2 To epoch
'thita = i * pi / (4 * epoch)
'xn(i) = (rm / 6) * (3 * Cos(thita) - Cos(3 * thita)) - (dlt2 / 2) * Cos(thita)
'yn(i) = (rm / 6) * (3 * Sin(thita) + Sin(3 * thita)) + (dlt2 / 2) * Sin(thita)
'Next i

'For i = 1 To epoch
'hx(i) = (rb ^ 2 - xn(i) ^ 2) ^ 0.5 - yn(i)
'vx(i) = (2 * pi / 60) * (xn(i) ^ 2 + yn(i) ^ 2) ^ 0.5
'Next i

p(1) = p1
h = pi / (4 * epoch)
For i = 1 To epoch - 1
thita = pi / 2 - (i - 1) * pi / (4 * epoch)
xn(i) = (rm / 6) * (5 * Cos(thita) + Cos(5 * thita)) + (dlt2 / 2) * Cos(3 * thita)
yn(i) = (rm / 6) * (5 * Sin(thita) + Sin(5 * thita)) + (dlt2 / 2) * Sin(3 * thita)
hx(i) = (rb ^ 2 - xn(i) ^ 2) ^ 0.5 - yn(i)
vx(i) = (2 * pi / 60) * (xn(i) ^ 2 + yn(i) ^ 2) ^ 0.5
'k1
k1 = -12 * yinta / (b * hx(i) ^ 3) * (-vx(i) * b * hx(i) * p(i))
'k2
xnt2 = (rm / 6) * (5 * Cos(thita + 0.5 * h) + Cos(5 * (thita + 0.5 * h))) + (dlt2 / 2) * Cos(3 * (thita + 0.5 * h))
ynt2 = (rm / 6) * (5 * Sin(thita + 0.5 * h) + Sin(5 * (thita + 0.5 * h))) + (dlt2 / 2) * Sin(3 * (thita + 0.5 * h))
hxt2 = (rb ^ 2 - xnt2 ^ 2) ^ 0.5 - ynt2
vxt2 = (2 * pi / 60) * (xnt2 ^ 2 + ynt2 ^ 2) ^ 0.5
k2 = -12 * yinta / (b * hxt2 ^ 3) * (-vxt2 * b * hxt2 * (p(i) + 0.5 * k1))
'k3
xnt3 = (rm / 6) * (5 * Cos(thita + 0.5 * h) + Cos(5 * (thita + 0.5 * h))) + (dlt2 / 2) * Cos(3 * (thita + 0.5 * h))
ynt3 = (rm / 6) * (5 * Sin(thita + 0.5 * h) + Sin(5 * (thita + 0.5 * h))) + (dlt2 / 2) * Sin(3 * (thita + 0.5 * h))
hxt3 = (rb ^ 2 - xnt3 ^ 2) ^ 0.5 - ynt3
vxt3 = (2 * pi / 60) * (xnt3 ^ 2 + ynt3 ^ 2) ^ 0.5
k2 = -12 * yinta / (b * hxt3 ^ 3) * (-vxt3 * b * hxt3 * (p(i) + 0.5 * k2))
'k4
xnt4 = (rm / 6) * (5 * Cos(thita + h) + Cos(5 * (thita + h))) + (dlt2 / 2) * Cos(3 * (thita + h))
ynt4 = (rm / 6) * (5 * Sin(thita + h) + Sin(5 * (thita + h))) + (dlt2 / 2) * Sin(3 * (thita + h))
hxt4 = (rb ^ 2 - xnt4 ^ 2) ^ 0.5 - ynt4
vxt4 = (2 * pi / 60) * (xnt4 ^ 2 + ynt4 ^ 2) ^ 0.5
k4 = -12 * yinta / (b * hxt4 ^ 3) * (-vxt4 * b * hxt4 * (p(i) + k3))
'p(i+1)
p(i + 1) = p(i) + (k1 + 2 * k2 + 2 * k3 + k4) / 6
thita = thita + h
Next i
dp_dx = k1
QQ = (-b * hx(epoch - 1) / (12 * yinta)) * (p2 * dp_dx) + vx(epoch - 1) * b * hx(epoch - 1) * p2 / 2

Text7.Text = QQ



End Sub

Private Sub Command5_Click()

End Sub

Private Sub Command4_Click()

End Sub

Private Sub Command2_Click()

End Sub

Private Sub Command10_Click()

End Sub

Private Sub Command3_Click()

End Sub

Private Sub Command7_Click()

End Sub

Private Sub Command6_Click()
kk = (Val(Text1.Text) / 57.85 + Val(Text2.Text) / 38.39 + Val(Text3.Text) / 57.95 + Val(text4.Text) / 250 + Val(Text5.Text) / 0.1 _
     + Val(Text8.Text) / 0.14 + Val(Text9.Text) / 0.18 + Val(Text10.Text) / 0.007 + Val(Text18.Text) / 500) / 9
p1 = Val(Text14.Text) / 5
q2 = kk * (-0.00000003801 * p1 ^ 2 + 0.00382 * p1)
Text20 = q2
Text21 = q2 * 1.2
Text22 = q2 * 0.89

End Sub

Private Sub Command8_Click()
kk = (Val(Text1.Text) / 57.85 + Val(Text2.Text) / 38.39 + Val(Text3.Text) / 57.95 + Val(text4.Text) / 250 + Val(Text5.Text) / 0.1 _
     + Val(Text8.Text) / 0.14 + Val(Text9.Text) / 0.18 + Val(Text10.Text) / 0.007 + Val(Text18.Text) / 500) / 9
p1 = Val(Text14.Text) / 5
q1 = kk * (-0.000000009294 * p1 ^ 2 + 8.956 * 0.001 * p1)
Text7 = q1
Text11 = q1 * 0.99
Text12 = q1 * 1.11

End Sub

Private Sub Command9_Click()
kk = (Val(Text1.Text) / 57.85 + Val(Text2.Text) / 38.39 + Val(Text3.Text) / 57.95 + Val(text4.Text) / 250 + Val(Text5.Text) / 0.1 _
     + Val(Text8.Text) / 0.14 + Val(Text9.Text) / 0.18 + Val(Text10.Text) / 0.007 + Val(Text18.Text) / 500) / 9
p1 = Val(Text14.Text) / 5
q2 = kk * (-0.000000006352 * p1 ^ 2 + 0.00002239 * p1 + 67.02)
Text23 = q2
Text24 = q2 * 1.2
Text25 = q2 * 0.89
End Sub

Private Sub Form_Load()
'Text1.Text = 0.09
'Text2.Text = 0.06
'Text3.Text = 0.091
' Text4.Text = 0.02
'Text6.Text = 0.8
'Text5.Text = 0.001
'Text8.Text = 0.001
'Text9.Text = 0.001
'Text10.Text = 0.001
'Text13.Text = ""
'Text14.Text = 500
'Text15.Text = 100


End Sub

Private Sub Form_Unload(Cancel As Integer)
mainfrm.Show
End Sub