K = 1.5 H = 7.5 S = 10 V1 = 300 T2 = 5 V2 = 10 theta = np.radians(3)重新回答

好的，根据您提供的数值，我们可以重新编写代码：

import matplotlib.pyplot as plt
import numpy as np
import math

def calculate_pressure(K, H, S, V1, T2, V2, theta, d):
    if H >= S:
        P = (K * H) / (S * (d/V1 + T2) * V2 * math.tan(theta))
    else:
        P = (K * (S - (S - H) / 2)) / (S * (d/V1 + T2) * V2 * math.tan(theta))
    return P

K = 1.5
H = 7.5
S = 10
V1 = 300
T2 = 5
V2 = 10
theta = np.radians(3)

d_array = np.linspace(0, 1000, 10000)
P_array = np.array([calculate_pressure(K, H, S, V1, T2, V2, theta, d) for d in d_array])
Z_array = (d_array / V1 + T2) * V2 + d_array

mask = P_array >= 0.8

plt.plot(d_array[mask], Z_array[mask])

plt.xlabel('d')
plt.ylabel('Z')
plt.show()

这段代码会生成一个图形窗口，显示出 P 大于等于 0.8 时，Z 和 d 的关系图。您可以根据需要调整常量和数组的值，来得到不同的图形结果。