Modeling of thick composite structures for finite element analysis is relatively complicated. 2-D plane elements may cause inaccurate result since the plane stress condition cannot be applicable in these structures. Therefore a 3-D modeling should be used. However, the difficulty to model all the layers with different material properties and ply orientation arise in this case. In this paper, an equivalent modeling is proposed and numerically tested for analysis of thick composite structures. By grouping layers with same material and ply orientation, number of elements through the thickness is remarkably reduced and still the result is close enough to the one from a detail finite element model. MSC/NASTRAN and PATRAN are used for the analysis. The proposed modeling technique has been applied for analysis of composite rotor hub system designed by Korea Aerospace Research Institute(KARI). Using the proposed equivalent modeling technique, we could conduct stress analysis for the hub system and check the safety factor of each part.