In this work, we develop and implement a plasticity model suitable to simulate the stress–strain response of 3D-printed metals subjected to multi-axial cyclic loads. The model is formulated within the infinitesimal strains framework and makes use of Hill’s orthotropic yield condition (Hill 1948), the associative flow rule, and Chaboche’s isotropic hardening rule (Chaboche 1979). An advanced kinematic hardening rule is developed based on a multi-component rule by (Dafalias and Feigenbaum 2011). The original rule is extended by an orthotropic factor that makes the kinematic hardening rate directionally dependent. The Hill’s yield function is employed to capture orthotropy in the yield strength expected in 3D-printed metals, while the advanced kinematic hardening rule is expected to be needed to accurately predict multi-axial ratcheting and material behavior under multi-axial cyclic plastic loading. The finite element implementation of the model mainly follows (Marek et al. 2015) and is described in detail, including the initial value problem formulation, explicit Euler integration, and yield surface return projection. All the implementation procedures are provided in general form suitable for any displacement-based finite element code, and an implementation example is provided for the commercial software Abaqus/Standard using the UMAT interface. Tensor operations are implemented using the Mandel notation, and matrix counterparts of these operations are shown. Also, the numerical stability is analyzed in detail by spectral decomposition of the model’s constitutive Jacobian. Several numerical examples are provided to verify the implementation, stability, and robustness of the code.