This article tackles the problem of controlling articulated soft robots (ASRs), i.e., robots with either fixed or variable elasticity lumped at the joints. Classic control schemes rely on high-authority feedback actions, which have the drawback of altering the desired robot softness. The problem of accurate control of ASRs, without altering their inherent stiffness, is particularly challenging because of their complex and hard-to-model nonlinear dynamics. Leveraging a learned anticipatory action, iterative learning control (ILC) strategies do not suffer from these issues. Recently, ILC was adopted to perform position control of ASRs. However, the limitation of position-based ILC in controlling variable stiffness robots is that whenever the robot stiffness profile is changed, a different input action has to be learned. Our first contribution is to identify a wide class of ASRs, whose motion and stiffness adjusting dynamics can be proved to be decoupled. This class is described by two properties that we define: strong elastic coupling, relative to motors and links of the system and their connections; and homogeneity, relative to the characteristics of the motors. Furthermore, we design a torque-based ILC scheme that, starting from a rough estimation of the system parameters, refines the torque needed for the joint positions tracking. The resulting control scheme requires minimum knowledge of the system. Experiments on variable stiffness robots prove that the method effectively generalizes the iterative procedure with respect to the desired stiffness profile and allows good tracking performance. Finally, potential restrictions of the method, e.g., caused by friction phenomena, are discussed.