“…In order to couple the ubiquitiform characteristic of pore network with the strain softening model of quasi-brittle materials, the energy balance principle is employed again. According to references (Carpinteri et al., 2002, 2003, 2009; Ou et al., 2019), in the damage band, the required energy to generate new pore surface is equal to the area under the stress-elongation curve in an interval with length of Δ w j . Therefore, in the j -th softening subprocess, the energy balance equation is written as where A is the cross-sectional area of the PUM, i.e., A = δitalicmax2; trueσ¯j is the stress of the material in the j -th softening subprocess; Δ w j is the increment of elongation in the j -th softening subprocess at the stress trueσ¯j; γf is the fracture surface energy equal to half of the material fracture energy G f; Δ S pj is the incremental pore surface area in the j -th subprocess, i.e., Δ S pj = S pj − S pj -1 .…”