pythonsymy：简单的操作即使在一天之后也会忙着运行客观的_Python_Sympy_Integral

pythonsymy：简单的操作即使在一天之后也会忙着运行客观的

python

pythonsymy：简单的操作即使在一天之后也会忙着运行客观的,python,sympy,integral,Python,Sympy,Integral,我试图用符号来解一个常数系数的积分（a_w，b_IN_w，c_IN_w），它们由简单的代数表达式组成我试过的我已经试着运行下面给出的代码一整天了，但是当我一天后检查它时，它仍然在运行。我以前使用过Symphy，我知道它可能无法解决一些复杂的操作，在这些操作中，它会抛出某种指示问题的消息或错误。但是，在下面描述的情况下，即使在一天后程序仍在忙于运行，这对于简单表达式的问题来说似乎是不合理的。是否可以得到以下给定表达式的解（对于q\u IN\u w）在运行此问题之前，我使用conda将Symp

我试图用符号来解一个常数系数的积分（

a_w

，

b_IN_w

，

c_IN_w

），它们由简单的代数表达式组成

我试过的我已经试着运行下面给出的代码一整天了，但是当我一天后检查它时，它仍然在运行。我以前使用过Symphy，我知道它可能无法解决一些复杂的操作，在这些操作中，它会抛出某种指示问题的消息或错误。但是，在下面描述的情况下，即使在一天后程序仍在忙于运行，这对于简单表达式的问题来说似乎是不合理的。是否可以得到以下给定表达式的解（对于

q\u IN\u w

）

在运行此问题之前，我使用conda将Symphy更新为其最新版本

import sympy as sym

def deg_to_rad(theta_deg):
    from numpy import pi
    theta_rad = (pi/180)*theta_deg
    return theta_rad


r, a_w, a_o, a_g, b_IN_w, b_IN_o, b_IN_g, c_IN_2w, c_IN_2o, c_IN_2g, r_1, r_2, R, \
sigma_dia, IFT_ow, theta_IN_CA_deg, D_IN_ads_coeff, nablaP, mu_w, deltaP = \
sym.symbols('r, a_w, a_o, a_g, b_IN_w, b_IN_o, b_IN_g, c_IN_2w, c_IN_2o, c_IN_2g, r_1, r_2, R, \
sigma_dia, IFT_ow, theta_IN_CA_deg, D_IN_ads_coeff, nablaP, mu_w, deltaP')

l_IN_slip = sigma_dia/((sym.pi - deg_to_rad(theta_IN_CA_deg))**4)
W_IN_egy = IFT_ow*(1 + sym.cos(deg_to_rad(theta_IN_CA_deg)))
u_IN_s = (l_IN_slip*R*nablaP)/(2*mu_w)
u_IN_ads = (D_IN_ads_coeff/W_IN_egy)*deltaP
u_IN_s_eff = (u_IN_s - u_IN_ads)

b_IN_g = 0
b_IN_o = 2*(a_g - a_o)*(r_1**2)
b_IN_w = b_IN_o + 2*(a_o - a_w)*(r_2**2)

c_IN_2w = u_IN_s_eff - a_w*(R**2) - b_IN_w*sym.log(R)

q_IN_w = sym.integrate((a_w*(r**2) + b_IN_w*(sym.log(r)) + c_IN_2w)*(2*sym.pi*r), (r, r_2, R))

在这里，它在30秒内完成，但它的答案是

sym.pi

pi*R**4*a_w/2 - R**2*(pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**4 - 720*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**3 + 194400*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**2 - 23328000*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg + 1049760000*pi**4*D_IN_ads_coeff*deltaP*mu_w + pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**4 - 720*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**3 + 194400*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**2 - 23328000*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg + 1049760000*pi**4*IFT_ow*R**2*a_w*mu_w*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*R**2*a_w*mu_w - 524880000*IFT_ow*R*nablaP*sigma_dia*cos(pi*theta_IN_CA_deg/180) - 524880000*IFT_ow*R*nablaP*sigma_dia + 2*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) + 2*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*log(R) + pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4 - 1440*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) - 1440*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*log(R) - 720*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3 + 388800*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 388800*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*log(R) + 194400*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2 - 46656000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) - 46656000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*log(R) - 23328000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg + 2099520000*pi**4*IFT_ow*a_g*mu_w*r_1**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 2099520000*pi**4*IFT_ow*a_g*mu_w*r_1**2*log(R) + 1049760000*pi**4*IFT_ow*a_g*mu_w*r_1**2*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*a_g*mu_w*r_1**2 - 2*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) - 2*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*log(R) - pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) - pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4 + 1440*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) + 1440*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*log(R) + 720*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) + 720*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3 - 388800*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 388800*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*log(R) - 194400*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) - 194400*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2 + 46656000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) + 46656000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*log(R) + 23328000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) + 23328000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg - 2099520000*pi**4*IFT_ow*a_o*mu_w*r_1**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 2099520000*pi**4*IFT_ow*a_o*mu_w*r_1**2*log(R) - 1049760000*pi**4*IFT_ow*a_o*mu_w*r_1**2*cos(pi*theta_IN_CA_deg/180) - 1049760000*pi**4*IFT_ow*a_o*mu_w*r_1**2 + 2*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) + 2*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*log(R) + pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4 - 1440*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) - 1440*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*log(R) - 720*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3 + 388800*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 388800*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*log(R) + 194400*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2 - 46656000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) - 46656000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*log(R) - 23328000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg + 2099520000*pi**4*IFT_ow*a_o*mu_w*r_2**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 2099520000*pi**4*IFT_ow*a_o*mu_w*r_2**2*log(R) + 1049760000*pi**4*IFT_ow*a_o*mu_w*r_2**2*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*a_o*mu_w*r_2**2 - 2*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) - 2*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*log(R) - pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) - pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4 + 1440*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) + 1440*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*log(R) + 720*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) + 720*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3 - 388800*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 388800*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*log(R) - 194400*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) - 194400*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2 + 46656000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) + 46656000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*log(R) + 23328000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) + 23328000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg - 2099520000*pi**4*IFT_ow*a_w*mu_w*r_2**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 2099520000*pi**4*IFT_ow*a_w*mu_w*r_2**2*log(R) - 1049760000*pi**4*IFT_ow*a_w*mu_w*r_2**2*cos(pi*theta_IN_CA_deg/180) - 1049760000*pi**4*IFT_ow*a_w*mu_w*r_2**2)/(pi**3*IFT_ow*mu_w*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**3*IFT_ow*mu_w*theta_IN_CA_deg**4 - 720*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**3 + 194400*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**2 - 23328000*pi**3*IFT_ow*mu_w*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**3*IFT_ow*mu_w*theta_IN_CA_deg + 1049760000*pi**3*IFT_ow*mu_w*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**3*IFT_ow*mu_w) - pi*a_w*r_2**4/2 + r_2**2*(pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**4 - 720*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**3 + 194400*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg**2 - 23328000*pi**4*D_IN_ads_coeff*deltaP*mu_w*theta_IN_CA_deg + 1049760000*pi**4*D_IN_ads_coeff*deltaP*mu_w + pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**4 - 720*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**3 + 194400*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg**2 - 23328000*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*R**2*a_w*mu_w*theta_IN_CA_deg + 1049760000*pi**4*IFT_ow*R**2*a_w*mu_w*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*R**2*a_w*mu_w - 524880000*IFT_ow*R*nablaP*sigma_dia*cos(pi*theta_IN_CA_deg/180) - 524880000*IFT_ow*R*nablaP*sigma_dia + 2*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) + 2*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*log(R) + pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**4 - 1440*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) - 1440*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*log(R) - 720*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**3 + 388800*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 388800*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*log(R) + 194400*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg**2 - 46656000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) - 46656000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*log(R) - 23328000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*a_g*mu_w*r_1**2*theta_IN_CA_deg + 2099520000*pi**4*IFT_ow*a_g*mu_w*r_1**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 2099520000*pi**4*IFT_ow*a_g*mu_w*r_1**2*log(R) + 1049760000*pi**4*IFT_ow*a_g*mu_w*r_1**2*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*a_g*mu_w*r_1**2 - 2*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) - 2*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*log(R) - pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) - pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**4 + 1440*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) + 1440*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*log(R) + 720*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) + 720*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**3 - 388800*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 388800*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*log(R) - 194400*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) - 194400*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg**2 + 46656000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) + 46656000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*log(R) + 23328000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) + 23328000*pi**4*IFT_ow*a_o*mu_w*r_1**2*theta_IN_CA_deg - 2099520000*pi**4*IFT_ow*a_o*mu_w*r_1**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 2099520000*pi**4*IFT_ow*a_o*mu_w*r_1**2*log(R) - 1049760000*pi**4*IFT_ow*a_o*mu_w*r_1**2*cos(pi*theta_IN_CA_deg/180) - 1049760000*pi**4*IFT_ow*a_o*mu_w*r_1**2 + 2*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) + 2*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*log(R) + pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**4 - 1440*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) - 1440*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*log(R) - 720*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**3 + 388800*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 388800*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*log(R) + 194400*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg**2 - 46656000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) - 46656000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*log(R) - 23328000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**4*IFT_ow*a_o*mu_w*r_2**2*theta_IN_CA_deg + 2099520000*pi**4*IFT_ow*a_o*mu_w*r_2**2*log(R)*cos(pi*theta_IN_CA_deg/180) + 2099520000*pi**4*IFT_ow*a_o*mu_w*r_2**2*log(R) + 1049760000*pi**4*IFT_ow*a_o*mu_w*r_2**2*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**4*IFT_ow*a_o*mu_w*r_2**2 - 2*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*log(R)*cos(pi*theta_IN_CA_deg/180) - 2*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*log(R) - pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) - pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**4 + 1440*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*log(R)*cos(pi*theta_IN_CA_deg/180) + 1440*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*log(R) + 720*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) + 720*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**3 - 388800*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 388800*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*log(R) - 194400*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) - 194400*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg**2 + 46656000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*log(R)*cos(pi*theta_IN_CA_deg/180) + 46656000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*log(R) + 23328000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) + 23328000*pi**4*IFT_ow*a_w*mu_w*r_2**2*theta_IN_CA_deg - 2099520000*pi**4*IFT_ow*a_w*mu_w*r_2**2*log(R)*cos(pi*theta_IN_CA_deg/180) - 2099520000*pi**4*IFT_ow*a_w*mu_w*r_2**2*log(R) - 1049760000*pi**4*IFT_ow*a_w*mu_w*r_2**2*cos(pi*theta_IN_CA_deg/180) - 1049760000*pi**4*IFT_ow*a_w*mu_w*r_2**2)/(pi**3*IFT_ow*mu_w*theta_IN_CA_deg**4*cos(pi*theta_IN_CA_deg/180) + pi**3*IFT_ow*mu_w*theta_IN_CA_deg**4 - 720*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**3*cos(pi*theta_IN_CA_deg/180) - 720*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**3 + 194400*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**2*cos(pi*theta_IN_CA_deg/180) + 194400*pi**3*IFT_ow*mu_w*theta_IN_CA_deg**2 - 23328000*pi**3*IFT_ow*mu_w*theta_IN_CA_deg*cos(pi*theta_IN_CA_deg/180) - 23328000*pi**3*IFT_ow*mu_w*theta_IN_CA_deg + 1049760000*pi**3*IFT_ow*mu_w*cos(pi*theta_IN_CA_deg/180) + 1049760000*pi**3*IFT_ow*mu_w) + (2*pi*R**2*a_g*r_1**2 - 2*pi*R**2*a_o*r_1**2 + 2*pi*R**2*a_o*r_2**2 - 2*pi*R**2*a_w*r_2**2)*log(R) - (2*pi*a_g*r_1**2*r_2**2 - 2*pi*a_o*r_1**2*r_2**2 + 2*pi*a_o*r_2**4 - 2*pi*a_w*r_2**4)*log(r_2)

您没有在

deg___rad

中定义

pi

，我想应该是

sym.pi

对不起。

pi

实际上是numpy import pi中定义为

的数值，但我不记得在这里的代码中包含它。如果在函数deg__rad
中使用pi
而不是在代码中实现的pi
的数值，则您是正确的（根据您的评论进行编辑），然后它会在合理的时间内运行。但是，如果使用问题中所示的pi
的数值，那么代码就会无限期地运行。你能解释一下为什么会发生这种情况吗？不知道。我不是Symphy
专家。也许这与简化有关，因为Symphy知道pi
，而不是它是数值。也许可以检查Symphy bugtracker。如果它没有列出，请用一个最简单的例子打开一个新问题，其中发生了这种情况。我猜Symphy不能很好地使用numpy值，尽管它（pi
）作为一个浮点，因为scipy使用自己的约定来存储变量和常量。SymPy通常会更好地处理符号值。例如，SymPy.sin（pi）
将计算为0
，如果pi
是符号，但如果它是数字，它将只计算为接近但不完全等于0的值（从技术上讲，numpy.pi
不完全是π，所以这是正确的）。这种情况有很多。@asmuer感谢您的澄清！