
# -*- coding: UTF-8 -*-
import math
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from pylab import mpl
mpl.rcParams['font.sans-serif'] = ['SimHei']
mpl.rcParams['axes.unicode_minus'] = False
# ===================================================================
# 算法——井身结构设计方法
# ===================================================================
# def bonds():
# rho_d>=rho_p_max+Delta_rho
#
# rho_d_max>=max(rho_p_max+Delta_rho,rho_c_max)
#
# rho_bn_max = rho_d_max + S_g
# rho_ff = rho_f - S_f
# rho_bn_max<=rho_ff_min
#
# Delta_p=0.00981*(rho_d_max-rho_p_min)*D_n<=Delta_p_N(Delta_p_A)
#
# rho_ba_max = rho_d_max + D_m * S_k / D_x
# rho_ba_max<=rho_ffDmin
print(np.degrees(np.arctan(0.13)))
df0=pd.read_excel(r'C:\Users\Administrator\Desktop\upupup\2-petrol_des