def calcSupportReac(self):
netLoad = 0
netBendingMoment = 0
integralOfBendingMoment = []
for x in self.shearForceEq:
temp = x.sub(self.length)
netLoad = netLoad + temp
for x in self.bendingMomentEq:
temp = x.integrate()
integralOfBendingMoment.append(temp)
for x in integralOfBendingMoment:
temp = x.sub(self.length)
netBendingMoment = temp + netBendingMoment
if (self.supportType == 'pin joint'):
self.getDiscreteForce(
self.support1, (netBendingMoment - netLoad * self.support1) *
-1 / (self.support2 - self.support1))
self.getDiscreteForce(
self.support2, (netLoad * self.support2 - netBendingMoment) *
-1 / (self.support2 - self.support1))
if (self.supportType == 'cantilever'):
self.getDiscreteForce(0, -1 * netLoad)
self.getBendingMoment(0, -1 * netBendingMoment)
netBendingMoment = 0
def test_integrate():
assert x.integrate(x) == x ** 2 / 2
assert x.integrate((x, 0, 1)) == S(1) / 2
def test_integrate():
assert x.integrate(x) == x**2 / 2
assert x.integrate((x, 0, 1)) == S(1) / 2