ewdy / structural-design Goto Github PK

# create beam array and load with position, shear, moment and deflection values
beam = np.zeros((4, SAMPLE_SIZE))
beam[0, : ] = np.linspace(0, L.magnitude, SAMPLE_SIZE) # position
for i in range(SAMPLE_SIZE):
    beam[1, i] = shear_x(beam[0, i]*L.units, P, a, L).magnitude # stuff shear values
    beam[2, i] = moment_x(beam[0, i]*L.units, P, a, L).magnitude # stuff moment values
    beam[3, i] = deflection_x(beam[0, i]*L.units, P, a, L, EI).magnitude # stuff deflection values
# beam[0:]
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-8-3cac110d7e7a> in <module>()
      3 beam[0, : ] = np.linspace(0, L.magnitude, SAMPLE_SIZE) # position
      4 for i in range(SAMPLE_SIZE):
----> 5     beam[1, i] = shear_x(beam[0, i]*L.units, P, a, L).magnitude # stuff shear values
      6     beam[2, i] = moment_x(beam[0, i]*L.units, P, a, L).magnitude # stuff moment values
      7     beam[3, i] = deflection_x(beam[0, i]*L.units, P, a, L, EI).magnitude # stuff deflection values
<ipython-input-5-787cab1916c6> in shear_x(x, P, a, L)
     15     """
     16     b = L - a # derive b given a and L
---> 17     if x >= 0 and x <= a:
     18         V_x = (P*b)/L
     19     elif x > a and x <= L:
/usr/local/lib/python3.6/dist-packages/pint/quantity.py in <lambda>(self, other)
   1077     __lt__ = lambda self, other: self.compare(other, op=operator.lt)
   1078     __le__ = lambda self, other: self.compare(other, op=operator.le)
-> 1079     __ge__ = lambda self, other: self.compare(other, op=operator.ge)
   1080     __gt__ = lambda self, other: self.compare(other, op=operator.gt)
   1081 
/usr/local/lib/python3.6/dist-packages/pint/quantity.py in compare(self, other, op)
   1065                 return op(self._convert_magnitude_not_inplace(UnitsContainer()), other)
   1066             else:
-> 1067                 raise ValueError('Cannot compare Quantity and {0}'.format(type(other)))
   1068 
   1069         if self._units == other._units:
ValueError: Cannot compare Quantity and <class 'int'>