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

    高強(qiáng)Q690鋼柱高溫下軸心受壓局部穩(wěn)定設(shè)計(jì)方法

    2021-03-11 10:16王衛(wèi)永SHARHANAhmedAL-AZZANIHisham李翔
    土木建筑與環(huán)境工程 2021年1期
    關(guān)鍵詞:鋼柱有限元分析高溫

    王衛(wèi)永 SHARHAN Ahmed AL-AZZANI Hisham 李翔

    摘 要:局部屈曲是鋼結(jié)構(gòu)構(gòu)件的一種破壞模式,鋼結(jié)構(gòu)發(fā)生局部屈曲破壞時(shí),屈曲應(yīng)力小于鋼材的屈服強(qiáng)度。為了研究高溫下高強(qiáng)Q690鋼柱的局部穩(wěn)定性能,采用有限元軟件ABAQUS建立有限元模型,模型采用其他學(xué)者完成的Q460鋼柱軸心受壓局部屈曲試驗(yàn)進(jìn)行驗(yàn)證,考慮寬厚比、溫度、初始缺陷、殘余應(yīng)力和翼緣與腹板之間相互作用的影響,對(duì)高強(qiáng)Q690鋼柱進(jìn)行參數(shù)分析。研究結(jié)果表明:寬厚比對(duì)局部屈曲有顯著影響,寬厚比的增大導(dǎo)致試件極限承載力的降低;初始缺陷和殘余應(yīng)力對(duì)局部屈曲應(yīng)力有較大影響,且試件的極限承載力隨著溫度的升高而明顯下降。基于有限元分析結(jié)果提出了適用于高強(qiáng)Q690鋼柱高溫下的局部穩(wěn)定設(shè)計(jì)方法和寬厚比限值,并與GB 50017-2017、Eurocode 3和ANSI /AISC 360-10中的設(shè)計(jì)方法進(jìn)行了比較。

    關(guān)鍵詞:高強(qiáng)度鋼;鋼柱;局部屈曲;高溫;有限元分析

    1 Introduction

    High strength steel (HSS) has been widely adopted in the construction of high-rise structures such as high-rise buildings and long-span bridges on account of its light-weight, with the result that these structures have become novel in shape and complicated in function. In this era of rapid development, steel structures have numerous benefits such as high strength, light-weight, fast construction, and good deformability. Steel construction is required to suit the quick development of the infrastructures and the world's economy. In addition, the use of HSS can greatly contribute to the reduction of the cross-sectional size of structural members, reduce the amount of steel and the weight of the structures, leading to the reduction of the work efforts of welding and transportation as well as inducing high resistance to earthquake effects. In general, the use of HSS in construction is good for the environment and saves time[1].

    The advantages to the construction industry of adopting HSS such as Q690 steel have motivated researchers to study the performance of HSS compression members at ambient and elevated temperatures. Shi et al.[2] found that the local stability of an H-shaped member could be ensured by increasing the width of the steel plate. The width-to-thickness ratio of the steel plate in the H-section member should be larger than that of the mild steel, also known as low carbon steel, in order to make full use of the advantages of the strength of high-strength steel. Based on finite element analysis, Knobloch et al.[3] proposed a strain-based calculation method for the local stability of steel members at elevated temperatures. However, this method has not been verified by experimental results, making it risky to use in engineering. Wang et al.[4] suggested a calculation equation of the local buckling stress of Q235 and Q460 axial compression members at elevated temperatures where the limit of the width-to-thickness ratio was provided to restrict the local buckling of the member at elevated temperatures. Additionally, some compression tests were conducted on slender columns of welded sections with steel grades from 460 to 960 N/mm2[5-11]. The results showed that the sections with steel grades higher than 460 exhibited significantly higher buckling resistance than counterparts made of conventional steel types. The improved resistance can be attributed to higher yield strength and lower residual stress to yield strength ratios. It was also found that the end restraints and limited initial imperfections contributed to improvements in resistance. Chiew et al.[12] conducted an elevated temperature test to investigate the strength of reheated, quenched, and tempered (RQT) HSS at elevated temperatures. Post-fire residual strength was evaluated by means of post-fire residual strength tests. Li et al.[13] and Kang et al.[14] investigated the strength deterioration of high strength structural steel Q690 after fire. They found that the post-fire mechanical properties of the structural steel were significantly decreased when subjected to temperatures exceeding 600 ℃. Through these investigations on steel stub columns made of welded sections, it was concluded that the plate slenderness limits for yielding for conventional steel sections can also be applied to the steel sections made of high strength steel. However, sufficient deformation capacities cannot be guaranteed if the limits of plate slenderness for conventional steel compact sections were simply extended to high strength steel sections[15].

    There are many international codes for the design of structural steel at ambient and elevated temperatures; however, the majority could only be applied to the low and intermediate grades of steel. The Chinese standard for design of steel structure (GB 50017-2017)[16] provides the limits of the width-to-thickness ratio of the steel plate for the local buckling of H-section members to ensure that no local buckling occurs before the overall buckling. However, the standard design provisions do not consider high strength steel. The Eurocode 3 (EC3)[17] generally includes a more reliable and safe design formula, which can be applied to the design of the local buckling capacity of high-strength S690 steel I-section members. The American steel structure design code (AISC 360-10)[18] provides slenderness limits to control local buckling for column strength, as expressed in Eq. (1) and (2), where the slenderness limits (λp andλr) of thin-plate elements must be satisfied to prevent local buckling of the plate elements.

    where,fy is the yield strength of steel; Es is the elastic modulus of steel.

    In the EC3[17] provisions, a reduction factor of 0.85 is adopted for the elastic modulus and the yield strength to treat the local buckling phenomena of the steel element subjected to elevated temperatures. Conversely, the Chinese standard for the design of steel structures[19] as well as code for fire safety of steel structures in buildings do not include specifications and guidelines for the design of steel at elevated temperatures. Overall, the current design approach in codes including GB 50017-2017[16], EC3[17], and ANSI/AISC 360-10[18] is only applicable to steel with a yield strength of up to 460 MPa and could not be adopted for the design of HSS columns. Therefore, in practice, the application of high-strength steel members is restricted. However, EC3 Parts 1-12[20] extends its specification to steel grade up to 700 MPa. Furthermore, most steel design codes adopt an equivalent imperfection approach to consider the effects of residual stress and out-of-plane imperfections that is generally assumed to be 1% of the member length. According to the research by Li et al[21], this equivalent value for residual stress would lead to unconservative predictions for high strength steel members. Based on the aforementioned points and on the FE modeling results, a simplified design method and new width-to-thickness ratio limits are proposed for the studied high strength Q690 steel stub columns with H-shaped section.

    In this study, finite element modeling was performed using the software package ABAQUS to simulate the performance of high-strength Q690 steel H-shaped welded columns at elevated temperatures under axial compression. The numerical models were first adopted to replicate the experimental results carried out by Wang et al.[27]. Then, an extensive parametric analysis was conducted for the Q690 HSS welded columns. The effects of many parameters including the width-to-thickness ratio, temperature, initial imperfection, residual stress and interaction between flange and web on the behavior of the studied axially-compressed steel member at elevated temperatures were evaluated. The results of the parametric study were then compared with the results obtained from the design provisions in GB 50017-2017[16], EC3[17], and ANSI/AISC 360-10[18]. Based on the FE results, a new design method was recommended for local buckling design of high strength Q690 steel columns at elevated temperatures.

    2 Material properties of high strength Q690 steel

    The deterioration of the mechanical properties including the elastic modulus and yield strength is considered the main factor influencing the performance of steel structures in fire. As the temperature increases, the yield strength and elastic modulus of the steel are reduced. Wang et al.[22] obtained the mechanical properties of Q690 steel at elevated temperatures through a series of tensile coupon tests and proposed reduction coefficients for the nominal yield strength and elastic modulus of Q690 steel exposed to elevated temperatures. The tests results are summarized in Table 1 and depicted in Fig. 1, where the reduction factors for Q690 at elevated temperatures (20-900 ℃) steel can be obtained from Eq.(3) and (4), expressed as:

    where,fyT and fy are the yield strength of steel at elevated and normal temperatures; ET and E are the elastic modulus of steel at elevated and normal temperatures; Ts is the steel temperature. Fig. 1 depicts the comparison of the test results and the results obtained using the proposed equation of the reduction factors for yield strength and elastic modulus of Q690 steel at elevated temperatures.

    Each type of steel has a particular chemical structure, and high temperatures may affect different material properties. The variations between high strength steel and mild steel in the reduction factors are shown in Fig.2. It can be seen that the stiffness and strength of mild steel are greater than those of high strength steel. This could be related to the variety of the components and formulations of the chemical elements and the differences in the manufacturing techniques. Test data on the reduction factors of yield strength and elastic modulus for mild steel SM41[23], Q235[24], and high strength steel with grades of S460, Q460, S690[25], and Q690[22] were also compared.

    3 Finite element modeling

    This paper utilizes the finite element analysis software ABAQUS 6.14 to simulate the local buckling of high strength Q690 steel welded H-section columns at elevated temperatures under axial compression. The element section type, constitutive relationship, geometric imperfection, residual stress, and boundary conditions are used to build the FE model.

    3.1 Element type, mesh, and material modeling

    The proper selection of element type with appropriate degrees of freedom is essential to simulate the behavior of the stub columns. Various options are available in ABAQUS to model the steel structure. In this paper, the four-noded element quadrilateral S4R shell element with reduced integration, in which each node has six degrees of freedom, is adopted to represent the element type for the steel plate, where lateral shear and large strains are allowed. The number and quality of the finite element mesh directly affect the accuracy of the results and the computation time. Therefore, a suitable mesh is a necessary condition for obtaining accurate simulation results. Based on sensitivity study, the mesh size of the element between 12-20 mm was found suitable. The stress-strain response of the steel at elevated temperatures was simulated using the steel model proposed by Wang et al.[22] for the Q690 steel at elevated temperatures. The stress-strain relationship of the studied Q690 steel can be obtained as follows

    where, σ0.01 andσ0.2 are the static tensile proof stress of 0.01% and 0.2%; m is the exponent for the modified Romberg-Osgood model given in Table 2; n is the strain-hardening exponent; n0 is the original strain-hardening exponent; K is the modular coefficient in the Romberg-Osgood expression; σu is the ultimate stress; εp is the plastic strain; εpu is the plastic strain at the ultimate strength.

    3.2 Boundary conditions and load application

    In order to appropriately apply the constraints and loads for simulating the rigid condition of the ends of the member, rigid endplates were set at both ends of the members. This step was performed as follows: the reference points were first set directly above the endplate and the degrees of freedom of all nodes of the endplate were rigidly connected with the degrees of freedom of the reference point using the constraint coupling command. Hence, rigid ends were generated at both ends of the element and the reference point became the loading point. Then, the contact surface between the steel column and the endplate was defined as the Tie constraint, so that the two surfaces were firmly bonded together to avoid separation during the analysis and achieve similarity with the actual welding. Therefore, it was similar to the actual test and it was possible to avoid local damage to the section of the element under the concentrated load.

    The boundary conditions were fixed at both ends. The nodes at both ends were fixed against the displacement perpendicular to the axial direction and with no twist about the axial direction. The reference point at the top-end was fixed against displacement in all directions except for displacement in the axial direction, while the rotation was fixed in all directions except for a weak bending axis. The reference point at the bottom-end was fixed against the displacement and rotation in all directions. To ensure that the element was a statically indeterminate structure, the longitudinal displacement of one end of the element was restricted. The load was continuously applied according to the displacement control method. Because displacement control method easily converged, the calculation speed was fast, and the descending branch of the load-displacement curve could be obtained. The analysis type of the model was static risk analysis.

    3.3 Residual stresses

    The welding process is an uneven heating and cooling process and this uneven process causes residual stresses in the steel element. These residual stresses could adversely affect the ultimate capacity of the steel element. Therefore, residual stress is an important parameter, which needs to be considered in the FE modeling. Zhang[26] measured the residual stress of H-section Q690 steel and found that the measured distribution model was similar to that of mild steel section.

    Based on the self-balancing principle of cross-section residual stress, the distribution width of residual tensile stress and residual compressive stress can be obtained. The distribution range of the H-section and the ratio of residual stress to the nominal yield strength of steel at ambient temperatures are summarized in Table 3 and 4, respectively. Fig.3 illustrates the distribution range of the residual stress of high strength steel H-section.

    At elevated temperatures, the equation for predicting the residual stress reduction coefficient is reported in Reference [26]. The resulting expression is adopted to calculate the residual stress reduction coefficient for the H-shaped section of high strength Q690 steel as follows

    As shown in Table 5, the reduction coefficient does not reduce significantly when the temperature ranges from 20 to 400 ℃. However, the residual stress reduction coefficient decreases obviously with the increasing of the temperature beyond 400 ℃. Conversely, the reduction coefficient significantly increases when the temperature increases from 500 to 900 ℃. It is also noted that no much reduction is observed for steel with different grades including Q235, Q460, and Q690 at elevated temperatures.

    3.4 Verification of the FE model

    The verification of the current model was carried out using the experimental program of local buckling of high strength Q460 steel stub columns tested under axial load conducted by Wang et al.[27]. A comparison of the FE results and the test results is shown in Fig.4(a), where the load-displacement curves obtained from the FE model and the test data are illustrated. As shown in Fig.4(a), the load-displacement curve obtained from the finite element model is in very good agreement with the load-displacement curve obtained from the test of the specimen. It can also be seen from Fig.4(a) that at the beginning, the load-displacement curves generated by FE model and test showed a slight variation and they continued smoothly during the whole elastic force-displacement stage until they reached the local buckling stress of the test specimen. A slight difference is observed in the descending branch of both the FE generated and the test load-displacement curve because it is very complicated to simulate the exact interaction effect between materials in the FE numerical model and the practical test specimen.

    The failure modes of the specimen Q460-B-20 ℃ and Q460-B-650 ℃ steel columns, and a comparison of the column failure mode observed in the experiment with the failure mode obtained from the FE model is shown in Fig.4(b). It can be seen that the failure model showed by the test fits well with that predicted by the FE model. Due to the position of the initial imperfections in the columns, there are some differences in the location of buckling and the imperfections will result in the first onset of local buckling.

    4 The load bearing capacity of the high strength Q690 steel column

    The design approach for predicting the axial compressive bearing capacity at elevated temperatures is based on the design method of axially compressed steel columns at normal temperatures.

    Reference[19]provide an equation for calculating the critical stress of axially loaded steel members at elevated temperature, expressed as

    where,σE,T is the Euler critical stress of the element at elevated temperature, which can be obtained by (σE,T=π2ET/λ2), and e0 is the initial eccentricity of the elements, for the welded H-shaped, which can be obtained by (e0=0.300λ--0.035) in which, λ- is the slenderness ratio of the column and can be calculated by (λ-=λπ fyT/ET).

    The stability factor of the high strength Q690 steel column can be calculated using the following equation

    Fig.6(a) shows the ultimate stress versus the width-to-thickness ratio of the studied specimens at different temperatures. It can be seen that the ultimate stress can reach the yield strength when the width-thickness ratio decreases, whereas increasing the flange width-to-thickness ratio results in a reduction of the ultimate stress. The ultimate stress was also found to be reduced with the increasing temperature. Fig.6(b) shows the effects of flange width-to-thickness ratio on the buckling stress. The buckling stress decreases with the increasing flange width-to-thickness ratio as well as with the increasing temperature. Similarly, the web ultimate stress reduces with the increasing height-to-thickness ratio, as shown in Fig.7(a). Additionally, the buckling stress of the web reduces with the increasing height-to-thickness ratio, as depicted in Fig.7(b). When the temperature is between 20 ℃ and 550 ℃, the buckling stress decreases rapidly with the increasing plate width(height)-to-thickness ratio. When the temperature is than 550 ℃, the buckling stress decreases slowly with the increase of the width(height)-thickness ratio of the plate. In both cases, the width(height)-to-thickness ratio of the flange or web of the studied specimens significantly influence the ultimate stress and buckling stress of the studied specimens at elevated temperatures.

    5.2 Interaction of web and flange

    When the web or flange of the H-section of the steel member is partially buckled, the buckling stresses occurring in the two elements are not equal. Therefore, the buckled part can be constrained by the unbuckled part. In the engineering design, the H-shaped section of steel consisted of thinner web and thicker flange plates. Consequently, when the flange plate buckled first, the web had a limited effect on the flange plate, hence the restraining effect of the web on the flange can be ignored. On the contrary, when the web plate buckled first, the restraining effect of the flange on the web cannot be ignored. In this investigation, only web buckling is taken, so the restraining effect of the flange plate on the web is considered. Fig.8 depicts the variation of the ultimate stress versus the flange width-to-thickness ratio at elevated temperatures, where the web height-to-thickness ratio is kept at 30 mm. In order to ensure effective constraining effects of the flange on the web when determining the width-to-thickness ratio of the flange, the buckling stress of the web should be greater than that of the flange. As a result, the buckling stress of the web can be attained before the flange is buckled. It can be seen from Fig.8 that when the flange width-to-thickness ratio is 6, the web buckling stress is about 625.1 MPa and with increasing the width-to-thickness ratio of the flange to 20, the web buckling stress is found to decrease to 533.54 MPa, which is a decrease of approximately 14.7%. Similarly, at temperatures of 400~650 ℃, when the flange width-to-thickness ratio increases from 6 to 20, the buckling stress decreases by 17%. Overall, whether the axial compressive steel members are loaded at ambient or elevated temperatures, the retaining effects of the flange on the web should be considered when estimating the web buckling stress.

    5.3 Initial geometrical imperfection

    The initial imperfection of the steel plate is a key factor, which greatly impacts the ultimate bearing capacity of the steel columns. Fig.9 elucidates the local buckling stress of members with different geometrical imperfections. It considers the initial geometrical imperfection of 0.1t (t is the plate thickness). According to the research results of Burgess et al.[29], the thermal expansion of the element at elevated temperatures has little effect on the initial geometric imperfection of the steel element and hence it can be considered that the geometric imperfections at normal temperatures and elevated temperatures are similar. The Chinese code for check & accept specification for steel structure[30] stipulates that the initial geometric imperfection amplitude of the element is 1% of the flange plate gross width. Kaitiala et al.[31] adopted an initial imperfection range similar to that in [30]. The results of buckling stress versus the width-to-thickness ratio of FE models are compared with the calculated ones and it can be observed that both results have good agreement. Hence, the initial imperfections can be used to analyze the buckling stress of the steel specimens.

    6 Simplified design method

    6.1 Local buckling stress at elevated temperatures

    The local buckling stress coefficient φ can be obtained from Eq.(11), while the slenderness ratio (width-to-thickness ratio) of the steel plate can be calculated from Eq.(12), expressed as

    where, fyT is the yield strength of the steel at elevated temperatures, ET is the modulus of elasticity at elevated temperatures, k is the local buckling coefficient related to the boundary conditions. For flange buckling k= 0.425, for web buckling k=4, v is Poisson's ratio of steel v= 0.3.

    Fig.10 illustrates the relationship between the coefficient of local buckling and the slenderness ratio (φλ- curves) at elevated temperatures. It can be seen from Fig. 10 that at elevated temperatures, the local buckling coefficient φ decreases with the increasing slenderness ratio (width-to-thickness ratio) and vice versa. It is also found that when the temperature is between 300 and 700 ℃, the curves are very close. These curves can be used to describe the φλ- relationship.

    Based on the buckling stress results obtained by finite element analysis, the equation for calculating the local buckling stress of the member between 300 ℃ and 700 ℃ was obtained using data fitting analysis. For any H-shaped section, the local buckling stress of the member can be calculated by Eq.(13).

    The Q690 axial compression members at elevated temperatures are

    where, λ- is the regularized width-to-thickness ratio of the plate, which can be calculated according to Eq.(12). In order to verify the reliability of the calculation results of the fitting equation, the buckling stresses calculated by Eqs. (13) to (17) are compared in Fig. 11. For the Q690 axial compression member, the results obtained from Eqs. (15) to (17) are lower than the FE generated results.

    6.2 Width-to-thickness limit

    To ensure that the steel member at elevated temperatures

    under axial compression satisfies the requirements of local buckling, it is important to select appropriate width(height)-to-thickness ratios of both web and flange elements. The relationship between the width(height)-to-thickness ratio and the slenderness ratio should be obtained through equaling the ultimate stress with the local buckling stress and the ultimate stress with the overall buckling stress of the steel element at elevated temperatures. Here, the relationships between the flange width-to-thickness ratio and the web height-to-thickness ratio with the slenderness limits can be expressed, respectively, as

    GB 50017-2017[16] provides the limit values of the web height-thickness ratio and the flange width-to-thickness ratio of axial compression members at ambient temperatures as follows.

    Fig. 12 compares the results of the width-to-thicknessratio and the height-to-thickness ratio with the slenderness limits at normal and elevated temperatures for the flange and web, respectively. In both cases, when the slenderness ratio of the specimen is less than 25, the limit value of the height-to-thickness ratio of the web at elevated temperatures is less than the limit of the height-to-thickness ratio at normal temperatures, while when the slenderness ratio of the component is greater than 25, the plate height-to-thickness ratio limit is larger than that at room temperature. This means that when the length of the element is relatively small, the element at high temperature is more likely to undergo local buckling, whereas, when the length is larger, the element at normal temperature is more likely to undergo local buckling.

    7 Existing international codes of practice

    7.1 Assessment of international codes

    The assessment of the suitability of the current code of practice in designing the studied specimens is conducted through comparison with the FE-generated results, as shown in Fig.13, where the width-to-thickness ratio of the flange and the height-to-thickness ratio of the web are plotted against the slenderness limit in Figs.13(a) and (b), respectively. From Fig.13 (a), it can be seen that GB 50017-2017 provided a more safer estimate for the ultimate stress of the H-shaped section of the web than both Eurocode 3 and ANSI/AISC 360-10 in predicting the ultimate strength, due to their higher design values than the Chinese standard. However, as shown in Fig. 13 (b), when λpf was smaller, the curves obtained from the design codes showed close results with FEM results, especially for the flange connected with quite a slender web. In addition, all design curves obtained from the standards became highly conservative when λpf was relatively large.

    8 Conclusions

    This paper utilized the ABAQUS finite element model to obtain the local buckling behavior of the welded H-section of high strength Q690 steel columns under axial compression. The FE models were first verified against experimental results. An extensive parametric study was conducted using the verified FE model over a wide range of parameters including elevated temperature, width-to-thickness ratio of flange, and height-to-thickness ratio of web, residual stress, and initial imperfections. A comparison with existing codes of practice was conducted to assess the suitability of current codes in predicting the compressive capacities of the studied specimens. A new design method was proposed. Based on these data, the following conclusions were drawn:

    1)At elevated temperatures, the elastic modulus and strength of high strength Q690 deteriorated faster than those of mild steel.

    2)With increases in the flange width-to-thickness ratio and the web height-to-thickness ratio, the ultimate bearing capacity of the member increases, but they have less effect on the ultimate bearing capacity of the members at elevated temperatures. The greater the temperature, the weaker the influence of other factors on the ultimate load-capacity, implying that the temperature of the member plays a leading role in affecting the ultimate bearing capacity of the members.

    3)Whether the axial compressive steel members are loaded at an ambient or elevated temperature, the restraining effect of the flange on the web cannot be ignored and should be considered when estimating the web buckling stress.

    4)Using the initial geometric imperfection of 0.1t (t is the plate thickness) or the initial geometric imperfection of b/100, the buckling stresses obtained by analysis are basically the same.

    5)When the length of the element is relatively small, the element at high temperature is more likely to undergo local buckling, whereas, when the length is larger, the element at normal temperature is more likely to undergo local buckling

    6)A simplified local buckling design method was proposed and recommended for designing the high strength steel section under axial compression. In addition, the current structural steel design codes were evaluated. It was found that GB 50017-2017 tends to provide a more safer estimate for the ultimate stress of the H-shaped section of the web than both EC3 and ANSI/AISC 360-10.

    Acknowledgements

    The authors would like to acknowledge the support from the National Program on Key Research and Development Project (Grant No.2016YFC0701203),and the Natural Science Foundation of China (Grant No.51878096).References:

    [1] ROMEIJN A. Steel-Concrete Bridges-III [M]. Faculty of Civil Engineering and Geosciences, Department Design and Construction-Section Structural and Building Engineering, 2006.

    [2] SHI G, XU K L, LIN C C. Finite element analysis and design method study on the local buckling of 460 megapascal high strength steel I-section columns under axial compression [J]. Industrial Construction, 2016, 46(7): 22-31.

    [3] KNOBLOCH M, FONTANA M. Strain-based approach to local buckling of steel sections subjected to fire [J]. Journal of Constructional Steel Research, 2006, 62(1/2): 44-67.

    [4] WANG W Y, YANG X C, WANG B, et al. Design approach for local stability of axially compressed steel members at elevated temperature [J]. Journal of Chongqing University, 2015, 38(3): 47-57.(in Chinese)

    [5] CHIEW S P, ZHAO M S, LEE C K. Mechanical properties of heat-treated high strength steel under fire/post-fire conditions [J]. Journal of Constructional Steel Research, 2014, 98: 12-19.

    [6] WANG Y B, LI G Q. Experimental and numerical study on the behavior of axially compressed high strength steel box-columns [J]. Engineering Structures, 2014, 58: 79-91.

    [7] BAN H Y, SHI G, SHI Y J, et al. Overall buckling behavior of 460 MPa high strength steel columns: Experimental investigation and design method [J]. Journal of Constructional Steel Research, 2012, 74: 140-150.

    [8] ZHOU F, TONG L W, CHEN Y Y. Experimental and numerical investigations of high strength steel weldedH-section columns [J]. International Journal of Steel Structures, 2013, 13(2): 209-218.

    [9] RASMUSSEN K J R, HANCOCK G J. Tests of high strength steel columns [J]. Journal of Constructional Steel Research, 1995, 34(1): 27-52.

    [10] BAN H Y, SHI G, SHI Y J, et al. Experimental investigation of the overall buckling behaviour of 960 MPa high strength steel columns [J]. Journal of Constructional Steel Research, 2013, 88: 256-266.

    [11] SHI G, BAN H Y, BIJLAARD F S K. Tests and numerical study of ultra-high strength steel columns with end restraints [J]. Journal of Constructional Steel Research, 2012, 70: 236-247.

    [12] WANG Y B, LI G Q, CHEN S W, et al. Experimental and numerical study on the behavior of axially compressed high strength steel columns with H-section [J]. Engineering Structures, 2012, 43: 149-159.

    [13] LI G Q, LYU H B, ZHANG C. Post-fire mechanical properties of high strength Q690 structural steel [J]. Journal of Constructional Steel Research, 2017, 132: 108-116.

    [14] KANG L, SUZUKI M, GE H B, et al. Experiment of ductile fracture performances of HSSS Q690 after a fire [J]. Journal of Constructional Steel Research, 2018, 146: 109-121.

    [15] YUAN B. Local buckling of high strength steel W-shaped sections [D]. McMaster University, 2018.

    [16] Code for Design of Steel Structures: GB 50017-2017 [S]. Beijing: China Architecture & Building Press, China Plan Press, 2017.

    [17] Eurocode 3:Design of steel structures. Part 1.2: General rules-Structural fire design: EN 1993-1-2 [S]. Brussels: European Committee for Standardization (ECS), 2005.

    [18] Specifications for Structural Steel Buildings: ANSI/AISC 360-10 [S]. Chicago, 2010.

    [19] Technical code for fire safety of steel structure in buildings: CECS 200:2006 [S]. Beijing: China Plan Press, 2006.

    [20] Eurocode 3: Design of Steel Structures-Part 1-12: Additional Rules for the Extension of EN 1993 up to Steel Grades S700: European 1993-1-12 [S]. Brussels, 2007.

    [21] LI T J, LI G Q, WANG Y B. Residual stress tests of welded Q690 high-strength steel box- and H-sections [J]. Journal of Constructional Steel Research, 2015, 115: 283-289.

    [22] WANG W Y, WANG K, KODUR V, et al. Mechanical properties of high-strength Q690 steel at elevated temperature [J]. Journal of Materials in Civil Engineering, 2018, 30(5): 04018062.

    [23] LI G Q, ZHANG X J. Experimental studies of the material properties of SM41 steel at elevated temperatures [J]. Industrial Construction, 2001, 31: 57-59.

    [24] TAN W. Experiments and research of steel material properties at elevated temperature [J]. Industrial Construction, 2010, 30: 61-63, 67.

    [25] WANG W Y, KODUR V. Introduction[M]//Material Properties of Steel in Fire Conditions. Amsterdam: Elsevier, 2020: 1-27.

    [26] ZHANG J. Study on residual stress and load capacity of welded high-strength Q690 steel column after fire exposure [D]. Chongqing: Chongqing University, 2019.

    [27] WANG W Y, KODUR V, YANG X C, et al. Experimental study on local buckling of axially compressed steel stub columns at elevated temperatures [J]. Thin-Walled Structures, 2014, 82: 33-45.

    [28] LI X, WANG W Y. Local bucking of Q690 high strength steel columns subjected to axial compression at elevated temperature [J]. Progress in steel building structures, 2020.

    [29] BURGESS I W, OLAWALE A O, PLANK R J. Failure of steel columns in fire [J]. Fire Safety Journal, 1992, 18(2): 183-201.

    [30] Check & accept specification for steel structure work: GB 50205-2001 [S]. Beijing: China Plan Press, 2003.

    [31] KAITILA O. Imperfection sensitivity analysis of lipped channel columns at high temperatures [J]. Journal of Constructional Steel Research, 2002, 58(3): 333-351.

    (編輯 章潤紅)

    猜你喜歡
    鋼柱有限元分析高溫
    全球高溫
    望火興嘆
    跨多層支護(hù)結(jié)構(gòu)的地下室鋼柱施工技術(shù)
    跨多層支護(hù)結(jié)構(gòu)的地下室鋼柱施工技術(shù)
    “烤”問城市
    自錨式懸索橋鋼箱梁頂推施工階段結(jié)構(gòu)分析
    探究跨多層支護(hù)結(jié)構(gòu)的地下室鋼柱施工技術(shù)
    隨機(jī)振動(dòng)載荷下發(fā)射裝置尾罩疲勞壽命分析
    有限元分析帶溝槽平封頭的應(yīng)力集中
    飛機(jī)起落架支撐桿強(qiáng)度有限元分析
    亚洲一区中文字幕在线| 黑人巨大精品欧美一区二区mp4| 一区二区三区高清视频在线| 国产又爽黄色视频| 黑丝袜美女国产一区| 久久香蕉激情| 中文字幕av电影在线播放| 成年版毛片免费区| 美女免费视频网站| 丝袜美腿诱惑在线| 两人在一起打扑克的视频| 国产精品野战在线观看| 91字幕亚洲| 91精品三级在线观看| 又大又爽又粗| 午夜久久久久精精品| 激情在线观看视频在线高清| 久久精品aⅴ一区二区三区四区| 日本免费a在线| 极品人妻少妇av视频| 后天国语完整版免费观看| 久久伊人香网站| 欧美乱码精品一区二区三区| 老熟妇仑乱视频hdxx| 精品乱码久久久久久99久播| 丝袜人妻中文字幕| 久久影院123| 欧美精品亚洲一区二区| 99re在线观看精品视频| 电影成人av| 色av中文字幕| 国产精品乱码一区二三区的特点 | 午夜老司机福利片| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 国产亚洲欧美在线一区二区| 亚洲人成伊人成综合网2020| 欧美乱色亚洲激情| 久久久久亚洲av毛片大全| 亚洲av日韩精品久久久久久密| 午夜精品久久久久久毛片777| 久久人人爽av亚洲精品天堂| 黄色成人免费大全| 午夜福利,免费看| 嫁个100分男人电影在线观看| 日本三级黄在线观看| 久久久久国产精品人妻aⅴ院| 国产成人系列免费观看| 亚洲av第一区精品v没综合| 首页视频小说图片口味搜索| 欧美日韩精品网址| 国产精品二区激情视频| 69精品国产乱码久久久| 高潮久久久久久久久久久不卡| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩一区二区三| 中国美女看黄片| 国产精品久久电影中文字幕| 免费观看精品视频网站| 日韩免费av在线播放| 免费高清视频大片| 真人一进一出gif抽搐免费| 人人澡人人妻人| 国产一级毛片七仙女欲春2 | 99精品久久久久人妻精品| 一级毛片女人18水好多| av天堂在线播放| 国产av精品麻豆| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品av在线| 亚洲av成人不卡在线观看播放网| 51午夜福利影视在线观看| 成熟少妇高潮喷水视频| 久久精品91蜜桃| 亚洲第一青青草原| 国产麻豆成人av免费视频| 丝袜美足系列| 免费看a级黄色片| 黄色视频,在线免费观看| 国产麻豆69| 亚洲午夜理论影院| 国产野战对白在线观看| 久久亚洲真实| 怎么达到女性高潮| 国产亚洲精品第一综合不卡| av网站免费在线观看视频| 日韩欧美在线二视频| 亚洲色图av天堂| 国产av一区二区精品久久| 18美女黄网站色大片免费观看| 黄色丝袜av网址大全| 中文亚洲av片在线观看爽| 午夜福利,免费看| 欧美色视频一区免费| 啪啪无遮挡十八禁网站| 黄色a级毛片大全视频| 每晚都被弄得嗷嗷叫到高潮| 香蕉丝袜av| 精品午夜福利视频在线观看一区| 中文亚洲av片在线观看爽| 亚洲熟妇中文字幕五十中出| 女人精品久久久久毛片| 在线观看66精品国产| 欧美老熟妇乱子伦牲交| 午夜老司机福利片| 97人妻天天添夜夜摸| 亚洲视频免费观看视频| 午夜老司机福利片| 欧美精品亚洲一区二区| 亚洲中文av在线| 久久这里只有精品19| 欧美乱码精品一区二区三区| 欧美午夜高清在线| 真人一进一出gif抽搐免费| 一级毛片高清免费大全| 亚洲国产精品sss在线观看| av视频在线观看入口| 午夜成年电影在线免费观看| 少妇粗大呻吟视频| 两个人视频免费观看高清| 久久伊人香网站| 日本黄色视频三级网站网址| 俄罗斯特黄特色一大片| www.熟女人妻精品国产| 老司机靠b影院| 亚洲av日韩精品久久久久久密| 成人国语在线视频| 亚洲三区欧美一区| 19禁男女啪啪无遮挡网站| 日本vs欧美在线观看视频| 日韩国内少妇激情av| 亚洲片人在线观看| av福利片在线| 午夜激情av网站| 国产精品免费一区二区三区在线| 日本vs欧美在线观看视频| 日韩精品中文字幕看吧| 亚洲精品美女久久av网站| 精品无人区乱码1区二区| 亚洲专区字幕在线| 亚洲av美国av| 性欧美人与动物交配| 777久久人妻少妇嫩草av网站| 久久中文看片网| 精品熟女少妇八av免费久了| 精品福利观看| 淫秽高清视频在线观看| 欧美一级毛片孕妇| 老司机靠b影院| 免费看美女性在线毛片视频| 亚洲国产中文字幕在线视频| 999久久久国产精品视频| 欧美日韩亚洲综合一区二区三区_| www.精华液| 精品国产超薄肉色丝袜足j| 久久久久久国产a免费观看| 无人区码免费观看不卡| 九色亚洲精品在线播放| 国产精品爽爽va在线观看网站 | 久久国产精品人妻蜜桃| 亚洲少妇的诱惑av| 禁无遮挡网站| 日韩国内少妇激情av| 18禁国产床啪视频网站| 很黄的视频免费| www日本在线高清视频| 国产私拍福利视频在线观看| 美女 人体艺术 gogo| 天天一区二区日本电影三级 | 精品高清国产在线一区| 日本欧美视频一区| 午夜福利视频1000在线观看 | 日韩av在线大香蕉| 一级毛片高清免费大全| 国产欧美日韩一区二区三| 热99re8久久精品国产| 欧美国产精品va在线观看不卡| 亚洲激情在线av| 女警被强在线播放| 女同久久另类99精品国产91| 久久中文字幕人妻熟女| 久久久久久久久免费视频了| 男人舔女人下体高潮全视频| 午夜老司机福利片| 精品国产亚洲在线| 久久精品国产亚洲av高清一级| 久久人妻av系列| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看 | 亚洲黑人精品在线| 亚洲 欧美一区二区三区| 美女大奶头视频| 在线观看免费日韩欧美大片| 在线十欧美十亚洲十日本专区| 亚洲成人国产一区在线观看| 可以免费在线观看a视频的电影网站| 一本大道久久a久久精品| 久久婷婷成人综合色麻豆| 他把我摸到了高潮在线观看| 国产成年人精品一区二区| 一级,二级,三级黄色视频| 精品无人区乱码1区二区| 91精品国产国语对白视频| 丁香欧美五月| 人人妻人人澡人人看| 日韩精品中文字幕看吧| 嫩草影院精品99| 日韩欧美在线二视频| 午夜福利高清视频| www.www免费av| 韩国av一区二区三区四区| 亚洲成人国产一区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 女人被躁到高潮嗷嗷叫费观| 久久天躁狠狠躁夜夜2o2o| 国产成人精品久久二区二区免费| 国产私拍福利视频在线观看| 校园春色视频在线观看| 精品少妇一区二区三区视频日本电影| 国产精品一区二区三区四区久久 | 午夜福利在线观看吧| 一级毛片女人18水好多| 亚洲视频免费观看视频| 亚洲免费av在线视频| 久久香蕉精品热| 欧美激情高清一区二区三区| 国产成人免费无遮挡视频| 一区福利在线观看| 女人精品久久久久毛片| 麻豆久久精品国产亚洲av| 日本a在线网址| 精品国内亚洲2022精品成人| 一进一出好大好爽视频| 精品午夜福利视频在线观看一区| 欧美一区二区精品小视频在线| 久久久久久久久久久久大奶| 老司机深夜福利视频在线观看| 长腿黑丝高跟| 国产精品久久电影中文字幕| 在线视频色国产色| 国产成人av教育| av免费在线观看网站| 成人亚洲精品一区在线观看| 精品国产美女av久久久久小说| 亚洲欧美激情综合另类| 老司机靠b影院| 啪啪无遮挡十八禁网站| 欧美 亚洲 国产 日韩一| www.999成人在线观看| 午夜久久久在线观看| 亚洲av成人不卡在线观看播放网| 69精品国产乱码久久久| 夜夜看夜夜爽夜夜摸| 性色av乱码一区二区三区2| 窝窝影院91人妻| 国产精品久久久人人做人人爽| 女人精品久久久久毛片| 日韩欧美国产一区二区入口| 精品熟女少妇八av免费久了| 身体一侧抽搐| 禁无遮挡网站| 麻豆成人av在线观看| 极品教师在线免费播放| 精品高清国产在线一区| 国产午夜精品久久久久久| 国产av一区在线观看免费| 在线视频色国产色| 夜夜夜夜夜久久久久| 亚洲成人久久性| 久久久久久久久久久久大奶| 色老头精品视频在线观看| 久久精品人人爽人人爽视色| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 日韩欧美一区视频在线观看| 欧美另类亚洲清纯唯美| 久久国产乱子伦精品免费另类| 国产欧美日韩一区二区精品| 国产午夜福利久久久久久| 日韩av在线大香蕉| 757午夜福利合集在线观看| 国产精品爽爽va在线观看网站 | 在线观看舔阴道视频| 十八禁网站免费在线| 欧洲精品卡2卡3卡4卡5卡区| 男女下面进入的视频免费午夜 | 成人18禁在线播放| 中文字幕人成人乱码亚洲影| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩综合在线一区二区| 少妇粗大呻吟视频| 久9热在线精品视频| 国产亚洲精品久久久久久毛片| 狠狠狠狠99中文字幕| 午夜日韩欧美国产| 午夜福利免费观看在线| 男男h啪啪无遮挡| 深夜精品福利| 精品第一国产精品| 一进一出好大好爽视频| 丰满人妻熟妇乱又伦精品不卡| 欧美av亚洲av综合av国产av| 精品欧美国产一区二区三| 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 国产精品免费视频内射| 叶爱在线成人免费视频播放| 国产亚洲精品综合一区在线观看 | 一进一出好大好爽视频| svipshipincom国产片| 美女大奶头视频| 一区二区三区精品91| 精品久久久久久久人妻蜜臀av | 亚洲精品久久成人aⅴ小说| 亚洲aⅴ乱码一区二区在线播放 | 人妻丰满熟妇av一区二区三区| 国产av一区二区精品久久| 成人三级做爰电影| 日韩 欧美 亚洲 中文字幕| 欧美黄色淫秽网站| 色在线成人网| 久久久久久亚洲精品国产蜜桃av| 亚洲伊人色综图| 99精品久久久久人妻精品| 久久久久久免费高清国产稀缺| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 欧美成人性av电影在线观看| 色婷婷久久久亚洲欧美| 日日干狠狠操夜夜爽| 搡老熟女国产l中国老女人| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 日韩欧美三级三区| 丰满人妻熟妇乱又伦精品不卡| 国产精品香港三级国产av潘金莲| 精品久久久久久成人av| 久久精品aⅴ一区二区三区四区| 欧美一级毛片孕妇| 男人舔女人下体高潮全视频| 午夜久久久在线观看| 国产不卡一卡二| 中出人妻视频一区二区| 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| av天堂在线播放| 久久亚洲真实| 高潮久久久久久久久久久不卡| 亚洲人成网站在线播放欧美日韩| 国产不卡一卡二| 国产精品久久视频播放| 老司机在亚洲福利影院| av超薄肉色丝袜交足视频| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 亚洲av第一区精品v没综合| 老司机深夜福利视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 超碰成人久久| а√天堂www在线а√下载| 又黄又粗又硬又大视频| 亚洲精品中文字幕一二三四区| 中文字幕精品免费在线观看视频| av视频在线观看入口| 久久久久国产精品人妻aⅴ院| 欧美成人性av电影在线观看| 午夜激情av网站| 一级片免费观看大全| 亚洲av五月六月丁香网| 又大又爽又粗| 老汉色∧v一级毛片| 久久九九热精品免费| 成人特级黄色片久久久久久久| 一级毛片女人18水好多| 国产一区二区在线av高清观看| 日日夜夜操网爽| 国产一区二区在线av高清观看| 国产精品野战在线观看| 国产av精品麻豆| av电影中文网址| 久久人妻福利社区极品人妻图片| 久久精品亚洲精品国产色婷小说| 日本a在线网址| 国产精品国产高清国产av| 色播亚洲综合网| 国产精品98久久久久久宅男小说| a在线观看视频网站| 日本黄色视频三级网站网址| 99精品久久久久人妻精品| 日韩免费av在线播放| 国产成人av教育| 亚洲少妇的诱惑av| 精品一区二区三区四区五区乱码| 日韩成人在线观看一区二区三区| 亚洲五月天丁香| 久久久久久国产a免费观看| 91精品国产国语对白视频| 91成年电影在线观看| 亚洲在线自拍视频| 美女高潮到喷水免费观看| 国产成人精品久久二区二区免费| 日本免费a在线| 女性被躁到高潮视频| 免费人成视频x8x8入口观看| 一二三四在线观看免费中文在| 一个人观看的视频www高清免费观看 | √禁漫天堂资源中文www| 丰满的人妻完整版| 欧美大码av| 亚洲熟妇中文字幕五十中出| 色播在线永久视频| 最近最新中文字幕大全免费视频| 一级a爱视频在线免费观看| 中国美女看黄片| 日本黄色视频三级网站网址| 女警被强在线播放| av天堂在线播放| 国产97色在线日韩免费| 日本 av在线| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 99精品在免费线老司机午夜| 亚洲九九香蕉| 国产精品98久久久久久宅男小说| tocl精华| 国产男靠女视频免费网站| 久久久国产成人精品二区| 巨乳人妻的诱惑在线观看| 国产精品久久电影中文字幕| 亚洲第一青青草原| 波多野结衣巨乳人妻| 日韩欧美免费精品| 国产精品98久久久久久宅男小说| 后天国语完整版免费观看| www.自偷自拍.com| 精品午夜福利视频在线观看一区| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 香蕉丝袜av| 精品福利观看| 制服诱惑二区| 人人妻人人澡欧美一区二区 | 亚洲男人天堂网一区| 99精品在免费线老司机午夜| 人人妻人人澡欧美一区二区 | 18禁国产床啪视频网站| 成人欧美大片| 一夜夜www| 欧美精品亚洲一区二区| 国产熟女午夜一区二区三区| 最好的美女福利视频网| 国产亚洲欧美精品永久| 国产精品二区激情视频| 国产熟女xx| 国产成人精品无人区| 亚洲欧美激情在线| 精品国产一区二区久久| 一级a爱片免费观看的视频| 日本在线视频免费播放| 99精品欧美一区二区三区四区| 国产激情欧美一区二区| 叶爱在线成人免费视频播放| 久久国产精品影院| 久久香蕉激情| 搡老岳熟女国产| 国产高清videossex| 欧美日本中文国产一区发布| 熟妇人妻久久中文字幕3abv| 1024视频免费在线观看| 亚洲全国av大片| 欧美中文日本在线观看视频| 国产欧美日韩一区二区三| 久久午夜综合久久蜜桃| 国产免费av片在线观看野外av| 亚洲男人的天堂狠狠| 亚洲国产看品久久| 欧美中文综合在线视频| 国产精品香港三级国产av潘金莲| 欧美日韩亚洲国产一区二区在线观看| 免费少妇av软件| 校园春色视频在线观看| 亚洲国产看品久久| 91麻豆av在线| 国产野战对白在线观看| 在线av久久热| 午夜影院日韩av| 十分钟在线观看高清视频www| 亚洲伊人色综图| 母亲3免费完整高清在线观看| 日本一区二区免费在线视频| 欧美亚洲日本最大视频资源| 久久九九热精品免费| 成人特级黄色片久久久久久久| 久久久久九九精品影院| 黄网站色视频无遮挡免费观看| 丝袜人妻中文字幕| 精品国产亚洲在线| 国产av又大| 亚洲欧美日韩另类电影网站| av视频在线观看入口| 在线观看舔阴道视频| 亚洲九九香蕉| 搡老熟女国产l中国老女人| 天天添夜夜摸| 午夜福利影视在线免费观看| 久久久精品欧美日韩精品| 精品高清国产在线一区| 视频区欧美日本亚洲| 怎么达到女性高潮| 女生性感内裤真人,穿戴方法视频| 看片在线看免费视频| 久久久久久国产a免费观看| 一级毛片高清免费大全| 十八禁网站免费在线| 一级毛片女人18水好多| 欧美日本中文国产一区发布| 久久国产亚洲av麻豆专区| 啦啦啦 在线观看视频| 波多野结衣巨乳人妻| 999久久久国产精品视频| 亚洲av电影在线进入| 看黄色毛片网站| 午夜影院日韩av| 欧美黄色淫秽网站| 老熟妇仑乱视频hdxx| 可以在线观看的亚洲视频| 亚洲中文字幕一区二区三区有码在线看 | 99re在线观看精品视频| 亚洲欧美激情综合另类| 午夜精品国产一区二区电影| 免费不卡黄色视频| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 亚洲电影在线观看av| 亚洲五月天丁香| 搡老熟女国产l中国老女人| 美女高潮到喷水免费观看| 国产精品一区二区在线不卡| 在线永久观看黄色视频| 欧美大码av| 女生性感内裤真人,穿戴方法视频| 国产激情欧美一区二区| 亚洲狠狠婷婷综合久久图片| 国产熟女xx| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费| 在线观看免费午夜福利视频| 99国产极品粉嫩在线观看| 在线av久久热| 欧美在线黄色| 欧美国产日韩亚洲一区| 18禁国产床啪视频网站| 脱女人内裤的视频| 黄色视频,在线免费观看| www日本在线高清视频| 中亚洲国语对白在线视频| 麻豆成人av在线观看| tocl精华| 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲国产毛片av蜜桃av| 国产成人av教育| 国产高清视频在线播放一区| 久久国产乱子伦精品免费另类| 一本久久中文字幕| 天天躁夜夜躁狠狠躁躁| 很黄的视频免费| 亚洲欧美日韩无卡精品| x7x7x7水蜜桃| 久久国产精品男人的天堂亚洲| 最近最新中文字幕大全电影3 | av欧美777| 最近最新中文字幕大全电影3 | 欧美激情 高清一区二区三区| 国产99久久九九免费精品| 在线天堂中文资源库| 黑人操中国人逼视频| 成在线人永久免费视频| 国产亚洲精品综合一区在线观看 | xxx96com| 99热只有精品国产| 午夜免费鲁丝| 丝袜美腿诱惑在线| 女性被躁到高潮视频| 亚洲人成电影观看| 久久人人97超碰香蕉20202| 亚洲第一欧美日韩一区二区三区| 日韩精品免费视频一区二区三区| 久久人妻av系列| 天堂动漫精品| 亚洲av片天天在线观看| 露出奶头的视频| 国产1区2区3区精品| 黄片小视频在线播放| 国产激情久久老熟女| 波多野结衣av一区二区av| 国产精品美女特级片免费视频播放器 | 香蕉丝袜av| 看黄色毛片网站| 在线天堂中文资源库| 国产精品一区二区免费欧美| 国产av一区二区精品久久| 成人三级做爰电影| 日韩成人在线观看一区二区三区| 久久这里只有精品19| 老熟妇仑乱视频hdxx| 精品欧美国产一区二区三| 在线观看免费视频日本深夜| 亚洲av熟女| 欧美 亚洲 国产 日韩一| 免费高清在线观看日韩|