“…These numerical methods that aimed to efficient achievement of modern compute essentially required finite element discretization of the inhomogeneous region, such as finite element modeling (FEM) 20 , Mechanics of structure genome 21,22 , Fast Fourier transforms 23 and so on. Different from the limited description of physical events of the materials by phenomenological model, the analytical approaches that provide a good reference for revealing the intrinsic properties of materials allow to correlate between the overall behavior of composites and the individual constituents, for instance, elasticity solution [24][25][26][27] , Eshelby theory 28 and energy method 29,30 , which is limited by the complex situation and boundary conditions. Morever, traditional micromechanics is designed to obtain the homogenized characterization at the smallest length scale of the materials 31 , unconsciously ignoring the scale correlation of material deformation mechanism 32,33 .…”