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