Sophisticated analytical solution methods for the stress concentration problem in multilayered composites with interference-fit bolts have been developed on the basis of layer-related solutions. The model is based on a circular or elliptical isotropic bolt integrated in a fibre- or textile-reinforced multilayered composite plate. For the simulation of the interference between bolt and composite plate a modelling strategy has been proposed featuring the application of a defined temperature load to a metallic elastic inclusion. For analysing the stress concentration effects in the surrounding plate, complex-valued displacement functions in combination with the method of conformal mappings for the interface between inclusion and plate and a combination of boundary collocation and least squares method for the outer boundary are used. For the verification of the developed calculation methods, a number of experimental and numerical studies have been carried out and the decay behaviour of the distortions have been compared for different radians. For all combinations of multilayered composite plate and interference-fit bolt investigated so far, a good correlation of the analytically calculated and the numerically determined results can be observed.