Research Article  Open Access
Simplified Material Solution for Orthotropic Symmetrical GFRP Laminates for Structural Facades
Abstract
GFRP (glassfiberreinforced polymer), as a composite material, possesses many favorable properties including high strength and low weight and is amenable to unique processing methods; therefore, it is a potential freeform surface material. However, the complex design theory owing to anisotropy limits its application. Thus, a simplified material solution becomes significant. In this study, the strength and stiffness of orthotropic symmetrical GFRP laminates are derived theoretically, and a simplified material solution is proposed to simplify the anisotropy as isotropy. Then, using the numerical simulation of an actual orthotropic symmetrical GFRP laminate freeform facade structure, the effectiveness of the simplified material solution is analyzed and evaluated. This solution can provide guidance for similar GFRP facades and further promote the application of GFRP in engineering.
1. Introduction
Recent times have seen the construction of a number of specialshaped buildings around the world. However, the realization of these complex architectural shapes poses great challenges for engineers. Traditional facade materials like glass, metal, and concrete cannot be used to construct freeform facades without significant construction time and costs considering the current technical conditions. GFRP (glassfiberreinforced polymer), as composite material, possesses many favorable properties including high strength and low weight and is amenable to unique processing methods. It seems to be a potential solution for freeform facades.
GFRP represents a series of composite materials that employs a polymer as the matrix and glass fiber as the reinforcing material. Although it is a relatively new material compared to traditional materials such as steel and concrete, it has been applied in civil engineering ever since its invention. However, its applications were limited until recent years. Of late, GFRP has been used to fabricate bridge decks, structural reinforcements, bars, and decoration materials.
As reinforcement material, Kurt [1] tested the properties of plastic pipes filled with concrete and used to strengthen concrete columns in 1977. Although the plastic pipes are not strictly fiberreinforced polymers (FRP, GFRP is a typical FRP), this could be regarded as a study that pioneered composite material reinforcement. Since then, researchers from Japan and the USA worked continuously on FRP applications, and multiple achievements were made in terms of reinforcement [2]. In 1988, carbonfiberreinforced polymer (CFRP) stranded wires were applied as prestressed reinforcements on the new Shinmiya bridge in Japan [3, 4]. In 2003, B. Täljsten studied the strengthening concrete beams for shear with CFRP sheets and briefly concluded how to design it [5]. Zadeh and Nanni presented methodologies for the design of concrete columns in their articles and highlight some important aspects affecting computations [6, 7]. Based on a large database collected from the literature, Jumaa and Yousif proposed three prediction models to estimate the shear capacity of FRPRC members using artificial neural networks, gene expression programming, and nonlinear regression analysis [8]. Currently, composite bridge construction is the major application of FRP structures and also the most technologymature fullFRP structure form now available. Composite bridges were first built in the 1980s, where FRP was employed as deck and beam materials. In 1983, an FRP honeycomb box girder highway bridge was built in Miyun, Beijing [9], which was the first FRP highway bridge. FRP pedestrian bridges are highly praised due to the lightweight and the ease in fabrication and installation. Many innovative FRP pedestrian bridges have been constructed in Asia, Europe, and North America [10, 11]. FRP is also applied as link slabs of bridge decks. In 2012, Saber and Ashok investigated the use of glassfiberreinforced polymers grid for reinforcement in link slabs for jointless bridge decks [12].
FRP was first employed as a building facade material in the 1950s, when composite materials were just beginning to find civilian applications. Although several creative FRP house prototypes were developed, only a few of them were applied successfully. In 1956, the Maison en Plastique, which was regarded as the first ever structure with loadbearing FRP members, was set up in Paris [13]. The structure was designed for the Salon des Arts; inspired by the snail, the house could be enlarged using addition modules around the circular core. In the same year, the famous Monsanto House of the Future [14], whose Cshaped surface was made of GFRP, was constructed in Disneyland, Orlando. Futuro, 1968 [15], possessed a distinctive GFRPsaucer shape and was the prototype of many subsequent GFRP structures. In 2005, the Yitzhak Rabin Center was built in Tel Aviv, Israël. To commemorate Yitzhak Rabin’s contribution to peace, the roof of the great hall and library are designed to resemble dove wings using composite materials. Similarly, the Mobile Art Chanel Contemporary Art Container was designed by Zaha Hadid Architecture and started its travel in 2008 [16]. The facade of the freeform “mobile” structure is made of FRP; the ease of processing and low weight characteristics are extremely valuable in the fabrication of this structure. Low cost, no heavy machinery requirements, and convenience of assembly, disassembly, and transport make FRP freeform facade structures desirable.
The inherent anisotropy of FRP is not very difficult for researchers. However, it is too complicated for applications in actual engineering. It results in complicated design theories that, in turn, require extensive calculations, thus directly limiting its applications. In recent years, some simplified calculation methods were proposed for composite materials. In 2014, a simplified reliabilitybased determination procedure on the strengthening ratio of a deteriorated concrete girder with CFRP strips was proposed [17]. In 2018, Ozbakkaloglu et al. proposed a simple yet powerful designoriented model to predict the axial stressstrain behavior of FRPconfined concrete in circular sections [18–20]. In 2019, Zhang et al. proposed a simplified theoretical solution for the evaluation of the equivalent bending stiffness of a unique hybrid FRPaluminum spatial decktruss structure [21]. However, simplified design methods for pure FRP members in engineering are still rare.
As facade material or loadbearing single shell members, the utilization value of orthotropic symmetrical GFRP laminates is the highest in terms of cost, process, and performance. The orthotropic symmetrical GFRP laminates in this paper are reinforced equally in two mutually perpendicular directions, such as and . This study proposes a simplified material solution to transform the anisotropy of orthotropic GFRP laminates into isotropy.
This study involves the theoretical derivation of the strength and stiffness of orthotropic symmetrical GFRP laminates and obtains the offaxis strength and stiffness. Bending experiments and corresponding finite element method analysis are carried out, and their results are compared to verify the consistency. Numerical simulation of an actual GFRP freeform facade is performed, in which the GFRP material and simplified material properties are employed, respectively, to examine the accuracy and effectiveness of the simplified solution. The simplified material solution is ultimately applied to the design of the actual GFRP freeform facade, and the construction of the freeform facade was completed during early December 2017.
2. Simplification Necessity
Theoretical research on the anisotropic properties of composite materials has always been the focus for composite material research. Although some theories are in good agreement with the actual performance of composites materials under specific conditions, their value in practical application is limited and experiments are always required to further examine FRP products. Generally, engineers are cautious when assuming composite material property values. Most current FRPrelated codes and standards in civil engineering mainly deal with structural reinforcement, and some of them involve FRP pultruded profiles [22–26]. These codes provide some regulations and criteria for FRP indexes, but there are no recommended or convenient calculation methods to compute these indexes.
Unidirectional reinforced GFRP can only bear load in the longitudinal direction as its properties in the transverse direction are too weak, and when applied to unidirectional members, it is almost identical to traditional isotropic materials. For multidirectional loadbearing members, orthotropic symmetrical GFRP laminates are mostly applied due to their relatively low cost, ease of processing, and relatively stable properties. The anisotropy of orthotropic symmetrical GFRP laminate is significant, which distinguishes it from conventional building materials. Generally, the stiffness of FRP is calculated based on the classical laminate theory, and failure criteria include the maximum stress criterion, maximum strain criterion, Tsai–Hill criterion, and Tsai–Wu criterion. The Tsai–Hill and Tsai–Wu criteria [27, 28] are widely applied for anisotropic materials, and codes such as Eurocode [24] use the Tsai–Hill criterion to verify ply level failure, and the maximum strain criterion for preliminary design to verify laminatelevel failure. However, the calculation of these stresses and strains is uncertain in the codes, and it is suggested that experiments should be conducted. Most designers and engineers are unfamiliar with anisotropic materials, not to mention the calculations regarding composite materials like GFRP. This is a direct factor that limits the application of GFRP in civil engineering.
In addition, even if designers are familiar with GFRP, it is still difficult to analyze orthotropic symmetrical GFRP laminate structures directly using composite material properties, especially when the calculation involves a large number of panels of various shapes as in a GFRP freeform facade. As a multilayer material, a single GFRP panel can be considered as a composite member, so that the calculation requirement is equivalent to that of a small structure. Therefore, the calculations pertaining to an entire structure is equivalent to several structures at least. For buildings or structures, macroscopic indexes like the maximum deformation are the most relevant factors, and the working environments of buildings are complicated so that a certain error in the calculation is always acceptable. Therefore, the practical significance of laminatelevel calculation is relatively limited, and it is necessary to propose a simplified material solution.
The simplified material solution should satisfy three requirements:(a)The calculated results pertaining to members or structures that apply simplified material properties are in good agreement with the original properties(b)The simplified material properties can be obtained conveniently through specific experiments(c)The simplified material solution should guarantee the safety of the structure, and so it should be a relatively conservative solution
Therefore, the target is to simplify orthotropic symmetrical laminates into the isotropic material, thus avoiding anisotropic properties, specifying appropriate stiffness and strength parameters for design reference, reducing design difficulty, and further promoting the application of GFRP in civil engineering.
3. Theoretical Derivation
3.1. Stiffness
According to the classical laminate theory, the loadstrain matrix of the laminate is shown in the following equation:where N is the inplane force vector, M is the inplane moment (torque) vector, A is the inplane stiffness matrix, B is the coupling stiffness matrix, D is the bending stiffness matrix, is the midplane strain vector, and is the midplane curvature vector.
An orthotropic symmetrical laminate is a special laminate having coupling stiffness matrix B equal to zero, and the 16 and 26 matrix components of A and D are also zero. Thus, no tensileshear coupling and tensilebending coupling exist. Therefore, calculations of orthotropic symmetrical laminates are actually the same as that for each ply, especially for laminates reinforced by bidirectional woven glass fabrics.
The thickness of a panel is far less than its length and width; therefore, the analysis of the panels could be assumed as plane stress problems: . The axes are shown in Figure 1, and the stressstrain relationships are shown in the matrix form as follows.
The axes directions and offaxis directions are as follows:where the flexibility matrix S can be expressed aswhere is the elastic modulus of 1axis, is the elastic modulus of 2axis, is the Poisson ratio of 1axis, is the Poisson ratio of 2axis, and is the shear modulus in 12.
The offaxis stiffness is lower than that in the reinforcement direction, and the load may be transmitted along the offaxis direction for some irregular shell members. Therefore, the offaxis properties actually determine the lower limit of material properties.
In Figure 1, x and y are offaxes, and the offaxis stressstrain relationship can be obtained through the transformation formula. Thus, equation (2) can be transformed into the following equation:where
Therefore, the offaxis stiffness properties can be expressed as follows:where is the elastic modulus of xaxis, is the elastic modulus of the yaxis, is the Poisson ratio of x, is the Poisson ratio of y, and is the shear modulus in xy.
For orthotropic symmetrical laminates, . The derivation of is shown in the following equation:
Therefore, when
is at the extremum. Equation (9) corresponds to the isotropic material. For orthotropic symmetrical laminates, is minimum and is given by . The simplification focuses on the elastic modulus as it is the main factor of panel deformation, and may be taken as the elastic modulus of a simplified isotropic material.
The stiffness properties of orthotropic symmetrical GFRP laminates applied in the actual structure described in chapter 4 are shown in Table 1.

Substituting the properties in Table 1 into equation (6), and the offaxis property variation curves with are shown in Figures 2–4.
When , , which is minimum; , which are maximum. In contrast with isotropic materials, Poisson’s ratio , which is greater than 0.5, is possible for anisotropic material, and it could be less than 0 as well.
3.2. Strength
Here, the Tsai–Hill criterion is applied as the failure criterion of GFRP laminate. The general form of the formula is given in following equation:
In the case of plane stress problems, . For orthotropic symmetrical GFRP laminates, the fibers are distributed equally in the 1 and 2 directions and , where X and Y are the axis strengths in longitudinal and transverse directions, respectively. Therefore, the strength parameters H, G, F, and N are confirmed and shown in the following equations:where Z is the strength in thickness direction and S is the shear strength in 12.
So, equation (10) can be expressed as equation (15) for orthotropic symmetrical laminates:
When the laminates are subjected to unidirectional load, , and can be transformed into (stress in offaxis) as follows:
Then, equation (15) can be converted into following equation:
The derivation of is shown in the following equation:
Therefore, is minimum: , and this strength can be a reference parameter for simplified material strength.
Strength properties of orthotropic symmetrical GFRP laminates applied in the actual structure described in chapter 4 are shown in Table 2.
 
where is the tension strength of 1axis (the same as 2axis), is the compression strength of 1axis, is the tension strength in thickness direction, and is the compression strength in thickness direction. 
Substitute the properties of Table 2 into equation (17), and the offaxis strength variation curves with are shown in Figures 5 and 6.
Appropriate failure criteria and strength parameters must be selected for the simplified isotropic material. The maximum stress criterion assumes that a material fails when the maximum principal stress exceeds the uniaxial tensile strength, and it is commonly applied among brittle materials like concrete. GFRP is essentially a brittle material and therefore, the maximum stress criterion might be appropriate for the simplified isotropic material. The stressstrain curves of brittle materials in uniaxial load experiments are always oblique lines terminated at the ultimate strain values, as shown in Figure 7(a), and the design strengths are commonly obtained by dividing the ultimate stresses by safety factors. Thus, many studies focus on the determination of this safety factor. However, another method to obtain the design strength of simplified isotropic material is proposed here.
(a)
(b)
Members like facade panels are subjected primarily to bending loads. Bending strength and bending modulus, which are obtained from bending tests, are common parameters for FRP. There are many factors that influence bending properties, including the type of polymer and fiber and the plystacking sequence. And thus, bending properties are a reflection of the overall material properties. A typical stressstrain curve of FRP in bending test is shown in Figure 7(b), which is totally different from the curve of brittle materials in Figure 7(a) and is indeed possible especially for highly orthotropic laminates. There is a long horizontal section, and it seems to be a curve of a ductile material.
The formulas provided in bending test standards, such as ASTM D790 [29], to calculate the bending stress and strain in Figure 7(b) are based on homogeneous elastic material. Actually, they are not totally appropriate for highly orthotropic laminates. The specific explanation is very complicated and different for different FRP. However, bending properties are applied as parameters of the simplified solution in this paper, and therefore, only a brief explanation is listed below, which could be applied for common GFRP, including the GFRP applied in this paper.
Within the elastic range of the curve in Figure 7(b), bending stresses are equivalent to the compressive and tensile stresses on the surface of the specimen, and the same applies to bending modulus. The elastic bending properties can be computed theoretically through formulas, and the practical elastic bending ultimate strength can be obtained at the end of the elastic range.
The compression strength of polymers (the matrix of GFRP) is generally higher than its tensile strength, and so, in bending experiments, the polymers in tensile side cracks and fails first (“crack point” in Figure 7(b)). Though the glass fiber of the tensile side can still bear the load, the bending stress formula is invalid. Therefore, the curve after the crack point in Figure 7(b) is a nominal stressstrain curve and the bending ultimate strength is a nominal strength, which is often applied as a qualitative criterion.
The elastic bending ultimate strength (stress on the crack point of Figure 7(b)) of GFRP at is generally much less than its ultimate strength. But when applied as facade panels, stresses in the GFRP panel under a common load are much lower than the GFRP ultimate strengths and even lower than the elastic bending ultimate strength. The design of a GFRP panel is mainly controlled by deformation in general. Therefore, the elastic bending ultimate strength of GFRP at may be applied directly as the design strength of the maximum stress criterion for simplified isotropic material.
3.3. Simplified Material Solution
The aim of this study is to simplify and idealize orthotropic symmetrical GFRP laminates into isotropic material. Theoretically, the stiffness and strength of orthotropic symmetrical GFRP laminates are both minimum at ; therefore, it is safe to consider properties as the basic properties of the simplified isotropic material; GFRP is a brittle material, and thus, the maximum stress criterion is employed as the failure criterion. GFRP facade panels are mainly subjected to bending loads. Hence, through bending tests at of orthotropic symmetrical GFRP laminates, the bending modulus can be employed as the elastic modulus of simplified isotropic material and elastic bending ultimate strength can be employed as the failure strength. An isotropic material has only two independent parameters, which often adopt the elastic modulus and Poisson’s ratio in engineering, so the shear modulus is not independent. Poisson’s ratio has a relatively small effect on deformation and can be assumed as 0.3 for common orthotropic symmetrical GFRP laminates.
A brief solution is shown in Table 3.

4. Overview of the Structure
In 2017, a structure with a freeform facade was built in Fuzhou, China. The architectural schematic is a sculpturelike structure representing three morning glories: a considerably freeform facade with strong visual impact. The facade model is shown in Figure 8. To realize the architecture, orthotropic symmetrical GFRP laminates were finally applied as facade materials, which could reduce construction time as well as cost. The simplified material solution was applied to calculate facade panel parameters in the design phase, and the structure was completed in December 2017, as shown in Figure 9. The simplified material solution will be evaluated by bending tests and numerical simulation of the actual structure.
The height of the structure is 9.8 m, and the surface area is approximately 1460 m^{2}. The structure is completely exposed to open air, covering the entire sunken plaza and protecting it from rain. The inner supporting structures are steel tubes and Tsection purlins. As Fuzhou is located near the coast and affected by Pacific typhoons in summer, the basic wind load is relatively large, and the wind load must include both wind suction and pressure as an open structure. According to the Load Code for the Design of Building Structures [30], the basic wind pressure is 0.7 kN/m^{2} in Fuzhou. In the design phase of the structure, a load combination of 2.8 kN/m^{2} was applied to evaluate the deformation of members, and a load combination of 4 kN/m^{2} was applied to evaluate the strength. In the numerical simulation of chapter 6 in this paper, a load of 4 kN/m^{2} is applied on the panel to calculate the deformation and stresses.
5. Bending Test
5.1. Test Scope and Specimen
The bending properties of orthotropic symmetrical GFRP laminates applied in this structure are tested to obtain the bending modulus and elastic bending ultimate strength at . The test procedure refers to ASTM D790: Standard Test Methods for Flexural Properties of Unreinforced and Reinforced Plastics and Electrical Insulating Materials [29]. Two batches of specimens were tested, and their material ingredients are shown in Table 4. The additive, ammonium polyphosphate, in Table 4 is a flame retardant, which is mixed in material B. A limited amount of additive has no influence on the material properties since FRP properties are mainly determined by the fiber. The axis properties of the two batches are shown in Table 5.
 
EGlass: alkalifree glass fiber. 

According to ASTM D790 [29], the specimen dimensions were prepared following Figure 10 and Table 6, and five specimens were prepared for each batch.
(a)
(b)

5.2. Test Result
The loading process and failed specimens are shown in Figures 11 and 12.
(a)
(b)
The tensile and compressive stresses and strain on specimen surfaces can be calculated using equations (19) and (20) within the elastic range:where is the load in the middle, D is the deformation in the middle, is the flexural stress, and is the flexural strain.
After the tensile side cracks, equations (19) and (20) become invalid and the results calculated from them become nominal stresses and strains. The curves of the tests are shown in Figure 13.
(a)
(b)
The bending elastic modulus can be obtained from the test data when and by equation (21). The bending elastic moduli of specimens are shown in Table 7.

The average elastic moduli are 9.86 GPa and 9.71 GPa, respectively, which is close to 9.074 GPa calculated theoretically but a little higher. The error is probably induced by two factors:(a)9.074 GPa is calculated based on the properties in Table 1. The properties in Tables 1 and 2 are provided by the GFRP manufacturer, and these properties are applied for design and calculation. Therefore, these properties might be conservative values and the real average values of these properties might be a little larger.(b)The test error of is given in Table 1. The test methods for shear modulus of FRP are much more complicated compared with methods for tensile and compressive properties, and the principles of the test methods induce errors themselves. However, the value of indeed influences the value of of orthotropic symmetrical GFRP laminates. The values of , , and are generally reliable. Therefore, assuming the values of , , and in Table 1 are accurate, and according to equation (6) and the bending elastic moduli which were just obtained, 3.3 GPa might be a more proper value for .
Yet, the little error does not influence the application of the simplified material solution.
In Figure 13, when the stresses are approximate 60 MPa, the tangent moduli of curves start to decline. Thus, 60 MPa could be adopted as the ultimate elastic bending strength. At last, the stresses stop rising and the curves almost become horizontal. The nominal bending ultimate strengths in Figure 13 are approximate 120 MPa and 140 MPa, respectively.
5.3. Numerical Simulation of the Bending Test
Numerical simulation of the bending test is conducted in Abaqus™. The multilayer and laminate material simulation functions of Abaqus™ are relatively mature and widely used in the simulation of composite materials. The laminate material parameters and model dimensions used in Abaqus™ are listed in Tables 5 and 6, respectively; The reinforcement direction is , and therefore, it is equivalent to test the bending properties of reinforced laminates.
The loads at the center of the model are 74.35 N and 61.66 N, respectively, and these are the forces when in the tests. The models are articulated at the same position as the bending tests. Geometric nonlinearity should be considered; however, as load and deformation are relatively small, it is not that important. The results of the calculation are shown in Figure 14.
(a)
(b)
The deformation at the center of materials A and B are 0.650 mm and 0.555 mm, respectively. The bending modulus can be calculated through the following equation:
The bending moduli of materials A and B at , obtained from the numerical simulation, are 8.753 GPa and 8.711 GPa, which are considerably close to the theoretical value of 9.074 GPa, and there is still a little error compared with 9.86 GPa and 9.71 GPa obtained from bending tests.
6. Numerical Simulation of the Actual Structure
6.1. Scope of the Simulation
The facade divisions and supporting purlins are shown in Figure 15. Every panel is a freeform panel that is different from each other. The maximum dimension of the panels is approximately . Fifteen panels, which are typical in position and shapes as shown in Figure 15, are selected for simulation in Abaqus™.
The fifteen selected panels include extremely twisted shapes as well as several flat ones. Simplified isotropic material properties and original orthotropic symmetrical GFRP laminate properties are given to the panels, respectively, following which the deformations and the stresses are calculated and the values of Tsai–Hill failure criterion of the laminate are evaluated. The specific analysis procedures are shown as follows.
Boundary conditions: panels are relatively thin and the elastic modulus is small; thus, the stiffness is small. Simultaneously, the panels are connected to purlins using structural adhesive. Therefore, the boundaries are assumed to be simply supported on all four sides.
Load: GFRP is light in weight and the panels are mainly subjected to wind loads. The wind load in Fuzhou is considerable, and the wind pressure and wind suction must be considered simultaneously. A load of 4 kN/m^{2} is applied in the numerical simulation after the combination.
Materials: this analysis focuses on the comparison of panels with different material properties, namely, the simplified isotropic material and original orthotropic symmetrical GFRP laminate. For GFRP laminate material, the panel depth of the structure is 4.8 mm, 0.6 mm for each layer and eight layers in total, and the material properties are the same as the material applied in the actual structure, as shown in Tables 1 and 2. The reinforcement directions of GFRP panels influence the deformation and stresses significantly. However, the shapes of freeform panels in Figure 15 are not regular. It is not realistic to control the reinforcement directions of every single panel during processing, especially for irregular shapes. Thus, the reinforcement directions of panels in the simulation are at the default setting of in the local coordinates directly. For simplified isotropic material, the depth of the panel is also 4.8 mm, while the elastic modulus varies as 9 GPa, 10 GPa, 11 GPa, 12 GPa, 13 GPa, 14 GPa, and 15 GPa, respectively, to evaluate errors, and Poisson’s ratio is assumed to be 0.3.
Analysis steps: the finite element method analysis must consider the geometric nonlinearity as the stiffness is typically low; the calculation includes the deformation and stresses of every mesh for all the isotropic materials and GFRP laminate material.
6.2. Numerical Simulation Result
The results of the deformations are listed in Table 8 and Figure 16. In Table 8, numbers in the first column are the panel number; values in the second column are the maximum deformations of every panel when adopting GFRP laminate property; values in the rest columns are the maximum deformations of every panel when adopting simplified isotropic materials with different elastic modulus.

The lines in Figure 16 are drawn by connecting the deformations in Table 8 to provide a more obvious contrast of the data. The bold line belongs to GFRP, and the others thin lines belong to simplified isotropic material with different elastic moduli.
The results of stresses in the panels are listed in Table 9. The numbers in the first column are the panel numbers. For simplified isotropic properties, since the influence of elastic modulus on the stresses is limited, only the largest absolute values of principal stresses in the panels when employing an elastic modulus of 13 GPa are listed in the second column; and for the GFRP laminate property, the largest absolute values of positiveaxis stresses and the maximum values of Tsai–Hill failure criterion are listed.
 
where S11 is the stress in 1axis, S22 is the stress in 2axis, and S12 is the shear stress in 12. 
6.3. Numerical Simulation Analysis
6.3.1. Strength
In Table 9, the largest principal stress of the simplified isotropic panels is 26.53 MPa and the largest positiveaxis stress of the GFRP laminate panels is 21.78 MPa. The stresses are far smaller than the ultimate strength and even smaller than the elastic bending ultimate strength of the GFRP, which is 60 MPa. The value of Tsa–Hill failure criterion should be between 0 and 1, where 1 means failure. The largest value of Tsai–Hill failure criterion in Table 9 is only 0.1764. In general, strengths are not the key factors for GFRP facade panels, since the GFRP strengths are large enough.
6.3.2. Stiffness
In Table 8 and Figure 16, the deformation tendency of the simplified isotropic material panels is consistent with laminate panels. Assume the GFRP laminate deformations are the accurate values, and therefore, subtraction of GFRP laminate deformations and isotropic materials deformations can be considered as errors induced by material simplification. The errors norms are shown in Table 10.

The norms of 12 GPa are the least in Table 10; therefore, simplified isotropic material with an elastic modulus of 12 GPa is the material most similar to the GFRP laminate. However, the elastic modulus should be smaller than 12 GPa for safety considerations. In Figure 16, the 9 GPa curve is higher than the GFRP curve, and the 10 GPa curve is almost higher than the GFRP curve but only a little lower at panel11. The theoretical elastic modulus of GFRP at is approximately 9 GPa, while the experimental value is approximately 10 GPa. Thus, it is reasonable to adopt the bending modulus of the orthotropic symmetrical GFRP laminate at as elastic modulus of simplified isotropic material.
In Table 8, the largest deformation is near 30 mm. The Technical code for building curtain wall [31] stipulates that the maximum deformation should be no more than 1/60 of the span. “1/60 of the span” is for glass, and it is the largest regulation. The regulation is different for different materials, and there is no regulation provided in the code for composite materials like GFRP. The spans of the panels are less than 2000 mm, so that 30 mm is almost 1/60. But the load applied in the numerical simulation is to evaluate the strength during the design phase of the structure. The load applied to evaluate deformation is smaller. And actually, the deformation ability of GFRP is much better than that of glass. Therefore, the stiffness of GFRP is enough considering function and safety problems.
6.4. PanelSpecific Analysis
In Table 8 and Figure 16, the deformation of panel8 is particularly presented. The deformation of GFRP laminate property is much less than the deformation of isotropic materials. The shape and reinforced direction of panel8 are shown in Figure 17.
The reinforcement direction could not be completely perpendicular to the edges owing to its freeform shape. Rotating the reinforcement direction of panel8 and recalculating its deformations, the largest deformation rises to 21.2 mm after a anticlockwise rotation of the reinforcement direction. This phenomenon indicates the customizability of GFRP properties. Under certain conditions, it is a mathematical extremum problem to obtain the most appropriate material properties and the best performance of the composite material. Considering panel8, its minimum deformation can be obtained through rotation of orthotropic reinforcement direction or by abandoning orthotropic properties and rearranging the reinforcement directions.
7. Discussion and Conclusion
This paper proposed a simplified material solution for orthotropic symmetrical GFRP laminates used in structural facades. Reducing the difficulty of design as well as computational complexity can help promote the application of GFRP in civil engineering. Theoretically, the elastic modulus and strength in of the orthotropic symmetrical GFRP laminates are the smallest and can be applied as simplified material indexes. The elastic modulus and strength of the simplified material could take on the value of the bending elastic modulus and elastic bending ultimate strength that is obtained directly through the bending test of GFRP laminates in . It is also proposed that the maximum stress criterion can be applied as the failure criterion of the simplified material model, and Poisson's ratio can take the value of 0.3.
This simplified material solution is convenient to use and recommended for similar structures as reference. The following conclusions can be obtained through this study:(1)The computational complexity and quantity of laminates are relatively large. The stiffness matrix of anisotropic laminates is highly complex, with an 8layer panel necessitating at least eight times the computational requirement. Coordination conditions between the layers further complicate the calculation, especially for freeform panels considering the recalculation and verification induced by design adjustments. In general, the analysis of integral structures by considering laminate material properties is very difficult, and the accuracy of the result might be doubtful. Even for single panels in this paper, multiple adjustments to element meshing and load steps have to be made to obtain reliable calculation results for some complicated panels when applying GFRP laminate properties. It is almost impossible to calculate GFRP facade structures by applying its original properties; thus, the simplified material solution obtained here becomes significant in promoting GFRP applications.(2)GFRP panels are stiffnesscontrolled, and the strength is large enough. The applied value of 4 kN/m^{2} is a relatively large load combination. However, the largest value of Tsai–Hill failure criterion is only 0.1764, while the largest principal tensile stress of simplified isotropic material is only 26.53 MPa, which is far smaller than the elastic bending ultimate strength, 60 MPa. The deformation is relatively large owing to the small stiffness of the GFRP. Therefore, the main design factor for the GFRP panels is deformation. In the Chinese code Technical code for building curtain wall, the regulation of glass facade, 1/60 of the span, is the largest deformation regulation, while the deformation ability of GFRP is much better than that of glass. Hence, it is suggested that the deformation regulation of GFRP could be larger than 1/60.(3)The simplified material solution is effective. Deformation of the simplified isotropic material is consistent with that of GFRP laminate material properties and slightly higher considering the safety. Although the elastic bending ultimate strength of 60 MPa seems to be much smaller than the ultimate strength of the GFRP, it is still sufficient for facade panels as design strength, even for extreme environmental conditions.
Data Availability
The experimental data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (Grant no. 51578404).
References
 C. E. Kurt, “Concrete filled structural plastic columns,” Journal of the Structural Division, vol. 104, no. 1, pp. 55–63, 1978. View at: Google Scholar
 C. E. Bakis, L. C. Bank, V. L. Brown et al., “Fiberreinforced polymer composites for constructionstateoftheart review,” Journal of Composites for Construction, vol. 6, no. 2, pp. 73–87, 2002. View at: Publisher Site  Google Scholar
 H. T. Nguyen, H. Masuya, T. M. Ha et al., “Longterm application of carbon fiber composite cable tendon in the prestressed concrete bridgeshinmiya bridge in Japan,” MATEC Web of Conferences, vol. 206, 2018. View at: Publisher Site  Google Scholar
 T. Yamashita, T. Kiuchi, H. Inukai et al., “PC bridge with new materialconstruction of Shinmiya Bridge,” Prestressed Concrete, vol. 31, no. 2, pp. 71–78, 1989. View at: Google Scholar
 B. Täljsten, “Strengthening concrete beams for shear with CFRP sheets,” Construction and Building Materials, vol. 17, no. 1, pp. 15–26, 2003. View at: Publisher Site  Google Scholar
 H. J. Zadeh and A. Nanni, “Design of RC columns using glass FRP reinforcement,” Journal of Composites for Construction, vol. 17, no. 3, pp. 294–304, 2012. View at: Publisher Site  Google Scholar
 H. J. Zadeh and A. Nanni, “Reliability analysis of concrete beams internally reinforced with fiberreinforced polymer bars,” Structural Journal, vol. 110, no. 6, pp. 1023–1032, 2013. View at: Publisher Site  Google Scholar
 G. B. Jumaa and A. R. Yousif, “Predicting shear capacity of FRPreinforced concrete beams without stirrups by artificial neural networks, gene expression programming, and regression analysis,” Advances in Civil Engineering, vol. 2018, Article ID 5157824, 16 pages, 2018. View at: Publisher Site  Google Scholar
 G. D. Tang, Y. Tang, and G. Z. Feng, “GRP bridges in China: research and practice,” Journal of Chengdu University of Science and Technology, vol. 6, pp. 69–80, 1995. View at: Google Scholar
 P. J. D. Mendes, J. A. O. Barros, J. M. SenaCruz, and M. Taheri, “Development of a pedestrian bridge with GFRP profiles and fiber reinforced selfcompacting concrete deck,” Composite Structures, vol. 93, no. 11, pp. 2969–2982, 2011. View at: Publisher Site  Google Scholar
 L. C. Hollaway, “A review of the present and future utilisation of FRP composites in the civil infrastructure with reference to their important inservice properties,” Construction and Building Materials, vol. 24, no. 12, pp. 2419–2445, 2010. View at: Publisher Site  Google Scholar
 A. Saber and A. R. Ashok, “Behavior of FRP link slabs in jointless bridge decks,” Advances in Civil Engineering, vol. 2012, Article ID 452987, 9 pages, 2012. View at: Publisher Site  Google Scholar
 A. Quarmby, The Plastics Architect, Pall Mall Press, London, UK, 1974.
 J. L. Meikle, American Plastic: A Cultural History, Rutgers University Press, New Brunswick, NJ, USA, 1995.
 M. Taanila, FuturoTomorrow’s House from Yesterday, Desura, Helsinki, Finland, 2004.
 D. Kendall, “The business case for composites in construction,” Reinforced Plastics, vol. 52, no. 7, pp. 20–27, 2008. View at: Publisher Site  Google Scholar
 M. Ju, H. Oh, and J. W. Sun, “Simplified reliability estimation for optimum strengthening ratio of 30yearold double Tbeam railway bridge by NSM techniques,” Mathematical Problems in Engineering, vol. 2014, Article ID 734016, 10 pages, 2014. View at: Publisher Site  Google Scholar
 A. F. Pour, T. Ozbakkaloglu, and T. Vincent, “Simplified designoriented axial stressstrain model for FRPconfined normal and highstrength concrete,” Engineering Structures, vol. 175, pp. 501–516, 2018. View at: Publisher Site  Google Scholar
 T. Ozbakkaloglu and J. C. Lim, “Axial compressive behavior of FRPconfined concrete: experimental test database and a new designoriented model,” Composites Part B: Engineering, vol. 55, pp. 607–634, 2013. View at: Publisher Site  Google Scholar
 A. F. Pour, A. Gholampour, J. Zheng, and T. Ozbakkaloglu, “Behavior of FRPconfined highstrength concrete under eccentric compression: tests on concretefilled FRP tube columns,” Composite Structures, vol. 220, pp. 261–272, 2019. View at: Publisher Site  Google Scholar
 D. Zhang, F. Li, F. Shao, and C. Fan, “Evaluation of equivalent bending stiffness by simplified theoretical solution for an FRPaluminum decktruss structure,” KSCE Journal of Civil Engineering, vol. 23, no. 1, pp. 367–375, 2019. View at: Publisher Site  Google Scholar
 A. Machida and T. Uomoto, Recommendation for Design and Construction of Concrete Structures Using Continuous Fiber Reinforcing Materials, Japan Society of Civil Engineers, Tokyo, Japan, 1997.
 A. H. Zureick, B. R. Ellingwood, A. S. Nowak et al., Recommended Guide Specification for the Design of Externally Bonded FRP Systems for Repair and Strengthening of Concrete Bridge Elements, Transportation Research Board, Washington, D.C., USA, 2010.
 L. Ascione, J. F. Caron, P. Godonou et al., Prospect for New Guidance in the Design of FRP Support to the Implementation, Harmonization and Further Development of the Eurocodes, European Commission, Brussels, Belgium, 2016.
 China Metallurgical Construction Association, Technical Code for Infrastructure Application of FRP Composites, China Planning Press, Beijing, China, 2011.
 Sichuan Department of Housing and UrbanRural Development, Code for Design of Strengthening Concrete Structure, China Architecture & Building Press, Beijing, China, 2014.
 R. Hill, The Mathematical Theory of Plasticity, Oxford University Press, England, UK, 1998.
 S. W. Tsai and E. M. Wu, “A general theory of strength for anisotropic materials,” Journal of Composite Materials, vol. 5, no. 1, pp. 58–80, 1971. View at: Publisher Site  Google Scholar
 ASTM, D79010: Standard Test Methods for Flexural Properties of Unreinforced and Reinforced Plastics and Electrical Insulating Materials, American Society for Testing Materials, West Conshohocken, PA, USA, 2010.
 Ministry of Housing and UrbanRural Development of the People’s Republic of China, Load Code for the Design of Building Structures, China Architecture & Building Press, Beijing, China, 2012.
 Shanghai Metal Construction Association, DGJ08562012: Technical Code for Building Curtain Wall, Shanghai Metal Construction Association, Shanghai, China, 2012.
Copyright
Copyright © 2019 Yuchao Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.