Discrete element methods (DEMs) provide new numerical means to study the behavior of soil-inclusion systems. In some cases, however, the classic DEM fails to model specific aspects of the inclusions. That is why a model based on spar elements is introduced, designed specifically for inclusions. In this model, the movement of the inclusion is considered as a dynamic process and is computed step by step in the same way as in the DEM. The model can be coupled with a DEM code, thus enabling one to simulate the interaction between an inclusion and a disk assembly. Contact laws at the contacts between disks and spar elements describe the interface constitutive behavior. Finally, the results obtained by simulating a geosynthetic anchorage in two different ways are reported. In the first case the inclusion is represented by disks, while in the last case it is represented by spar elements. The comparison shows that spar elements are much more versatile and can simplify the calibration of the discrete models used to simulate soil-inclusion systems. CE Database subject headings: Anchorage; Composite structures; Discrete elements; Soils; Geosynthetics; Granular materials; Interfaces.