The finite element (FE) method was used to calculate the axial and radial stress distributions as a function of axial distance, z, from the centre, and radius, r, in an elastic fibre surrounded by a plastic matrix. Plastic deformation of the matrix was considered to exert a uniform interfacial stress, τ, along half the length of the fibre. Axisymmetric models were created for uniform cylindrical, ellipsoidal, paraboloidal and conical fibres characterised by an axial ratio, q, and half length, L. Young's modulus for the material of the fibre and L were arbitrarily assigned values of unity, since they act as scaling factors; q also acts as a scaling factor but was assigned a value of 10 to create models with a fibrous appearance. For the cylindrical fibre, the axial stress increased linearly from the end towards the centre; the radial stress was more evenly distributed. At the other extreme, the conical fibre showed a uniform distribution of axial and radial stress. Results for ellipsoidal and paraboloidal fibres were intermediate between these two extremes. In general, the effect of taper is to lower peak stress at the fibre centre and make the stress distribution throughout the fibre more even. These results are in good agreement with recent analytical theories for the axial distribution of surface radial stress and axial stress along the fibre axis. However, FE models have the advantage of predicting full three-dimensional stress distributions.