Moderator: Rhiannon
Code: Selecteer alles
from sympy import *
R, P, t, T0, T, p0, p, cv, A1, A2, l1, l2, x, V0, V, n = symbols('R, P, t, T0, T, p0, p, cv, A1, A2, l1, l2, x, V0, V, n')
eq0 = Eq(R, 8.31)
eq1 = Eq(P, 75)
eq2 = Eq(t, 4*60)
eq3 = Eq(T0, 310)
eq4 = Eq(p0, 10**5)
eq5 = Eq(cv, 20.77) #J/molK
eq6 = Eq(A1, 40*10**-2)
eq7 = Eq(A2, 70*10**-2)
eq8 = Eq(l1, 0.25)
eq9 = Eq(l2, 0.25)
eq10 = Eq(V0, A1*l1 + A2*l2)
eq11 = Eq(V, A1*(l1-x) + A2*(l2+x))
eq12 = Eq(n, p0*V0/(R*T0))
eq13 = Eq(n, p*V/(R*T))
eq14 = Eq(P*t, x*(A2-A1)*p0+cv*n*(T-T0))
solve([eq0, eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8, eq9, eq10, eq11, eq12, eq13, eq14], dict=True)
volgens mij zie je dan over het hoofd dat de volume toe of afname als functie van de verplaatsing van de zuigers met staaf ertussen voor de kleine zuiger dus kleiner is als voor de grote zuiger. het systeem kan dus altijd een nieuwe positie vinden waarbij het volume verandert is en hoort bij de temperatuur en druk. Dus dan is druk weer overal gelijk en p0 dus staat er geen kracht op elke zuiger, dus ook geen kracht op de staaf.
Dat eerste snap ik best wel. Ik verwijs naar Eq11 uit mijn uitwerkingHansH schreef: ↑zo 23 jun 2024, 10:12volgens mij zie je dan over het hoofd dat de volume toe of afname als functie van de verplaatsing van de zuigers met staaf ertussen voor de kleine zuiger dus kleiner is als voor de grote zuiger. het systeem kan dus altijd een nieuwe positie vinden waarbij het volume verandert is en hoort bij de temperatuur en druk. Dus dan is druk weer overal gelijk en p0 dus staat er geen kracht op elke zuiger, dus ook geen kracht op de staaf.
Code: Selecteer alles
eq11 = Eq(V, A1*(l1-x) + A2*(l2+x))
Dat dit een oplossing is, snap ik. Met mijn uitwerking kom ik een vergelijking te weinig uit. Deze oplossing wordt als ik geen rekenfouten heb gemaakt ook afgedekt in mijn uitwerking. Maar waarom is dit dé oplossing? Andere oplossing zijn vanuit evenwichtstandpunt ook mogelijk lijkt mij. Heb je dan ook nog de tweede hoofdwet nodig om aan te tonen dat dit dé oplossing is?Xilvo schreef: ↑zo 23 jun 2024, 10:28 Het volume is 0,275 m3.
Dan vind ik bij p0=105 Pa en T=310 K n=10,67.
Met cv=20,76 J/mol.K en E=18000 J kom ik op een temperatuurstijging van 81,24 K, dat is voordat de zuigers gaan verplaatsen.
Dan wordt de nieuwe druk 1,262E5 Pa.
Met γ=cp/cv=1,4 en C=p.Vγ wordt het nieuwe volume bij p0 0,3247 m3
Dat geeft een zuigerverplaatsing van 0,166 m. Mooi in overeenstemming met wat ukster vindt.
Als de druk binnen hoger is dan de druk buiten, dan is de kracht naar rechts op de rechtse zuiger groter dan de kracht naar links op de linkse zuiger. Dus bewegen de zuigers naar links tot de druk binnen gelijk is aan de druk buiten.
Dat klopt bij nader inzien. Op basis daarvan kan ik p=p0 als vergelijking toevoegen aan mijn stelsel.HansH schreef: ↑zo 23 jun 2024, 11:43 als er een drukverschil is tussen binnen en buiten dan is er een kracht op zowel de linker als rechter zuiger die tegengesteld is maar vanwege het grotere oppervlak van de rechter zuiger is die kracht netto altijd naar rechts (links) als de druk binnen groter (kleiner) is dan buiten. dus zal de constructie dan net zo lang naar rechts (links) bewegen tot er geen netto kracht meer is. Dus dan moet het drukverschil ook 0 zijn.