1. import numpy as np
2. import matplotlib.pyplot as plt
3.  
4. # Veriler
5. F0 = 4787.875  # N
6. F1 = 57454.5 # N
7. P0 = 476000  # Pa
8. P1 = 5240000  # Pa
9. V0 = 1.2441e-3  # m³
10. V1 = 1.131e-4  # m³
11.  
12. # İzotermal sabit
13. iso_constant = P0 * V0
14.  
15. # Politropik sabit
16. poly_constant = P0 * (V0 ** 1.35)
17.  
18. # Piston alanı (m²)
19. A = 6.94e-3
20.  
21. # Vuruş (stroke) değerleri (m)
22. strokes = np.linspace(0, 0.163, 100)
23.  
24. # Hacim değerlerini hesaplayalım (m³)
25. volumes = V0 + strokes * (V1 - V0) / 0.163
26.  
27. # İzotermal basınç değerlerini hesaplayalım (Pa)
28. iso_pressures = iso_constant / volumes
29.  
30. # Politropik basınç değerlerini hesaplayalım (Pa)
31. poly_pressures = poly_constant / (volumes ** 1.35)
32.  
33. # Yük (force) değerlerini hesaplayalım (N)
34. iso_loads = iso_pressures * A
35. poly_loads = poly_pressures * A
36.  
37.  
38. # Grafik oluşturma
39. plt.plot(strokes, iso_loads, marker='o', label='Isothermal')
40. plt.plot(strokes, poly_loads, marker='x', label='Polytropic')
41. plt.title('Load-Stroke Graph')
42. plt.xlabel('Stroke (m)')
43. plt.ylabel('Load (N)')
44. plt.legend()
45. plt.grid(True)
46. plt.show()
47.