An investigation of the stability properties of stellar models describing massive stars is motivated observationally by the necessity to explain the observed Humphreys - Davidson (HD) limit and the variability of the most massive stars known, i.e. the existence of luminous blue variables (LBVs). Theoretically, a determination of the upper mass limit for stable stellar objects together with its physical explanation and interpretation is of fundamental interest.