code: def test_polar_moment_Beam3D():
      l, E, G, A, I1, I2 = symbols("l, E, G, A, I1, I2")
      I = [I1, I2]
  
      b = Beam3D(l, E, G, I, A)
      assert b.polar_moment() == I1 + I2
  
imports: from sympy.core.function import expand
  from sympy.core.numbers import Rational, pi
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol, symbols
  from sympy.sets.sets import Interval
  from sympy.simplify.simplify import simplify
  from sympy.physics.continuum_mechanics.beam import Beam
  from sympy.functions import SingularityFunction, Piecewise, meijerg, Abs, log
  from sympy.testing.pytest import raises
  from sympy.physics.units import meter, newton, kilo, giga, milli
  from sympy.physics.continuum_mechanics.beam import Beam3D
  from sympy.geometry import Circle, Polygon, Point2D, Triangle
  from sympy.core.sympify import sympify
label_desc: Use image node_id 14 for calling a global function with example usage: test_polar_moment_Beam3D() without return types