2021-07-28

油管问题

import numpy as np
import pandas as pd

#初始值
T = 0.288   #单向阀开启时间
P[1] = 100  #初始压强
rou = 0.85  #初始密度
dt = 0.001  #时间步长
Tall = 30000  #总模拟时间
t0 = 0      #喷油开始时刻
V = 12500 * pi  #油管体积
m = rou * V  #初始质量
P_I = 160   #入口处的压强
rou_I = 0.87143  #160MPa时的密度
Pe = 100    #期望压强
C = 0.85    #流量系数
A = pi * 0.7 * 0.7  #小孔面积
Itimes = 1  #进油次数
Etimes = 1  #出油次数
I_prev = 0  #前一时刻的进油量
E_prev = 0  #前一时刻的出油量

Tin_start = 0  #进油开始时刻
Tin_end = T  #进油结束时刻

Tout_start = t0  #第一次喷油开始时刻
Tout_1 = tout_start + 0.2  #出油速率第一段
Tout_2 = Tout_1+2  #出油速率第二段
Tout_end = Tout_2 + 0.2  #出油速率第三段，完成一个周期的出油

t0 = np.arange(0,97.6,0.001)
for i in range(1,n1+1,1):
    ou1[i] = [0 + (i-1) * (T+10) , T + (i-1) * (T+10)]
for j in range(1,n2+1,1):
    ou2[j] = [t0 + 100 * (j-1),2.4 + t0 + 100 * (j-1)]

for i in range(1,n1+1,1):
    for t in range(ou1[i]):
        Q1[t] = C * A * sqrt(2 * dP / rou_I)
        I[t] = Q1[t]
for j in range(1,n2+1,1):
    for t in range(ou2[j]):
        if t1 in range(0,0.2,0.001):
            Q2[t] = 100 * t1
        if t1 in range(0.2,2.2,0.001):
            Q2[t] = 20
        if t1 in range(2.2,2.4,0.001):
            Q2[t] =240 - 100 * t1
        E[t] = Q2[t]
    t1 = t - to + 100 * (j-1)

m[0] = rou[0] * V
m[t+0.001]-m[t] = (rou_I * I[t] - rou[t] * E[t])
dm / dt = rou_I * I - rou[t] * V