A new plastic–damage constitutive model based on the combination of damage mechanics and classical plastic theory was developed to simulate the failure of concrete. In order to explain different material behaviors of concrete under tensile and compressive loadings, the plastic yield criterion, the different kinematic hardening rule for tension and compressive and the isotropic flow rule were established in the effective stress space. Meanwhile, two different empirical damage evolution equations were adopted: one for compression and the other for tension. A multi-axial damage influence factor was also introduced to fully describe the anisotropic damage of concrete. Finally, the model response was compared with a wide range of experiment results. The results showed that the model could well describe the nonlinear behavior of concrete in a complex stress state.

Considering that ANSYS software will automatically quit or the computer will freeze when generating random aggregate models of concrete by using some existing methods that are based on the ANSYS parametric design language (APDL), a new method of random aggregate placement using the ESEL command in APDL and the rotation of the local coordinate system is proposed in this paper. According to this method, a multiscale macroscopic and mesoscopic finite element model of the No. 9 non-overflow dam section of Shimantan dam is constructed. In addition, considering that most of the damage models adopted by the existing mesoscale simulation of concrete damage and fracture cannot take into account the interaction between aggregates, interfacial transition zone (ITZ), and mortar, an improved anisotropic temperature damage model is proposed in this paper. The aggregate placement simulation results show that the method presented in this paper can quickly generate two-dimensional (2D) random concrete aggregates, and the generation of three-dimensional (3D) aggregates can also be completed in a very short time, which can greatly improve the aggregate generation efficiency. Moreover, the aggregate shape generated by this method is very close to the real concrete aggregate shape. The crack propagation simulation results show that the sudden rise and fall of temperature can cause damage in the mortar and ITZ of concrete inside the dam body, which is the main reason for the generation of macroscopic through-cracks in the No. 9 non-overflow dam section of Shimantan dam during the operation period. Finally, it can be learned from the results that the method presented in this paper is reasonable and feasible, and can be extended to the crack propagation simulation of some other concrete gravity and arch dams.