# Applying Stress Boundary Conditions in Commercial Finite Element Analysis Codes

I am trying to replicate a finite element analysis given in a research paper titled On the Detection of Stress Singularities in Finite Element Analysis 1 by G.B.Sinclair et. al.

The geometry of the problem is given in Figure 1.

and boundary conditions are $$(x, y) = (D/2 < x < L + D/2) \hspace{0.5in} \sigma_y = -p, \tau_{xy} = 0$$ and the left-most bolt hole is clamped.

Typically, in finite element analysis codes such as Patran/MSC Nastran, ANSYS, and FEMAP, we can apply displacement, force, and distributed loads as boundary conditions. The boundary conditions given in this problem specify stresses. How can I apply them in commercial FEA codes?