The article describes a shock safety modeling method for low-voltage electric devices, based on using a Bayesian network. This method allows for taking into account all possible combinations of the reliability and unreliability states for the shock protection elements under concern. The developed method allows for investigating electric shock incidents, analysing and assessing shock risks, as well as for determining criteria of dimensioning shock protection means, also with respect to reliability of the particular shock protection elements. Dependencies for determining and analysing the probability of appearance of reliability states of protection as well as an electric shock risk are presented in the article.