Arterial stenosis is a vascular pathology which leads to serious cardiovascular diseases. Blood flow through a constriction generates sound and vibration due to fluctuating turbulent pressures. Generated vibro-acoustic waves propagate through surrounding soft tissues and reach the skin surface and may provide valuable insight for noninvasive diagnostic purposes. Motivated by the aforementioned phenomena, vibration of constricted arteries is investigated employing computational models. The flow-induced pressure field in an artery is modeled as broadband harmonic pressure loading based on previous studies in the literature and applied on the inner artery wall. Harmonic analysis is performed for determining radial velocity responses on the outer surface of the models. Results indicate that stenosis severities higher than 70% lead to significant increase in response amplitudes, especially at high frequencies between 250 and 600 Hz. The findings agree well with experimental and theoretical results in the literature considering bending mode frequencies, amplitude scales, and mainly excited frequency ranges. It is seen that artery vibration is sensitive to the phase behavior of pressure loading but its effect becomes less significant with the presence of surrounding tissue. As the surrounding tissue thickness increases, radial velocity response amplitudes decrease but the effect of changes in tissue elastic modulus is more pronounced.