“…where K is the spring stiffness coefficient, u i the unit vector, C the body center and a coefficient in (0, 1), Δl tot the volume differences per area unit and Δl i the length changes of spring ith (Zhengdong and Shuyuan, 2011). Vassilev (2000Vassilev ( , 2001, Provot (1995), Louchet et al (1995), , Vassilev and Spanlang (2000), Bayraktar et al (2007), Connor and Stevens (2003), Dochev and Vassilev (2003), Zhibin and Zhanli (2006), Haiyan and Zhaofeng (2008), Qing et al (2010), Liu et al (2012a, b), Fuhrmann et al (2003), Mesit et al (2007), Cho et al (2010), Rusinko and Swan (2012), Huang et al (2014), Zhang and Yuen (2001), Kang et al (2000), Smith (2011), Goldenthal (2010), Yang and Shang (2013)…”