• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Compressive behavior of hybrid steel-polyvinyl alcohol fiber-reinforced concrete containing fly ash and slag powder: experiments and an artificial neural network model*#

    2021-09-21 05:28:10FangyuLIUWenqiDINGYafeiQIAOLinbingWANGQiyangCHEN

    Fang-yu LIU, Wen-qi DING, Ya-fei QIAO, Lin-bing WANG, Qi-yang CHEN

    Compressive behavior of hybrid steel-polyvinyl alcohol fiber-reinforced concrete containing fly ash and slag powder: experiments and an artificial neural network model*#

    Fang-yu LIU1, Wen-qi DING2,3, Ya-fei QIAO??2,3, Lin-bing WANG??1, Qi-yang CHEN4

    2, Shanghai 200092, China Key Laboratory of Geotechnical and Underground Engineering, Ministry of Education, Shanghai 200092, China ?E-mail: yafei.qiao@#edu.cn; wangl@vt.edu

    Understanding the mechanical behavior of hybrid fiber-reinforced concrete (HFRC), a composite material, is crucial for the design of HFRC and HFRC structures. In this study, a series of compression experiments were performed on hybrid steel- polyvinyl alcohol (PVA) fiber-reinforced concrete containing fly ash and slag powder, with a focus on the fiber content/ratio effect on its compressive behavior; a new approach was built to model the compression behavior of HFRC by using an artificial neural network (ANN) method. The proposed ANN model incorporated two new developments: the prediction of the compressive stress-strain curve and consideration of 23 features of components of HFRC. To build a database for the ANN model, relevant published data were also collected. Three indices were used to train and evaluate the ANN model. To highlight the performance of the ANN model, it was compared with a traditional equation-based model. The results revealed that the relative errors of the predicted compressive strength and strain corresponding to compressive strength of the ANN model were close to 0, while the corresponding values from the equation-based model were higher. Therefore, the ANN model is better able to consider the effect of different components on the compressive behavior of HFRC in terms of compressive strength, the strain corresponding to compressive strength, and the compressive stress-strain curve. Such an ANN model could also be a good tool to predict the mechanical behavior of other composite materials.

    Experiments; Artificial neural network (ANN); Hybrid fiber-reinforced concrete (HFRC); Compressive behavior; Stress-strain curve

    1 Introduction

    Fiber-reinforced concrete (FRC), a typical composite material, has been successfully and widely used in civil infrastructure, such as tunnels (de la Fuente et al., 2012; Liao et al., 2015; Gong et al., 2017), bridges (Stoll et al., 2000), concrete structures (Sim et al., 2005), and pavements (Belletti et al., 2008). This is because FRC greatly improves tensile behavior (Xu et al., 2011), flexural behavior (Nguyen et al., 2013), and durability (Liu et al., 2019a) compared with concrete. The compression behavior and basic mechanical characteristics of concrete can also be significantly affected by the addition of fibers (Zhou et al., 2018, 2019). Such enhanced compression behavior is important for many concrete structures, especially in the joints of precast structures where stress is concentrated (Ding et al., 2020). To design a suitable FRC, different types of fibers can be used, including steel, polyvinyl alcohol (PVA), polypropylene (PP), basalt, and synthetic fibers (Sim et al., 2005; Soutsos et al., 2012; Liu et al., 2019a, 2020a). The mechanical behavior of FRC is related to the type of fiber, the geometry of fiber, and the total fiber content (Lawler et al., 2003, 2005; Liu et al., 2020b).

    According to the number of fiber types combined, FRC can be categorized as single fiber- reinforced concrete or hybrid fiber-reinforced concrete (HFRC), which includes two or more types of fibers. HFRC generally has better mechanical properties than single fiber-reinforced concrete. This is because the hybridization of different fibers can maximize their contribution to the improvement in the mechanical properties of the concrete mixture (Lawler et al., 2005). These improvements are strongly dependent on the fiber type, fiber volume content, hybrid ratio, fiber aspect ratio, and mechanical properties of the fibers (Lawler et al., 2003, 2005; Blanco et al., 2016; Liu et al., 2019b; Zhou et al., 2019). For instance, steel fibers, a type of macrofiber, could create a bridge effect to prevent the development of a major crack because of their greater length, leading to improvements in tensile strength, crack resistance, toughness, and flexural post-cracking creep behavior (Zhao et al., 2015; Serna Ros et al., 2016; Pujadas et al., 2017; Liu et al., 2020b); PVA fibers, a type of microfiber, could introduce many noticeable microcracks and delay the propagation of macrocracks, thereby improving the toughness of concrete (Lawler et al., 2005; Pujadas et al., 2014; Liu et al., 2019b). Therefore, understanding the effects of these fiber features is important for the design and modeling of HFRC.

    To describe the compression behavior of HFRC, many researchers have built models based on the constitutive models of plain concrete (Chi et al., 2014a, 2014b; Zhou et al., 2019). The effects of many features of fibers, such as their length, diameter, and volume, can be introduced into constitutive models via statistical analysis. However, as HFRC exhibits complex responses in terms of strain-softening/hardening, volumetric dilatancy, and pressure sensitivity, these constitutive models are difficult to develop from a purely mathematical perspective (Chi et al., 2014a). They also have limitations for describing the complex compressive behavior of HFRC, especially in capturing the effects of different components.

    To overcome these limitations, deep learning methods have recently been applied to predicting the mechanical properties of concrete (Mashhadban et al., 2016) and designing the mix compositions of civil materials (A?ikgen? et al., 2015). Artificial neural network (ANN), as one of the basic deep learning methods, has shown its strong application ability in civil engineering. It can be used to analyze the relationship between different types of data, produce a correct or nearly correct response to a large database, and make a good prediction for new data. The ANN method has been successfully used to model the compressive, tensile, and flexural strengths of FRC (Karahan et al., 2008; Mashhadban et al., 2016), predict the mix composition of steel fiber-reinforced concrete (A?ikgen? et al., 2015), describe the stress-strain model for fiber reinforced polymer (FRP)-confined concrete (Jiang et al., 2020), and predict the workability of concrete incorporating metakaolin and fly ash (Bai et al., 2003). Since more factors can be considered in an ANN model, it can make more accurate predictions. Moreover, an ANN model is implicit and does not need an experimental data-based formula. However, all the above studies focused on the key characteristics of the mechanical behavior of FRC. Few studies have been carried out to predict the whole stress-strain curve of HFRC using the ANN method. The authors have developed an ANN model to capture the tensile stress-strain relationship of HFRC (Liu et al., 2020a) and attempted to extend it to describe the compressive behavior in this work.

    In summary, the compression behavior of HFRC, especially the stress-strain relationship, is important for many types of concrete structures and their modeling. Conventional constitutive models have limitations in describing such complex compression behaviors. Therefore, this study aimed to reveal the compressive behavior of HFRC containing fly ash and slag powder and develop a new approach to model its compressive stress-strain curve by using the ANN method, and to highlight the capability of the ANN model by investigating the effect of different HFRC components on the compressive behavior of HFRC. To accurately represent the effect of different components of HFRC, the proposed ANN model was used to consider many important factors, including the fiber characteristics (volume, weight, length, diameter, aspect ratio, and reinforcement index), the mechanical properties of plain concrete (elastic modulus, compressive strength, and strain corresponding to compressive strength), and the composition of the concrete mixture (water binder ratio, cement, fly ash, slag powder, water, coarse aggregate, and fine aggregate).

    In the following sections, a series of compression tests on HFRC containing fly ash and slag powder were first introduced to reveal its compressive behavior. Then relevant published data were collected to extend the size and diversity of the database. After that, the ANN model was built and trained based on the database and an equation-based model was also proposed. Finally, a comparison between these two models was carried out via the description of the experimental data, followed by the sensitivity analysis of ANN model.

    2 Compressive experiments

    2.1 Materials

    The selected cement was Class 42.5R ordinary Portland cement, which had a specific surface area of 368 m2/kg. The specific surface area of fly ash was 242 m2/kg compared with 436 m2/kg for slag powder. The fly ash had a SiO2mass fraction of 54.2% compared with 36.9% for the slag powder. The maximum particle size of the fine aggregate was 5 mm, and the coarse aggregate had a continuous size of 5–20 mm. Other properties of the cement, fly ash, and slag powder can be found in (Liu et al., 2019a). The mass proportions of the concrete mixture were cement:fly ash:slag powder:fine aggregate:coarse aggregate:water:superplasticizer=225:75:75:785:1024:165:3.2, and the water-to-binder ratio was set to 0.44.These were designed according to the Chinese Codes (CAECS, 2010; MHURD, 2011a, 2011b).

    As a component of HFRC, fibers significantly affect its mechanical properties, which warrant detailed investigation. In this study, steel fiber and PVA fiber were used, and the detailed properties of these two fibers can be found in (Liu et al., 2019a). To highlight the effects of hybridization and fiber content on the compressive behavior of HFRC, various volume-based hybrid ratios of steel fiber to PVA fiber, ranging from 1:3 to 3:1, were tested, and the total fiber volume content ranged from 1% to 3%. Steel and PVA fiber combinations with fiber volume contents of 0.5%, 1.0%, and 1.5% were tested, based on published values shown to improve the mechanical properties of HFRC (Hai et al., 2016; Feng et al., 2018; Liu et al., 2019b, 2020b). Concrete reinforced with either steel or PVA fiber alone was also tested simultaneously, as well as plain concrete, collectively acting as the control group. Therefore, there are 12 groups of specimens in total and they are named as SP, wheremeans the volume content of steel fibers andmeans the volume content of PVA fibers. For example, S1.0P1.5means the specimen that consists of 1.0% steel fibers and 1.5% PVA fibers. C means the plain concrete.

    2.2 Specimen and experiment program

    All concrete specimens were mixed according to the literature, to ensure the uniform dispersion of the fibers (CAECS, 2010; MHURD, 2011a; Li and Liu, 2016; Liu et al., 2019b, 2020b). During mixing of the HFRC samples, more steel fibers had little effect on the workability of HFRC, but increasing the PVA fibers could worsen the workability due to their high-water absorption properties when the proportions of the concrete mix remained constant, as observed in other studies (Hossain et al., 2013; Tabatabaeian et al., 2017). The geometry of the specimen was 150 mm×150 mm×300 mm. A high stiffness testing machine (Fig. 1) with a capacity of 3000 kN was used to carry out the compressive tests. To obtain a better descent part of the compressive stress-strain curve, the test was controlled by a displacement control, loading at a rate of 0.1 mm/min before reaching the peak value, and at a rate of 0.06 mm/min during the post-peak stage. The loading force was recorded by a force sensor. The strain in the middle part of the specimen was measured by an extensometer (Fig. 1) with a gauge length of 100 mm, attached to the side of the specimen.

    2.3 Experimental results

    Fig. 2 presents typical measured compressive stress-strain curves. The compressive stress increased quickly and almost linearly before reaching the compressive strength at a strain of about 0.002. The compressive stress then dropped nonlinearly as the strain increased. According to Figs. 2a–2c, the peak value of the curve increased as the content of steel fibers increased, but it appeared to decrease when the volume of steel fibers was too high. In comparison, with an increase in PVA fibers, the peak value reduced slightly, and the gap between the different curves widened during the post-peak stage.

    Fig. 1 Compressive test instrument and sensors

    Fig.2 Compressive stress-strain curves: (a) S0.5P0.5, S1.0P0.5, and S1.5P0.5; (b) S0.5P1.0, S1.0P1.0, and S1.5P1.0; (c) S0.5P1.5, S1.0P1.5, and S1.5P1.5; (d) S0.5P0.5, S0.5P1.0, and S0.5P1.5; (e) S1.0P0.5, S1.0P1.0, and S1.0P1.5; (f) S1.5P0.5, S1.5P1.0, and S1.5P1.5

    Table 1 summarizes the average compressive strength and its corresponding strain for all test specimens. The compressive strength and strain corresponding to compressive strength (SCCS) depended on the fiber volume and the hybrid ratio. The single fiber-reinforced concrete (S1.0, P1.0) had a higher compressive strength than the plain concrete (C). For the same fiber volume of 1.0%, the compressive strength of the steel fiber-reinforced concrete (32.97 MPa) was similar to that of the PVA fiber-reinforced concrete (33.05 MPa). Compared with the single steel fiber-reinforced concrete, HFRC had a higher compressive strength owing to the fiber hybridization effect. However, when the total fiber volume was too high, the compressive strength appeared to decrease. Among all test specimens, the sample with a steel fiber volume of 1.0% and PVA fiber volume of 0.5% had the highest compressive strength, which was 22.3% higher than that of the plain concrete, as the strain corresponding to compressive strength decreased.

    Table 1 Compressive test results

    *tis the ratio of the compressive strength of FRC to that of plain concrete (C)

    When the steel fiber volume was constant (e.g. at 0.5%, 1.0%, or 1.5%), the compressive strength decreased with increasing PVA fiber volume (Fig. 3b). Taking the samples with a steel fiber volume of 1.0% as an example, the measured compressive strength decreased from 38.33 MPa to 34.32 MPa as the PVA fiber volume increased from 0.5% to 1.5%. In comparison, increasing the steel fiber volume increased the compressive strength when the PVA fiber volume was fixed as a constant value (Fig. 3a). However, when the steel fiber volume was too high (such as above 1.0%), the compressive strength tended to decrease when the steel fiber volume increased. For instance, when the PVA fiber volume was 0.5%, the compressive strength increased by 2.64 MPa as the steel fiber volume increased from 0.5% to 1.0%, but its value decreased by 3.66 MPa when the steel fiber volume continued to increase to 1.5%. We conclude that steel fiber plays a more important role in increasing the compressive strength of HFRC. Zhou et al. (2019) reported a similar phenomenon.

    Fig. 3 Effect of fiber type on the compressive strength of HFRC: (a) steel fiber effect; (b) PVA fiber effect

    Fig. 4 Failure images of HFRC, single fiber-reinforced concrete, and plain concrete: (a) S0.5P0.5, S0.5P1.0, and S0.5P1.5; (b) S1.0P0.5, S1.0P1.0, and S1.0P1.5; (c) S1.5P0.5, S1.5P1.0, and S1.5P1.5; (d) S1.0, P1.0, and C

    3 Relevant data from the literature

    The 36 concrete samples tested in this study were insufficient to represent the diversity of HFRC. Therefore, more HFRC data from the literature were collected (Table 2). Although these samples had different cementitious materials, this seems to have had no clear influence on the hybridization of steel and PVA fibers (Liu et al., 2020b). Therefore, all these data, along with those obtained from the above tests were used to prepare the database in this study. Table 3 summarizes the mean, median, and standard deviation of each feature for this database. The size of the whole dataset is 44 135. Compared with the compressive stress and the compressive strength, these three indices of the strain corresponding to compressive strength are much lower.

    4 Performance of the ANN model

    The ANN model that described the compressive behavior was built based on the previous work (Liu et al., 2020a) and its performance was investigated in this section. The details about the ANN model can be found in theelectronic supplementary materials. Moreover, a traditional equation-based model was also introduced for comparison with the ANN model. Finally, the performance of the ANN model was discussed.

    4.1 Indices for evaluating the model

    Apart from the loss function, several different comprehensive indices were used to examine the performance of the ANN model. They included the root mean squared error (RMSE), mean absolute error (MAE), and the coefficient of determination (2), which were identified and calculated by Eqs. (1)–(3). The closer2is to 1 and the smaller the values of the other errors, the better the performance of the ANN model.

    Table 2 Steel-PVA fiber hybridization of HFRC from the literature

    Table 3 Summary of the dataset for the ANN model (size=44 135)

    4.2 ANN model training and evaluation

    The database in Section 3 was used to train and evaluate the ANN model. The loss curves of the training dataset and test dataset are shown in Fig. 5a. Before 480 epochs, the loss values of the training dataset and test dataset varied greatly and showed an obvious decrease. After 480 epochs, the loss values declined gently and tended to converge.This is because the learning parameter was set as 0.01 before 480 epochs and then changed to 0.001.The lowest flat part of the curves can be observed after 960 epochs when the learning rate was changed to 0.0001. The higher learning rate could enable the model to quickly move towards a minimum of the loss function, but would also cause the loss value to oscillate easily, which is suitable for use in the initial stage. By comparison, although a smaller learning rate would slow the forward speed of the model, it could enable the model to converge faster, which is often useful after several epochs.

    The RMSE, MAE, and2values of the training and test datasets are shown in Fig. 5b. With the increase in epochs, both the RMSE and MAE of the training and test datasets decreased linearly in the initial stage, then gradually reduced until they converged after 480 epochs. The2of the training and test datasets increased quickly then remained stable after 480 epochs, with both sets having final values above 0.998. All these indices reveal the good performance of the ANN model.

    4.3 Equation-based model

    For comparison with the ANN model, an equation-based model was also designed based on the constitutive model of MHURD (2011a). The compressive stress-strain curve of HFRC can be described by

    whereis the compressive stress;is the strain;cis the elastic modulus;c,ris the compressive strength;c,ris the strain corresponding to compressive strength (SCCS);cis the damage evolution parameter, which is dependent mainly on the degree of plasticity, defined as the ratio oftoc,r;cis the descending shape parameter, which determines the slope of the descending branch during the post-peak stage;c,andare internal variables. Based on Zhou et al. (2019), thecof HFRC was assumed to be identical to that of the corresponding plain concrete in this study.

    c,r,c,r, andcare the main three parameters, and can be determined by back-fitting the experimental data.c,randc,rcan be easily obtained from the peak value of each compressive stress-strain curve;cneeds to be calculated by fitting each compressive stress-strain curve using Eq. (5). To describe the effects of steel fiber and PVA fiber, the relationship between these three parameters and the characteristic parameters of steel fiber or PVA fiber need to be formulated. In this study, a polynomial regression was used as a first approximation, as indicated in Eqs. (6)–(8).

    Fig. 5 Training and evaluation results of the ANN model: (a) loss curves; (b) RMSE, MAE, and R2

    wherecHis the compressive strength of HFRC;cHis the SCCS of HFRC;cHis the descending shape parameter of HFRC;cpis the compressive strength of plain concrete;cpis the SCCS of plain concrete;Sis the reinforcement index of steel fibers;Pis the reinforcement index of PVA fibers.S,S, andSare the length, diameter, and volume, respectively, of the steel fibers;P,P, andPare the length, diameter, and volume, respectively, of the PVA fibers.

    SandPwere introduced into Eqs. (6)–(8) to consider the hybridization of steel and PVA fibers because the hybridization effect in HFRC appeared obvious in the above experiments. The back-fitting results are illustrated in Fig. 6. Eqs. (6)–(8) reproduced the trend well, although the value of2was not big.

    As indicated in Eq. (6), the addition of PVA fibers decreased the compressive strength of HFRC, while adding a lower volume of steel fibers increased the compressive strength. However, when the fiber volume content of the steel fibers is too high (such as 1.0%), more steel fibers may result in a decrease in the compressive strength. This phenomenon is consistent with the test results in Section 2. Eq. (7) shows that adding both steel and PVA fibers increased the SCCS, and the hybridization of steel and PVA fibers had a positive influence on the SCCS, which is similar to the test results. According to Eq. (8), PVA fibers could significantly influencec, reducing the decreasing speed of stress during the post-peak stage.

    To describe the compressive stress-strain curve of HFRC,cp,cp, andcof the plain concrete were firstly obtained and used as the reference values in Eqs. (6)–(8). Secondly,cH,cH, andcHwere calculated with the fiber information using Eqs. (6)–(8). Then, the compressive stress-strain curve of HFRC was obtained using Eqs. (4) and (5) by replacingc,r,c,r, andcwith these parameters.

    Fig. 6 Back-fitting results of key parameters in the equation-based model: (a) compressive strength, fcH; (b) SCCS, εcH; (c) descending shape parameter, αcH

    4.4 Model performance analysis

    Fig. 7 presents the simulation results from the ANN model and the equation-based model, as well as the experimental results of the compressive stress-strain curve of several HFRC samples. Both the ANN model and the equation-based model reproduced the experimental data well, but the ANN model showed a better capacity. For instance, in Fig. 7a, the equation-based model could predict accurately the stress-strain curve before the peak value, but the difference during the post-peak value increased with the strain. This may be caused by the error of the descending shape parameter induced by Eq. (8); in Fig. 7d, the equation-based predicted compressive stress-strain curve is parallel with the tested curve, owing to the underestimation of the compressive strength by Eq. (6). However, for these two cases, the ANN model reproduced the stress-strain curve well, and almost overlapped the tested curves.

    Table 4 summarizes the RMSE, MAE, and2values of the ANN model and the equation-based model. The ANN model had lower errors and a higher2than the equation-based model. Therefore, we conclude that the ANN model can perform better than the equation-based model to reproduce the compressive stress-strain curve.

    In this paper, we also introduce the relative error (Zhou et al., 2019) to evaluate the performance of the ANN model and equation-based model in terms of predicting the compressive strength and SCCS. The relative error is calculated by

    4.5 Sensitivity analysis

    Since the ANN model is implicit, it is difficult to understand the relationship between input features and output features. Hence, a sensitivity analysis has been proposed to illustrate how input features affect output features for an ANN model (Saltelli, 2002; Shojaeefard et al., 2013). Among several sensitivity analysis methods, the use of the “PaD” method (partial derivatives) makes it easy to calculate the

    Table 4 RMSE, MAE, and R2 of the ANN model and the equation-based model

    Fig. 8 Relative errors of the ANN model and the equation-based model in terms of the compressive strength (a) and SCCS (b)

    Fig. 9 Contributions of each input feature to the compressive strength (a) and the strain corresponding to compressive strength of HFRC (b)

    Fig. 9 shows the contribution of each input feature to the compressive strength and the SCCS for HFRC samples. The redder the color, the higher the contribution of the specified input feature. The highest contribution appears to be the compressive strength and the SCCS of plain concrete. There are four main influence parts, including steel fibers, PVA fibers, mechanical properties of plain concrete, and components of HFRC. Fig. 10 reveals the contribution of these four main parts to the compressive strength and the SCCS, which is the sum of their sub-input features. According to Figs. 9 and 10, the mechanical properties of plain concrete make the main contribution to the compressive strength and the SCCS for all HFRC samples. The contribution of steel fibers is higher than that of PVA fibers in terms of the compressive strength and the SCCS for all HFRC samples. This is similar to the experimental results.

    Fig. 10 Contribution of steel fibers, PVA fibers, mechanical properties of plain concrete, and components of HFRC to the compressive strength (a) and the SCCS (b)

    4.6 Discussion

    According to the analysis above, it is clear that the ANN model developed in this study is better able to describe the compressive behavior of HFRC than the equation-based model. The performance of the equation-based model depends on the determination of key parameters, including the compressive strength, SCCS, and descending shape parameter, calculated using Eqs. (6)–(8). If one of these parameters is not well determined, the predicted compressive stress-strain curve would not fit well with the experimental results. Also, only limited factors are considered in Eqs. (6)–(8), and many other components, such as the cementitious material, that could significantly influence the mechanical properties of concrete, are not directly considered in these equations. The2value of Eqs. (6)–(8) seems relatively low, indicating that the equation-based model could reproduce well the stress-strain curve for some HFRC samples while having obvious error for others, as shown in Figs. 6 and 7. However, the equation-based model is simple and can be easily used in the calculation of the structural components and structures, resulting in the wide application of such types of models. In summary, the equation-based model needs further development, especially the methodology to consider the effects of different compositions.

    By comparison, the ANN model can consider more factors. Thus, it is able to extract more important features and has more powerful nonlinear fitting abilities, resulting in a better capacity to predict the compressive stress-strain curve of HFRC. Moreover, the ANN model could be applied not only to HFRC, but also to other composite materials (such as ultra-high-performance concrete) and even material that contains different composite components, which are difficult to model by an equation-based model. On the other hand, the ANN model needs a large database to improve its performance, especially when it is used to predict the stress-strain curve. To consider these factors, more and credible databases are necessary, and this work is being carried out by our group. Furthermore, the geometry of steel fibers and PVA fibers may also affect the compressive behavior, but cannot be considered by the current ANN model. However, the influence of these factors could be clarified via the ANN model when the diversity of the database is big enough. Also, the training of the ANN model needs time. For example, about 1 h was necessary for training the ANN model in this study based on GeForce RTX 2070 SUPER GPU. When adding more datasets, the model needs a longer time for training to improve its performance. However, these problems may be solved by developments in computer science. As a first approximation, the presented work demonstrates the suitability of the ANN model for describing the compressive behavior of HFRC.

    In summary, the choice of an equation-based model or an ANN model should be related to actual needs. The ANN model is more suitable for the accurate design of composite material, such as HFRC, while the equation-based model can be easily used in structural analysis.

    5 Conclusions

    This paper presents a comprehensive investigation on the compressive behavior of HRFC containing fly ash and slag powder via experiments and modelling. A series of experiments were carried out to investigate the hybrid effect of steel and PVA fiber on the compressive behavior of HFRC containing fly ash and slag powder. Meanwhile, an ANN model to describe the compressive behavior was also built and trained with the experimental data and relevant data from the literature. The performance of the ANN model was compared to an equation-based model. Based on the above work, the following conclusions can be drawn:

    1. The proposed ANN model has a better capacity than the equation-based model to reproduce the compressive behavior of hybrid steel-PVA fiber concrete containing fly ash and slag powder in terms of the compressive stress-strain curve, the compressive strength, and the SCCS. The ANN model can consider many factors affecting the mechanical behavior of HFRC. In this study, 23 factors were considered, including features of fibers, concrete composition, and mechanical behavior of plain concrete. Note that more factors may be included when the database is large enough. Therefore, the ANN model is a good tool for designing composite materials, and is not limited to HFRC.

    2. Both the experimental results and the ANN model simulation results revealed that steel fibers play a more important role in increasing the compressive strength of HFRC than PVA fibers, while both contribute to the deformation. Increasing the steel fiber volume could increase the compressive strength, while increasing the PVA fiber volume could reduce it. However, when the steel fiber volume is too high (such as above 1.0%), the compressive strength tends to decrease with increasing steel fiber volume. The hybridization effect of steel and PVA fibers could enable HFRC to shift from a brittle break to a ductile break.

    Contributors

    Fang-yu LIU, Wen-qi DING, and Ya-fei QIAO designed the research. Fang-yu LIU,Ya-fei QIAO, and Qi-yang CHEN processed the corresponding data. Fang-yu LIU and Ya-fei QIAO wrote the first draft of the manuscript. Wen-qi DING and Lin-bing WANG helped to organize the manuscript. Fang-yu LIU, Wen-qi DING, Ya-fei QIAO, and Lin-bing WANG revised and edited the final version.

    Conflict of interest

    Fang-yu LIU, Wen-qi DING, Ya-fei QIAO, Lin-bing WANG, and Qi-yang CHEN declare that they have no conflict of interest.

    List of electronic supplementary materials

    Data S1 Framework of the ANN model

    A?ikgen? M, Ula? M, Alyama? KE, 2015. Using an artificial neural network to predict mix compositions of steel fiber-reinforced concrete., 40(2):407-419. https://doi.org/10.1007/s13369-014-1549-x

    Arulmurugan P, 2018. Evaluation of FRC beams using steel and pva fibres in concrete., 6(1):121-124.

    Bai J, Wild S, Ware JA, et al., 2003. Using neural networks to predict workability of concrete incorporating metakaolin and fly ash., 34(11-12): 663-669. https://doi.org/10.1016/S0965-9978(03)00102-9

    Belletti B, Cerioni R, Meda A, et al., 2008. Design aspects on steel fiber-reinforced concrete pavements., 20(9):599-607. https://doi.org/10.1061/(asce)0899-1561(2008)20:9(599)

    Blanco A, Pujadas P, de la Fuente A, et al., 2016. Influence of the type of fiber on the structural response and design of FRC slabs., 142(9): 04016054. https://doi.org/10.1061/(ASCE)ST.1943-541X.0001515

    CAECS (China Association for Engineering Construction Standardization), 2010. Standard Test Methods for Fiber Reinforced Concrete, CECS 13:2009. National Standards of the People’s Republic of China (in Chinese).

    Chi Y, Xu LH, Yu HS, 2014a. Constitutive modeling of steel-polypropylene hybrid fiber reinforced concrete using a non-associated plasticity and its numerical implementation., 111:497-509. https://doi.org/10.1016/j.compstruct.2014.01.025

    Chi Y, Xu LH, Zhang YY, 2014b. Experimental study on hybrid fiber-reinforced concrete subjected to uniaxial compression., 26(2):211-218. https://doi.org/10.1061/(ASCE)MT.1943-5533.0000764

    de la Fuente A, Pujadas P, Blanco A, et al., 2012. Experiences in Barcelona with the use of fibres in segmental linings., 27(1): 60-71. https://doi.org/10.1016/j.tust.2011.07.001

    Ding WQ, Gong YF, Qiao YF, et al., 2020. Experimental investigation on mechanical behavior of segmental joint under combined loading of compression-bending-shear., 98: 103346. https://doi.org/10.1016/j.tust.2020.103346

    Feng J, Sun WW, Zhai HZ, et al., 2018. Experimental study on hybrid effect evaluation of fiber reinforced concrete subjected to drop weight impacts., 11(12): 2563. https://doi.org/10.3390/ma11122563

    Gong CJ, Ding WQ, Mosalam KM, et al., 2017. Comparison of the structural behavior of reinforced concrete and steel fiber reinforced concrete tunnel segmental joints., 68: 38-57. https://doi.org/10.1016/j.tust.2017.05.010

    Hai R, Liu JX, Zhang ML, et al., 2016. Performance of hybrid steel-polyvinyl alcohol fiber reinforced ultra high performance concrete., (5):95-97 (in Chinese). https://doi.org/10.3969/j.issn.1002-3550.2016.05.025

    Hossain KMA, Lachemi M, Sammour M, et al., 2013. Strength and fracture energy characteristics of self-consolidating concrete incorporating polyvinyl alcohol, steel and hybrid fibres., 45:20-29. https://doi.org/10.1016/j.conbuildmat.2013.03.054

    Jiang KJ, Han Q, Bai YL, et al., 2020. Data-driven ultimate conditions prediction and stress-strain model for FRP-confined concrete., 242:112094. https://doi.org/10.1016/j.compstruct.2020.112094

    Karahan O, Tanyildizi H, Atis CD, 2008. An artificial neural network approach for prediction of long-term strength properties of steel fiber reinforced concrete containing fly ash., 9:1514- 1523.https://doi.org/10.1631/jzus.A0720136

    Lawler JS, Wilhelm T, Zampini D, et al., 2003. Fracture processes of hybrid fiber-reinforced mortar., 36(3):197-208. https://doi.org/10.1007/BF02479558

    Lawler JS, Zampini D, Shah SP, 2005. Microfiber and macrofiber hybrid fiber-reinforced concrete., 17(5):595-604. https://doi.org/10.1061/(asce)0899-1561(2005)17:5(595)

    Li HY, Liu G, 2016. Tensile properties of hybrid fiber-reinforced reactive powder concrete after exposure to elevated temperatures., 10(1):29-37. https://doi.org/10.1007/s40069-016-0125-z

    Liao L, de la Fuente A, Cavalaro S, et al., 2015. Experimental and analytical study of concrete blocks subjected to concentrated loads with an application to TBM-constructed tunnels., 49:295-306. https://doi.org/10.1016/j.tust.2015.04.020

    Liu FY, Ding WQ, Qiao YF, 2019a. An experimental investigation on the integral waterproofing capacity of polypropylene fiber concrete with fly ash and slag powder., 212:675-686. https://doi.org/10.1016/j.conbuildmat.2019.04.027

    Liu FY, Ding WQ, Qiao YF, 2019b. Experimental investigation on the flexural behavior of hybrid steel-PVA fiber reinforced concrete containing fly ash and slag powder., 228:116706. https://doi.org/10.1016/j.conbuildmat.2019.116706

    Liu FY, Ding WQ, Qiao YF, et al., 2020a. An artificial neural network model on tensile behavior of hybrid steel-PVA fiber reinforced concrete containing fly ash and slag power., 14: 1299-1315.https://doi.org/10.1007/s11709-020-0712-6

    Liu FY, Ding WQ, Qiao YF, 2020b. Experimental investigation on the tensile behavior of hybrid steel-PVA fiber reinforced concrete containing fly ash and slag powder., 241:118000. https://doi.org/10.1016/j.conbuildmat.2020.118000

    Mashhadban H, Kutanaei SS, Sayarinejad MA, 2016. Prediction and modeling of mechanical properties in fiber reinforced self-compacting concrete using particle swarm optimization algorithm and artificial neural network., 119:277-287. https://doi.org/10.1016/j.conbuildmat.2016.05.034

    MHURD (Ministry of Housing and Urban-Rural Development of the People’s Republic of China), 2011a. Code for Design of Concrete Structure, GB 50010-2010. National Standards of People’s Republic of China (in Chinese).

    MHURD (Ministry of Housing and Urban-Rural Development of the People’s the Republic of China), 2011b. Specification for Mix Proportion Design of Ordinary Concrete, JGJ 55-2011. National Standards of the People’s Republic of China (in Chinese).

    Nguyen DL, Kim DJ, Ryu GS, et al., 2013. Size effect on flexural behavior of ultra-high-performance hybrid fiber-reinforced concrete., 45(1):1104-1116. https://doi.org/10.1016/j.compositesb.2012.07.012

    Pizarroso J, Portela J, Mu?oz A, 2020. NeuralSens: sensitivity analysis of neural networks. arXiv: 2002.11423.https://arxiv.org/abs/2002.11423

    Pujadas P, Blanco A, Cavalaro S, et al., 2014. Plastic fibres as the only reinforcement for flat suspended slabs: experimental investigation and numerical simulation., 57:92-104. https://doi.org/10.1016/j.conbuildmat.2014.01.082

    Pujadas P, Blanco A, Cavalaro S, et al., 2017. The need to consider flexural post-cracking creep behavior of macro- synthetic fiber reinforced concrete., 149:790-800. https://doi.org/10.1016/j.conbuildmat.2017.05.166

    Saltelli A, 2002. Sensitivity analysis for importance assessment., 22(3):579-590. https://doi.org/10.1111/0272-4332.00040

    Serna Ros P, Martí-Vargas JR, Bossio ME, et al., 2016. Creep and residual properties of cracked macro-synthetic fibre reinforced concretes., 68(4):197-207. https://doi.org/10.1680/macr.15.00111

    Shen QZ, 2013. Mechanical Properties Research for Steel-PVA Hybrid Fibers Reinforced Concrete. MS Thesis, Chang’an University, Xi’an, China (in Chinese).

    Shojaeefard MH, Akbari M, Tahani M, et al., 2013. Sensitivity analysis of the artificial neural network outputs in friction stir lap joining of aluminum to brass., 2013:574914. https://doi.org/10.1155/2013/574914

    Sim JS, Park C, Moon DY, 2005. Characteristics of basalt fiber as a strengthening material for concrete structures., 36(6-7):504-512. https://doi.org/10.1016/j.compositesb.2005.02.002

    Soutsos MN, Le TT, Lampropoulos AP, 2012. Flexural performance of fibre reinforced concrete made with steel and synthetic fibres., 36:704-710. https://doi.org/10.1016/j.conbuildmat.2012.06.042

    Stoll F, Saliba JE, Casper LE, 2000. Experimental study of CFRP-prestressed high-strength concrete bridge beams., 49(2):191-200. https://doi.org/10.1016/S0263-8223(99)00134-8

    Tabatabaeian M, Khaloo A, Joshaghani A, et al., 2017. Experimental investigation on effects of hybrid fibers on rheological, mechanical, and durability properties of high-strength SCC., 147:497-509. https://doi.org/10.1016/j.conbuildmat.2017.04.181

    Xu LH, Xu HR, Chi Y, et al., 2011. Experimental study on tensile strength of steel-polypropylene hybrid fiber reinforced concrete., 4(3):911- 916. https://doi.org/10.1166/asl.2011.1740

    Yang KH, 2011. Tests on concrete reinforced with hybrid or monolithic steel and polyvinyl alcohol fibers., 108(6):664-672.

    Zhao GY, di Prisco M, Vandewalle L, 2015. Experimental investigation on uniaxial tensile creep behavior of cracked steel fiber reinforced concrete., 48(10):3173-3185. https://doi.org/10.1617/s11527-014-0389-1

    Zhou Y, Xiao Y, Gu AQ, et al., 2018. Dispersion, workability and mechanical properties of different steel-microfiber-reinforced concretes with low fiber content., 10(7):2335. https://doi.org/10.3390/su10072335

    Zhou Y, Xiao Y, Gu AQ, et al., 2019. Orthogonal experimental investigation of steel-PVA fiber-reinforced concrete and its uniaxial constitutive model., 197:615-625. https://doi.org/10.1016/j.conbuildmat.2018.11.203

    Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering)

    ISSN 1673-565X (Print); ISSN 1862-1775 (Online)

    www.jzus.zju.edu.cn; link.springer.com

    E-mail: jzus_a@zju.edu.cn

    https://doi.org/10.1631/jzus.A2000379

    # Electronic supplementary materials: The online version of this article (https://doi.org/10.1631/jzus.A2000379) contains supplementary materials, which are available to authorized users

    ? Zhejiang University Press 2021

    TU528.572

    Aug. 23, 2020;

    Revision accepted Nov. 11, 2020;

    Crosschecked Sept. 3, 2021

    *Project supported by the National Natural Science Foundation of China (Nos. 51978515 and 52090083), the Shanghai Sailing Program (No. 19YF1451400), and the Shanghai Municipal Science and Technology Major Project (No. 2017SHZDZX02), China

    日韩人妻精品一区2区三区| 亚洲精品美女久久久久99蜜臀 | 777米奇影视久久| 成人毛片60女人毛片免费| 80岁老熟妇乱子伦牲交| 国产成人欧美| 无限看片的www在线观看| 国产日韩一区二区三区精品不卡| 操出白浆在线播放| 51午夜福利影视在线观看| 免费看不卡的av| 国产乱人偷精品视频| 老司机影院成人| 啦啦啦中文免费视频观看日本| 十分钟在线观看高清视频www| 日韩人妻精品一区2区三区| 黄色一级大片看看| 中文字幕制服av| 亚洲四区av| 黄色怎么调成土黄色| 亚洲精品乱久久久久久| 香蕉国产在线看| 叶爱在线成人免费视频播放| 国产一卡二卡三卡精品 | av在线老鸭窝| 久久久久久免费高清国产稀缺| 亚洲欧美清纯卡通| 精品一区二区三区av网在线观看 | 在线天堂中文资源库| 亚洲人成77777在线视频| 国产精品久久久久久精品古装| 看免费av毛片| 少妇人妻久久综合中文| 午夜福利网站1000一区二区三区| h视频一区二区三区| 老司机亚洲免费影院| 777米奇影视久久| 中文字幕人妻丝袜制服| 日韩 欧美 亚洲 中文字幕| 十八禁高潮呻吟视频| 亚洲精品久久午夜乱码| 少妇人妻 视频| 叶爱在线成人免费视频播放| 国产日韩欧美在线精品| 在线 av 中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 一本大道久久a久久精品| 中国三级夫妇交换| 免费观看人在逋| 美女中出高潮动态图| 免费黄色在线免费观看| 国产精品三级大全| 1024视频免费在线观看| 亚洲久久久国产精品| 99久久人妻综合| 一边摸一边做爽爽视频免费| 一区二区三区乱码不卡18| 高清不卡的av网站| 在线观看免费高清a一片| 色视频在线一区二区三区| 国产亚洲最大av| xxx大片免费视频| 最近最新中文字幕免费大全7| av福利片在线| 最黄视频免费看| 色婷婷av一区二区三区视频| 成人手机av| 视频在线观看一区二区三区| 在线看a的网站| 最近中文字幕2019免费版| 伦理电影大哥的女人| 中文字幕人妻丝袜一区二区 | 最近最新中文字幕免费大全7| 欧美日韩一级在线毛片| 只有这里有精品99| 国产黄色免费在线视频| 亚洲一区二区三区欧美精品| 电影成人av| 欧美精品人与动牲交sv欧美| 精品国产一区二区三区久久久樱花| 日韩熟女老妇一区二区性免费视频| 日本黄色日本黄色录像| 国产精品二区激情视频| 亚洲视频免费观看视频| a级片在线免费高清观看视频| 水蜜桃什么品种好| 久久精品国产综合久久久| 18禁裸乳无遮挡动漫免费视频| 99热国产这里只有精品6| 最新的欧美精品一区二区| 欧美精品一区二区大全| 你懂的网址亚洲精品在线观看| 亚洲中文av在线| 在线天堂中文资源库| 免费看不卡的av| 一区二区日韩欧美中文字幕| 免费少妇av软件| 伊人久久国产一区二区| 在线亚洲精品国产二区图片欧美| 超碰97精品在线观看| 亚洲欧美精品综合一区二区三区| www日本在线高清视频| 一本久久精品| 精品卡一卡二卡四卡免费| 黄色怎么调成土黄色| 9色porny在线观看| 国产高清国产精品国产三级| 国产片内射在线| 欧美中文综合在线视频| 亚洲av福利一区| 操美女的视频在线观看| 亚洲国产av影院在线观看| 成人国产麻豆网| 免费在线观看黄色视频的| 99精品久久久久人妻精品| 别揉我奶头~嗯~啊~动态视频 | 男人爽女人下面视频在线观看| 9色porny在线观看| 久久这里只有精品19| 又大又黄又爽视频免费| 欧美精品一区二区免费开放| 国产精品偷伦视频观看了| 国产精品偷伦视频观看了| 国产黄频视频在线观看| 日韩大码丰满熟妇| 日韩一本色道免费dvd| 久久精品久久久久久噜噜老黄| 日本wwww免费看| 美女主播在线视频| 精品国产乱码久久久久久小说| 日韩伦理黄色片| 亚洲欧洲国产日韩| 日韩不卡一区二区三区视频在线| 亚洲av成人精品一二三区| 亚洲国产av影院在线观看| 性少妇av在线| 国产97色在线日韩免费| 最新在线观看一区二区三区 | 岛国毛片在线播放| a级片在线免费高清观看视频| 亚洲,欧美精品.| 亚洲,欧美精品.| 国产精品二区激情视频| 欧美日本中文国产一区发布| 伦理电影大哥的女人| 90打野战视频偷拍视频| 国产精品一区二区精品视频观看| 观看美女的网站| 啦啦啦中文免费视频观看日本| 丁香六月欧美| 久久国产亚洲av麻豆专区| 久久精品国产综合久久久| 99久久99久久久精品蜜桃| 久久韩国三级中文字幕| 欧美乱码精品一区二区三区| 亚洲四区av| 日本欧美国产在线视频| 免费在线观看视频国产中文字幕亚洲 | 久久av网站| 91精品三级在线观看| 女人久久www免费人成看片| 久久国产精品大桥未久av| 精品国产乱码久久久久久男人| 国产免费福利视频在线观看| 久久久久久久久久久久大奶| 亚洲一区二区三区欧美精品| 香蕉国产在线看| 久久久久精品久久久久真实原创| 精品一区二区三区四区五区乱码 | 国产在线免费精品| 色网站视频免费| 在线看a的网站| 汤姆久久久久久久影院中文字幕| 国产一区二区三区av在线| 中文字幕人妻熟女乱码| 亚洲伊人久久精品综合| xxx大片免费视频| 国产人伦9x9x在线观看| 一区二区三区四区激情视频| 国产免费一区二区三区四区乱码| 国产精品熟女久久久久浪| 亚洲国产精品成人久久小说| 麻豆精品久久久久久蜜桃| 婷婷成人精品国产| 韩国精品一区二区三区| 午夜福利视频在线观看免费| 在线观看免费视频网站a站| av视频免费观看在线观看| 午夜福利,免费看| 成人免费观看视频高清| 汤姆久久久久久久影院中文字幕| 久久人妻熟女aⅴ| 国产成人精品久久久久久| 在线看a的网站| 满18在线观看网站| 精品免费久久久久久久清纯 | 免费久久久久久久精品成人欧美视频| 精品一区二区三区四区五区乱码 | 午夜福利乱码中文字幕| 欧美97在线视频| 丝袜喷水一区| 欧美在线黄色| 精品国产一区二区三区四区第35| 欧美97在线视频| 午夜福利一区二区在线看| 国产av一区二区精品久久| 男女无遮挡免费网站观看| 美女中出高潮动态图| 精品一区二区三卡| 日韩 欧美 亚洲 中文字幕| h视频一区二区三区| 精品人妻在线不人妻| 久久久欧美国产精品| 亚洲第一青青草原| 欧美黄色片欧美黄色片| 亚洲久久久国产精品| 看免费成人av毛片| 久久精品人人爽人人爽视色| 国产免费现黄频在线看| 波野结衣二区三区在线| 99香蕉大伊视频| av网站在线播放免费| 99久久综合免费| www日本在线高清视频| 高清在线视频一区二区三区| 欧美黑人精品巨大| videosex国产| 日韩人妻精品一区2区三区| 久久精品久久久久久噜噜老黄| 美女扒开内裤让男人捅视频| 国产成人免费无遮挡视频| 中文乱码字字幕精品一区二区三区| 三上悠亚av全集在线观看| 人妻 亚洲 视频| 人体艺术视频欧美日本| 大香蕉久久成人网| 日韩大片免费观看网站| 欧美久久黑人一区二区| 精品国产国语对白av| 侵犯人妻中文字幕一二三四区| 十八禁高潮呻吟视频| 久久亚洲国产成人精品v| 中文字幕人妻丝袜制服| av免费观看日本| 2018国产大陆天天弄谢| 91精品国产国语对白视频| 侵犯人妻中文字幕一二三四区| 久久人人97超碰香蕉20202| 国产不卡av网站在线观看| 一级毛片黄色毛片免费观看视频| 最新的欧美精品一区二区| 日韩大片免费观看网站| 色视频在线一区二区三区| av女优亚洲男人天堂| 另类精品久久| 国产深夜福利视频在线观看| 宅男免费午夜| av.在线天堂| 人人妻人人澡人人看| 亚洲色图综合在线观看| 国产午夜精品一二区理论片| 一级,二级,三级黄色视频| 国产精品欧美亚洲77777| 亚洲国产精品成人久久小说| 国产一区有黄有色的免费视频| 你懂的网址亚洲精品在线观看| 色吧在线观看| 亚洲欧美一区二区三区久久| 自线自在国产av| 夫妻性生交免费视频一级片| 中文字幕人妻熟女乱码| 校园人妻丝袜中文字幕| 亚洲国产欧美日韩在线播放| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频 | 婷婷色av中文字幕| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品一区在线观看| 大香蕉久久成人网| 欧美日韩亚洲国产一区二区在线观看 | 国产人伦9x9x在线观看| 老司机深夜福利视频在线观看 | 亚洲精品,欧美精品| 大码成人一级视频| 午夜日韩欧美国产| 999精品在线视频| 国产午夜精品一二区理论片| 黑丝袜美女国产一区| av免费观看日本| 中文字幕最新亚洲高清| 亚洲精品日本国产第一区| 亚洲精品一区蜜桃| 欧美老熟妇乱子伦牲交| 纵有疾风起免费观看全集完整版| 亚洲美女视频黄频| 亚洲精品久久午夜乱码| 一级黄片播放器| 欧美av亚洲av综合av国产av | 国产高清国产精品国产三级| 久久久久网色| 国产成人系列免费观看| 亚洲成人免费av在线播放| 国产欧美亚洲国产| 只有这里有精品99| 亚洲精品国产一区二区精华液| 亚洲欧美清纯卡通| 黑人巨大精品欧美一区二区蜜桃| 国产伦理片在线播放av一区| 亚洲国产欧美一区二区综合| 亚洲成国产人片在线观看| 黄色毛片三级朝国网站| 国产乱来视频区| 一区二区日韩欧美中文字幕| 精品国产乱码久久久久久男人| 亚洲欧美成人综合另类久久久| 91老司机精品| 亚洲精品第二区| 精品人妻在线不人妻| 国产成人精品无人区| 女人精品久久久久毛片| 国产成人a∨麻豆精品| 亚洲国产欧美日韩在线播放| 久久av网站| 欧美亚洲 丝袜 人妻 在线| 熟女少妇亚洲综合色aaa.| 不卡av一区二区三区| 极品少妇高潮喷水抽搐| 女人久久www免费人成看片| 波多野结衣一区麻豆| 久久人人爽av亚洲精品天堂| 精品少妇一区二区三区视频日本电影 | 国产精品蜜桃在线观看| 欧美黑人欧美精品刺激| 国产精品偷伦视频观看了| 午夜91福利影院| 亚洲伊人色综图| 国产一区有黄有色的免费视频| 超碰成人久久| 日韩中文字幕欧美一区二区 | 一个人免费看片子| 国产熟女欧美一区二区| 成人国语在线视频| 国语对白做爰xxxⅹ性视频网站| 伊人久久国产一区二区| 久久精品人人爽人人爽视色| 亚洲欧洲国产日韩| 中文字幕高清在线视频| 国产爽快片一区二区三区| 亚洲中文av在线| 久久久久精品久久久久真实原创| 成人漫画全彩无遮挡| 狂野欧美激情性bbbbbb| 国产色婷婷99| 交换朋友夫妻互换小说| 国产成人啪精品午夜网站| 水蜜桃什么品种好| 考比视频在线观看| 久久免费观看电影| 男女高潮啪啪啪动态图| 99热全是精品| 免费在线观看完整版高清| 国产一区二区三区av在线| 校园人妻丝袜中文字幕| 如何舔出高潮| av有码第一页| 精品久久久久久电影网| 欧美乱码精品一区二区三区| 午夜福利视频在线观看免费| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲自偷自拍图片 自拍| 欧美精品高潮呻吟av久久| 美女视频免费永久观看网站| 日韩精品免费视频一区二区三区| 亚洲国产欧美一区二区综合| 男女国产视频网站| 一边摸一边做爽爽视频免费| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 国精品久久久久久国模美| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区 | 自线自在国产av| 老司机影院毛片| 国产一区亚洲一区在线观看| 久久久久精品性色| 亚洲成人国产一区在线观看 | av在线观看视频网站免费| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久久久免| 色94色欧美一区二区| 亚洲免费av在线视频| 五月天丁香电影| 久久久欧美国产精品| 日本色播在线视频| 99久久综合免费| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 亚洲国产精品一区三区| 如何舔出高潮| 亚洲成色77777| 亚洲精品乱久久久久久| 国产精品免费视频内射| 在线精品无人区一区二区三| 少妇被粗大猛烈的视频| 高清欧美精品videossex| 少妇被粗大的猛进出69影院| 国产一区二区三区av在线| 亚洲成人手机| 一级爰片在线观看| 色精品久久人妻99蜜桃| 黄色怎么调成土黄色| 精品亚洲乱码少妇综合久久| 51午夜福利影视在线观看| 黑人猛操日本美女一级片| 亚洲精品中文字幕在线视频| 国产成人精品在线电影| 亚洲综合色网址| 青春草国产在线视频| 国产精品.久久久| 亚洲国产欧美在线一区| 亚洲欧美激情在线| 十八禁网站网址无遮挡| 欧美日韩一级在线毛片| 免费在线观看视频国产中文字幕亚洲 | 亚洲av电影在线进入| 中国国产av一级| 亚洲av综合色区一区| 不卡视频在线观看欧美| 欧美日韩福利视频一区二区| 99九九在线精品视频| 一级爰片在线观看| 天天躁夜夜躁狠狠久久av| 两个人看的免费小视频| 中文字幕亚洲精品专区| 久久久久久久久久久免费av| 国产成人免费观看mmmm| 久久精品亚洲熟妇少妇任你| 国产免费福利视频在线观看| 亚洲av国产av综合av卡| 天堂8中文在线网| 在线观看人妻少妇| 中文字幕人妻熟女乱码| 亚洲欧美日韩另类电影网站| 人人妻人人爽人人添夜夜欢视频| 欧美日韩一级在线毛片| 女人爽到高潮嗷嗷叫在线视频| 国产色婷婷99| 又大又爽又粗| 中文字幕av电影在线播放| 97在线人人人人妻| 肉色欧美久久久久久久蜜桃| 国产精品人妻久久久影院| 欧美日韩精品网址| 99re6热这里在线精品视频| 看非洲黑人一级黄片| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 中文天堂在线官网| 国产成人av激情在线播放| 韩国精品一区二区三区| 久久精品国产亚洲av高清一级| 伦理电影免费视频| 亚洲精品国产av蜜桃| 伊人亚洲综合成人网| 免费在线观看视频国产中文字幕亚洲 | 国产精品 国内视频| 精品少妇内射三级| 亚洲欧美成人精品一区二区| 中国三级夫妇交换| 丝袜在线中文字幕| 一级片免费观看大全| 欧美乱码精品一区二区三区| 国产精品二区激情视频| 亚洲综合色网址| www.自偷自拍.com| 国产精品偷伦视频观看了| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 成人三级做爰电影| 18禁观看日本| 视频在线观看一区二区三区| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 国产男人的电影天堂91| 纯流量卡能插随身wifi吗| 电影成人av| 日韩欧美精品免费久久| 亚洲在久久综合| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| 婷婷色麻豆天堂久久| 天天躁夜夜躁狠狠久久av| 大码成人一级视频| 飞空精品影院首页| 一级片免费观看大全| 亚洲国产最新在线播放| 欧美人与善性xxx| 日本猛色少妇xxxxx猛交久久| 乱人伦中国视频| 亚洲美女黄色视频免费看| 欧美日韩av久久| 国产高清不卡午夜福利| 欧美久久黑人一区二区| 青春草视频在线免费观看| 日韩不卡一区二区三区视频在线| 亚洲熟女毛片儿| 男女免费视频国产| 国产免费一区二区三区四区乱码| 人人妻人人爽人人添夜夜欢视频| 一级片'在线观看视频| 国产乱人偷精品视频| 国产精品一国产av| 黄色视频在线播放观看不卡| 人体艺术视频欧美日本| 午夜日韩欧美国产| 黄色一级大片看看| 亚洲 欧美一区二区三区| 尾随美女入室| 男女高潮啪啪啪动态图| avwww免费| 在线观看三级黄色| 国产又爽黄色视频| 亚洲免费av在线视频| 久久久久久人妻| 国产精品成人在线| 少妇人妻精品综合一区二区| av.在线天堂| 亚洲国产欧美网| 国产在线视频一区二区| 久久ye,这里只有精品| 巨乳人妻的诱惑在线观看| 精品久久久久久电影网| 欧美 日韩 精品 国产| 只有这里有精品99| 水蜜桃什么品种好| 欧美xxⅹ黑人| 欧美乱码精品一区二区三区| 看免费成人av毛片| 精品一区二区三区av网在线观看 | 免费在线观看视频国产中文字幕亚洲 | 18在线观看网站| 久久久久精品性色| 亚洲第一区二区三区不卡| 亚洲国产欧美网| 大香蕉久久网| 日本爱情动作片www.在线观看| a级片在线免费高清观看视频| 午夜久久久在线观看| 搡老岳熟女国产| 日韩大片免费观看网站| 国产在线一区二区三区精| 国产欧美日韩综合在线一区二区| 免费看av在线观看网站| 又大又爽又粗| 一级a爱视频在线免费观看| 国产女主播在线喷水免费视频网站| 夜夜骑夜夜射夜夜干| 天天添夜夜摸| 色播在线永久视频| 久久久久久久久久久免费av| av国产精品久久久久影院| 在线精品无人区一区二区三| 精品第一国产精品| 亚洲精品一二三| svipshipincom国产片| 80岁老熟妇乱子伦牲交| 日本猛色少妇xxxxx猛交久久| 精品国产一区二区三区四区第35| 久久久久视频综合| tube8黄色片| 一本色道久久久久久精品综合| 亚洲少妇的诱惑av| 久久这里只有精品19| 韩国精品一区二区三区| 亚洲av成人精品一二三区| 午夜精品国产一区二区电影| 美女主播在线视频| 纯流量卡能插随身wifi吗| 美女主播在线视频| 中文字幕制服av| 国产欧美日韩综合在线一区二区| 中文字幕制服av| bbb黄色大片| 99久久人妻综合| 日韩一卡2卡3卡4卡2021年| 乱人伦中国视频| 99热网站在线观看| 午夜福利视频在线观看免费| 国产成人欧美在线观看 | 中国三级夫妇交换| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 美女午夜性视频免费| 日韩av免费高清视频| 亚洲精品久久成人aⅴ小说| 欧美少妇被猛烈插入视频| 午夜免费男女啪啪视频观看| 久久午夜综合久久蜜桃| 高清视频免费观看一区二区| 我的亚洲天堂| 黄色视频不卡| 大片免费播放器 马上看| 人体艺术视频欧美日本| 亚洲精品日本国产第一区| 亚洲精品日韩在线中文字幕| 精品国产超薄肉色丝袜足j| 最黄视频免费看| 999久久久国产精品视频| 精品久久久久久电影网| 免费黄频网站在线观看国产| 久久久久国产一级毛片高清牌| 国产精品香港三级国产av潘金莲 | 国产成人欧美在线观看 | 在线观看免费午夜福利视频|