Structural analysis for materials containing regularly spaced in-homogeneities is usually executed by using averaged material properties. For the homogenization process, a unit cell is defined and loaded somehow, and its response is investigated to evaluate the properties. The imposed loading conditions should accord to the behavior of unit cell immersed in the macroscopic structure in order to guarantee the accuracy of the effective properties. Each unit cell shows periodic variation of strain if the material is loaded uniformly, and in this study, direct implementation of this characteristic behavior is attempted on FE models of unit cell. Conventional finite element analysis tool can be used without any modification, and the boundary of unit cell is constrained in a way that the periodicity is satisfied. The proposed method is applicable to skew arrayed in-homogeneity problems. The flexibility matrix relating tonsorial stress and strain components in skewed rectilinear coordinate system is transformed so that the required engineering constants can be evaluated. Effective properties are computed for the materials with square and skew arrayed circular holes, and its accuracy is examined.