Moderator: Rhiannon
Code: Selecteer alles
from sympy import *
g, R, c_zuiger, cv_edelgas, cv_lucht, A, Vboven_f, V0, d_zuiger, Vonder_f, n1, n2, T1, T2, Tzuiger, p1, p2, U1, U2, U1_f, U2_f, T_f, m, ponder_f, pboven_f = symbols("g, R, c_zuiger, cv_edelgas, cv_lucht, A, Vboven_f, V0, d_zuiger, Vonder_f, n1, n2, T1, T2, Tzuiger, p1, p2, U1, U2, U1_f, U2_f, T_f, m, ponder_f, pboven_f")
eq1 = Eq(g, 9.81)
eq2 = Eq(R, 8.31) # J/(mol*K)
eq3 = Eq(c_zuiger, 210) # J/(kgK).
eq4 = Eq(cv_edelgas, 1.5*R) # J/(mol*K)
eq5 = Eq(cv_lucht, 20.77) # J/(mol*K)
eq6 = Eq(m, 0.5) # kg
eq7 = Eq(A, 100*(10**(-4))) # m**2
eq8 = Eq(Vboven_f, V0-d_zuiger*A) # m**3
eq9 = Eq(Vonder_f, V0+d_zuiger*A) # m**3
eq10 = Eq(n1 , 0.05) # mol
eq11 = Eq(n2 , 0.03) # mol
eq12 = Eq(T1 , 273-90) # K
eq13 = Eq(T2 , 273+46) # K
eq14 = Eq(Tzuiger , 273+100) # K
eq15 = Eq(p1*V0 , n1*R*T1) # J
eq16 = Eq(p2*V0 , n2*R*T2) # J
eq17 = Eq(p2 , p1 + m*g/A) # Pa
eq18 = Eq(pboven_f*Vboven_f , n1*R*T_f)
eq19 = Eq(ponder_f*Vonder_f , n2*R*T_f)
eq20 = Eq(ponder_f , pboven_f + m*g/A)
eq21 = Eq(U1 , n1*cv_edelgas*T1)
eq22 = Eq(U2 , n2*cv_lucht*T2)
eq23 = Eq(U1_f , n1*cv_edelgas*T_f)
eq24 = Eq(U2_f , n2*cv_lucht*T_f)
eq25 = Eq(U1+U2+m*c_zuiger*Tzuiger , U1_f+U2_f+m*c_zuiger*T_f+m*g*d_zuiger)
solve([eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9, eq10, eq11, eq12, eq13, eq14, eq15, eq16, eq17, eq18, eq19, eq20, eq21, eq22, eq23, eq24, eq25], dict=True)
