Conclusions
In this paper, a novel numerical implementation algorithm of the asymptotic homogenization theory for predicting the effective complex moduli of a composite in the frequency domain is proposed. To solve the local problem in the homogenization process, an equivalent harmonic analysis is performed, and a double-layer elements method is proposed. Through this algorithm, the viscoelastic response of the composites in the frequency domain can be obtained easily by using commercial software in a black-box role. In this way, the powerful modeling and analytical abilities of the commercial software are utilized for computing complicated unit cells. The viscoelastic properties of the particle-reinforced composites, the fiber-reinforced composites, and the octet-truss lattice materials are computed with this method to illustrate the effectiveness and simplicity of this new implementation.