notebook.community

Edit and run 





In [1]:



    


# Set up Physique for physical constants
import Physique

Useful fundamental constants and units conversions





In [2]:



    


from Physique import FundamentalPhysicalConstants as FPC



    











In [3]:



    


g_0pd = FPC[FPC["Quantity"].str.contains("gravity") ] 
# standard acceleration of gravity as a panda DataFrame
print(g_0pd)



    


















    






                             Quantity    Value Uncertainty    Unit
303  standard acceleration of gravity  9.80665        None  m s^-2

















In [4]:



    


g_0 = g_0pd["Value"].values[0]
print(type(g_0))
print(g_0)



    


















    






<class 'decimal.Decimal'>
9.80665

















In [5]:



    


# access the values you're interested in 
print(g_0pd.Quantity)
print(g_0pd.Value.get_values()[0]) 
print(g_0pd.Unit.get_values()[0])



    


















    






303    standard acceleration of gravity
Name: Quantity, dtype: object
9.80665
m s^-2

















In [6]:



    


from Physique import Conversions



    











In [7]:



    


Conversions[Conversions['Toconvertfrom'].str.contains("pound-force ")]



    


















    
Out[7]:











  
    

      
      Toconvertfrom
      to
      Multiplyby
    

  
  
    

      175
      foot pound-force (ft · lbf)
      joule (J)
      1.355818
    

    

      176
      foot pound-force per hour (ft · lbf/h)
      watt (W)
      0.0003766161
    

    

      177
      foot pound-force per minute (ft · lbf/min)
      watt (W)
      0.02259697
    

    

      178
      foot pound-force per second (ft · lbf/s)
      watt (W)
      1.355818
    

    

      340
      pound-force (lbf) 23
      newton (N)
      4.448222
    

    

      341
      pound-force foot (lbf · ft)
      newton meter (N · m)
      1.355818
    

    

      342
      pound-force foot per inch (lbf · ft/in)
      newton meter per meter (N · m/m)
      53.37866
    

    

      343
      pound-force inch (lbf · in)
      newton meter (N · m)
      0.1129848
    

    

      344
      pound-force inch per inch (lbf · in/in)
      newton meter per meter (N · m/m)
      4.448222
    

    

      345
      pound-force per foot (lbf/ft)
      newton per meter (N/m)
      14.59390
    

    

      346
      pound-force per inch (lbf/in)
      newton per meter (N/m)
      175.1268
    

    

      347
      pound-force per pound (lbf/lb) (thrust to mass...
      newton per kilogram (N/kg)
      9.80665
    

    

      348
      pound-force per square foot (lbf/ft2)
      pascal (Pa)
      47.88026
    

    

      349
      pound-force per square inch (psi) (lbf/in2)
      pascal (Pa)
      6894.757
    

    

      350
      pound-force per square inch (psi) (lbf/in2)
      kilopascal (kPa)
      6.894757
    

    

      351
      pound-force second per square foot (lbf · s/ft2)
      pascal second (Pa · s)
      47.88026
    

    

      352
      pound-force second per square inch (lbf · s/in2)
      pascal second (Pa · s)
      6894.757
    

    

      372
      psi (pound-force per square inch) (lbf/in2)
      pascal (Pa)
      6894.757
    

    

      373
      psi (pound-force per square inch) (lbf/in2)
      kilopascal (kPa)
      6.894757
    

  




















In [8]:



    


Conversions[Conversions['to'].str.contains("newton ")]



    


















    
Out[8]:











  
    

      
      Toconvertfrom
      to
      Multiplyby
    

  
  
    

      137
      dyne (dyn)
      newton (N)
      0.000010
    

    

      138
      dyne centimeter (dyn · cm)
      newton meter (N · m)
      1.0E-7
    

    

      238
      kilogram-force (kgf)
      newton (N)
      9.80665
    

    

      239
      kilogram-force meter (kgf · m)
      newton meter (N · m)
      9.80665
    

    

      247
      kilopond (kilogram-force) (kp)
      newton (N)
      9.80665
    

    

      250
      kip (1 kip= 1000 lbf)
      newton (N)
      4448.222
    

    

      251
      kip (1 kip= 1000 lbf)
      kilonewton (kN)
      4.448222
    

    

      300
      ounce (avoirdupois)-force (ozf)
      newton (N)
      0.2780139
    

    

      301
      ounce (avoirdupois)-force inch (ozf · in)
      newton meter (N · m)
      0.007061552
    

    

      302
      ounce (avoirdupois)-force inch (ozf · in)
      millinewton meter (mN · m)
      7.061552
    

    

      336
      poundal
      newton (N)
      0.1382550
    

    

      340
      pound-force (lbf) 23
      newton (N)
      4.448222
    

    

      341
      pound-force foot (lbf · ft)
      newton meter (N · m)
      1.355818
    

    

      342
      pound-force foot per inch (lbf · ft/in)
      newton meter per meter (N · m/m)
      53.37866
    

    

      343
      pound-force inch (lbf · in)
      newton meter (N · m)
      0.1129848
    

    

      344
      pound-force inch per inch (lbf · in/in)
      newton meter per meter (N · m/m)
      4.448222
    

    

      345
      pound-force per foot (lbf/ft)
      newton per meter (N/m)
      14.59390
    

    

      346
      pound-force per inch (lbf/in)
      newton per meter (N/m)
      175.1268
    

    

      347
      pound-force per pound (lbf/lb) (thrust to mass...
      newton per kilogram (N/kg)
      9.80665
    

    

      423
      ton-force (2000 lbf)
      newton (N)
      8896.443
    

    

      424
      ton-force (2000 lbf)
      kilonewton (kN)
      8.896443
    

  




















In [9]:



    


lbf2N = Conversions.loc[340,:];
print(lbf2N)



    


















    






Toconvertfrom    pound-force (lbf) 23
to                         newton (N)
Multiplyby                   4.448222
Name: 340, dtype: object

















In [10]:



    


print(lbf2N.Toconvertfrom)
print(lbf2N.to)
print(lbf2N.Multiplyby)



    


















    






pound-force (lbf) 23
newton (N)
4.448222

kinematics conversions; length, speed, acceleration





In [33]:



    


Conversions[Conversions['Toconvertfrom'].str.contains("foot ")];



    











In [34]:



    


ft_to_m = Conversions.loc[340,:];
print(ft_to_m)



    


















    






Toconvertfrom    pound-force (lbf) 23
to                         newton (N)
Multiplyby                   4.448222
Name: 340, dtype: object

















In [38]:



    


Conversions[Conversions['Toconvertfrom'].str.contains("mile ")];



    











In [37]:



    


mi_to_km = Conversions.loc[271,:];
print(mi_to_km)



    


















    






Toconvertfrom         mile (mi)
to               kilometer (km)
Multiplyby             1.609344
Name: 271, dtype: object

















In [39]:



    


mph_to_m_per_s = Conversions.loc[277,:];
print(mph_to_m_per_s)



    


















    






Toconvertfrom      mile per hour (mi/h)
to               meter per second (m/s)
Multiplyby                      0.44704
Name: 277, dtype: object

















In [43]:



    


Conversions[Conversions['to'].str.contains("meter per ")];



    











In [42]:



    


ft_per_s2_to_m_per_s2 = Conversions.loc[173,:];
print(ft_per_s2_to_m_per_s2)



    


















    






Toconvertfrom    foot per second squared (ft/s2)
to               meter per second squared (m/s2)
Multiplyby                                0.3048
Name: 173, dtype: object

mass and "weight" conversions





In [57]:



    


Conversions[Conversions['to'].str.contains("kilogram") & Conversions['Toconvertfrom'].str.contains("pound")]



    


















    
Out[57]:











  
    

      
      Toconvertfrom
      to
      Multiplyby
    

  
  
    

      48
      British thermal unitIT per pound (BtuIT/lb)
      joule per kilogram (J/kg)
      2326
    

    

      49
      British thermal unitth per pound (Btuth/lb)
      joule per kilogram (J/kg)
      2324.444
    

    

      50
      British thermal unitIT per pound degree Fahren...
      joule per kilogram kelvin [J/(kg · K)]
      4186.8
    

    

      51
      British thermal unitth per pound degree Fahren...
      joule per kilogram kelvin [J/(kg · K)]
      4184
    

    

      52
      British thermal unitIT per pound degree Rankin...
      joule per kilogram kelvin [J/(kg · K)]
      4186.8
    

    

      53
      British thermal unitth per pound degree Rankin...
      joule per kilogram kelvin [J/(kg · K)]
      4184
    

    

      334
      pound (avoirdupois) (lb) 22
      kilogram (kg)
      0.4535924
    

    

      335
      pound (troy or apothecary) (lb)
      kilogram (kg)
      0.3732417
    

    

      339
      pound foot squared (lb · ft2)
      kilogram meter squared (kg · m2)
      0.04214011
    

    

      347
      pound-force per pound (lbf/lb) (thrust to mass...
      newton per kilogram (N/kg)
      9.80665
    

    

      353
      pound inch squared (lb · in2)
      kilogram meter squared (kg ·m2)
      0.0002926397
    

    

      354
      pound per cubic foot (lb/ft3)
      kilogram per cubic meter (kg/m3)
      16.01846
    

    

      355
      pound per cubic inch (lb/in3)
      kilogram per cubic meter (kg/m3)
      27679.90
    

    

      356
      pound per cubic yard (lb/yd3)
      kilogram per cubic meter (kg/m3)
      0.5932764
    

    

      357
      pound per foot (lb/ft)
      kilogram per meter (kg/m)
      1.488164
    

    

      360
      pound per gallon [Canadian and U.K. (Imperial)...
      kilogram per cubic meter (kg/m3)
      99.77637
    

    

      361
      pound per gallon [Canadian and U.K. (Imperial)...
      kilogram per liter (kg/L)
      0.09977637
    

    

      362
      pound per gallon (U.S.) (lb/gal)
      kilogram per cubic meter (kg/m3)
      119.8264
    

    

      363
      pound per gallon (U.S.) (lb/gal)
      kilogram per liter (kg/L)
      0.1198264
    

    

      364
      pound per horsepower hour [lb/(hp · h)]
      kilogram per joule (kg/J)
      1.689659E-7
    

    

      365
      pound per hour (lb/h)
      kilogram per second (kg/s)
      0.0001259979
    

    

      366
      pound per inch (lb/in)
      kilogram per meter (kg/m)
      17.85797
    

    

      367
      pound per minute (lb/min)
      kilogram per second (kg/s)
      0.007559873
    

    

      368
      pound per second (lb/s)
      kilogram per second (kg/s)
      0.4535924
    

    

      369
      pound per square foot (lb/ft2)
      kilogram per square meter (kg/m2)
      4.882428
    

    

      370
      pound per square inch (not pound-force) (lb/in2)
      kilogram per square meter (kg/m2)
      703.0696
    

    

      371
      pound per yard (lb/yd)
      kilogram per meter (kg/m)
      0.4960546
    

  




















In [68]:



    


lbm_to_kg = Conversions.loc[334,:];
print(lbm_to_kg)



    


















    






Toconvertfrom    pound (avoirdupois) (lb) 22
to                             kilogram (kg)
Multiplyby                         0.4535924
Name: 334, dtype: object

Specific Impulse





In [60]:



    


F_avg_thrust_USC = 73500 # lbf; USC = United States Customary units



    











In [61]:



    


I_s_RP1 = 353 # secs., max I_sp, https://en.wikipedia.org/wiki/RP-1



    











In [62]:



    


F_avg_thrust_SI = F_avg_thrust_USC * lbf2N.Multiplyby
print(F_avg_thrust_SI)



    


















    






326944.317000

















In [63]:



    


m_dot_propellant_expelled_avg_SI = F_avg_thrust_SI / (I_s_RP1 * g_0)
print(m_dot_propellant_expelled_avg_SI)



    


















    






94.44487840960206383628591970

















In [64]:



    


T_burn = 3 * 60



    











In [65]:



    


m_propellant_tot_SI = m_dot_propellant_expelled_avg_SI * T_burn
print(m_propellant_tot_SI)



    


















    






17000.07811372837149053146555

















In [66]:



    


F_avg_thrust_USC = 5000 # lbf; USC = United States Customary units
F_avg_thrust_SI = F_avg_thrust_USC * lbf2N.Multiplyby
print(F_avg_thrust_SI)
m_dot_propellant_expelled_avg_SI = F_avg_thrust_SI / (I_s_RP1 * g_0)
print(m_dot_propellant_expelled_avg_SI)
T_burn = 6 * 60
m_propellant_tot_SI = m_dot_propellant_expelled_avg_SI * T_burn
print(m_propellant_tot_SI)



    


















    






22241.110000
6.424821660517147199747341476
2312.935797786172991909042931

















In [69]:



    


M_0_lbm = 57000
M_0_SI = M_0_lbm * lbm_to_kg.Multiplyby
print(M_0_SI)



    


















    






25854.7668000

















In [70]:



    


25854 - (17000 + 2313)



    


















    
Out[70]:








6541

















In [71]:



    


6541 / 25854.



    


















    
Out[71]:








0.2529976019184652

















In [72]:



    


print(650 * lbm_to_kg.Multiplyby)
print(1100 * lbm_to_kg.Multiplyby)



    


















    






294.8350600
498.9516400

















In [73]:



    


(6541 - 295) / (25854 - 295)



    


















    
Out[73]:








0.24437575805000195

















In [74]:



    


(6541 - 500) / (25854 - 500)



    


















    
Out[74]:








0.23826615129762563

















In [75]:



    


6541 / (25854 - 295)



    


















    
Out[75]:








0.25591768066043274

















In [ ]:

Content source: ernestyalumni/Propulsion

Similar notebooks:

notebook.community | gallery | about