diff --git a/model-optimizer/mo/utils/broadcasting.py b/model-optimizer/mo/utils/broadcasting.py
index c1d251a98c7356..d900514f033e0b 100644
--- a/model-optimizer/mo/utils/broadcasting.py
+++ b/model-optimizer/mo/utils/broadcasting.py
@@ -134,7 +134,10 @@ def explicit_broadcasting(input_value: np.array, target_shape: np.array, axes_ma
     :return: broadcasted value
     """
     res_shape, normalized_axes_mapping = explicit_shape_broadcasting(input_value.shape, target_shape, axes_mapping)
+    #TODO: Function 'expand_dims' should be replaced with 'numpy.expand_dims' if numpy version will be >=18.x in requirements.
     expand_dim_axis = set(np.arange(len(target_shape))) - set(normalized_axes_mapping)
-
-    input_expanded = np.expand_dims(input_value.copy(), axis=list(expand_dim_axis))
+    input_expanded = input_value.copy()
+    
+    for axis in sorted(list(expand_dim_axis)):
+        input_expanded = np.expand_dims(input_expanded, axis)
     return np.broadcast_to(input_expanded, res_shape)