The core of mechanical behaviors simulation for fiber reinforced plastic (FRP) laminated structures is the construction of both material components′ mechanical responses and the coupling mechanisms between them. Simulation ultimate aim is to obtain the evolution process of timerelated damage spatial field. Damage spatial field, materialproperty field and stressstrain distribution field interact with each other. The whole logical framework which describes the coupling relationship among three spatial fields can be defined as numerical analysis architecture. Within this architecture, material components′ mechanical responses, nonlinear numerical analysis, damage initiation criterion, damage evolution criterion, size effect and insitu effect are main modules. In this paper, several methods that can constitute these modules are discussed, and the promising ones are suggested to obtain a better simulation architecture. Finally, numerical examples under three different loading cases are conducted and compared with corresponding experimental data. The comparison results show the reasonability and effectiveness of the numerical analysis architecture.