1 def get_dx(self, component): 2 return [x - self.centroid[0] for x in component.x_start] 3 4 5 def get_dz(self, component): 6 return [x - self.centroid[1] for x in component.x_start] 7 8 9 def get_dP(self, xDist, zDist, V_x, V_z, area): 10 I_x = self.I_['x'] 11 I_z = self.I_['z'] 12 I_xz = self.I_['xz'] 13 denom = float(I_x * I_z - I_xz**2) 14 z = float() 15 for _ in range(len(xDist)): 16 z += float(-area * xDist[_] * (I_x * V_x - I_xz * V_z) / denom - 17 area * zDist[_] * (I_z * V_z - I_xz * V_x) / denom) 18 return z 19