“…The normalized relative displacements of the two springs in the directions normal and tangential to the boundary line (i.e., δ n / h and δ s / h , where h is the characteristic length of a concrete spring element [Figure 1a]) represent the tensile/compressive and shear strains of the concrete, respectively. The normal, shear and rotational stiffnesses are set to approximate the elastic properties of concrete at the continuum levels and the rotational stiffness K φ is suggested to be zero when either the normal or the shear spring reaches its maximum strength (Saito and Hikosaka, 1999). Through the virtual work principle, the stiffness matrix associated with the two‐particle assembly, describing the relationship between the force vector ( X 1 , Y 1 , M 1 , X 2 , Y 2 , M 2 ) and the displacement vector ( u 1 ,v 1 ,θ 1 ,u 2 ,v 2 ,θ 2 ) can be obtained (Bolander and Saito, 1998).…”