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

    A novel modification method for the dynamic mechanical test using thermomechanical analyzer for composite multi-layered energetic materials

    2023-03-28 08:37:34LeQiShilinZhngHoYunZhonglingZhonglingXio
    Defence Technology 2023年3期

    Le Qi , Shi-lin Zhng , Ho Yun , Zhong-ling M ,*, Zhong-ling Xio

    a School of Environment and Safety Engineering, The North University of China, Taiyuan, 030051, China

    b School of Chemical Engineering, The Nanjing University of Science and Technology, Nanjing, 210094, China

    Keywords:Thermo-mechanical analyzer (TMA)Dynamic mechanical analyzer (DMA)Modified Time-temperature equivalence principle

    ABSTRACT Herein, we present a thermo-mechanical analyzer (TMA) and dynamic mechanical analyzer (DMA) of composite multi-layered gun propellant, focusing on thermal expansion coefficients and dynamic thermomechanical properties. The linear thermal expansion coefficient of the prepared energetic material is determined as approx. 0.1800 × 10-4 - 0.2081 × 10-4 K-1. According to DMA test and dynamic thermomechanical properties, the glass transition temperature is also obtained. The tested value is within the range of 223.01-223.50 K,which indicates the lower limit of the energetic material.However,DMA tests reveal temperature changes, which occur due to thermal expansion. Moreover, the geometrical factor decreases with increasing temperature.Therefore,thermal expansion significantly affects the storage modulus and loss modulus.Additionally,the thermal expansion coefficient can be used to modify the storage and loss modulus.The results show that the proposed method provides effective and reliable modified results.

    1. Introduction

    Double-based (DB) gun propellant is produced from mixing nitrocellulose and nitroglycerine with small amounts of other auxiliary materials [1,2]. DB propellants are widely used for medium and small caliber weapons;hence,many studies have focused on improving their combustion performance [3,4]. The composite multi-layered flake gun propellant has good burning progressivity,which can effectively improve the muzzle velocities of the gun's projectile.Additionally,the mechanical properties of the composite multi-layered flake gun propellant are similar to that of general polymers [5,6]. Therefore, a testing method for general polymers can be used for this composite multi-layered flake gun propellant.

    Thermo-mechanical analyzer(TMA)is commonly used to study the thermal expansion of energetic materials. The thermal expansion coefficient of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine(HMX) has been measured over the temperature range of 58-253 K using TMA. Furthermore, the coefficient of thermal expansion (CTE) value is affected by the temperature range and heating rate, and is within the range of 51-66 × 10-6K-1[7].Thermo-mechanical analyzer (TMA) technique has been used to characterize the thermomechanical properties of carbon nanotubes and graphene nano-platelets[8].This approach has proven to be a simple and effective method to measure the glass transition temperature [9]. Victor et al. investigated the nitramine-based amorphous energetics using TMA, and found that the material's softening was associated with the glass transition at approx.323 K[10]. The relaxation transitions, surface tension, and porosity of nitrocellulose have been studied using TMA and dynamic mechanical analyzer (DMA), where three small relaxation transitions of nitrocellulose can be observed, which are contained in gun propellant at 233 K, 308 K, and 353 K [11,12]. CTE of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) based explosives has been measured using TMA within the temperature range of 300-566 K,in which the axial and radical CTE are 2.7 × 10-6K-1and 48 × 10-6K-1, respectively, which indicate a small degree of alignment of TATB molecules under uniaxial compression [13,14].

    The thermo-mechanical properties tests can be divided into TMA and DMA tests. Furthermore, the thermo-mechanical properties of DB gun propellants influence their proper functioning,which has been the critical focus of many studies. In our previous works, linear thermal expansion coefficient of composite multilayered flake gun propellant was measured using TMA method[3,15]. DMA method can also be used to study the thermal expansion coefficient. However, this method is complicated [16]. It is known that the glass transition temperature is related to phase change and can be evaluated using several techniques, which are attributed to the frequency effect.The glass transition temperature obtained using DMA and temperature modulated differential scanning calorimetry (TMDSC) are different [17-19]. The glass transition temperature of automotive coatings has been measured using TMA and DMA. Reports have shown that the correlation between the temperature of midpoint in DSC, the onset temperature of TMA,and the onset temperature of loss modulus using DMA at low frequencies is good. However, this correlation decreases at high frequencies [20,21]. Most glass transition temperatures of energetic materials have been studied using DMA method[22-24],including perfluorosulfonic acid ionomers.The graph between tanδ and temperature shows four distinct peaks. The glass transition temperature of H+from Nafion corresponds to a weak β-relaxation that is centered at 253 K[25].Furthermore,DMA method is widely used in the research of double-based solid propellant to study the glass transition temperature [26,27]. The dynamic mechanical properties of double-based rocket propellant have been studied at 353 K, 358 K, and 363 K, where the storage and loss modulus depend on temperature and time. However, the prediction of the propellant's artificially aged time has been determined to be unreliable above 333 K [28]. The viscoelastic parameters obtained from DMA measurement can be optimized using a genetic algorithm, which are matched with the experimental results at different temperatures, strain rates, and pressures [29]. Additionally,DMA technique can be used to predict the artificially aged time of the propellant by storage modulus [30].

    The purpose of this work is to provide a fundamental background on the method of modifying the dynamic mechanical test results using thermo-mechanical analyzer for energetic materials,which provides a modified reference for DMA method. TMA and DMA tests are used to describe the thermo-mechanical properties of the energetic material. Moreover, frequency sweep curves are obtained using DMA. According to the thermal expansion coefficient, modification of the storage and loss modulus is an effective and reliable method.

    2. Theoretical backgr ound

    In DMA the dynamic mechanical behavior of a material acts as a function of frequency and temperature. The sample mounted in DMA is subjected to a predefined mechanical oscillation program,defined by the frequency and amplitude,and temperature program.The dynamic displacement of the sample during the mechanical oscillation program is measured.The phase lag angle δ between the applied stimulus and the response is also obtained using the displacement measurement.Additionally,a force signal is produced based on the force exerted on the sample by the drive motor.

    The force and displacement are denoted as F and L,respectively.The complex modulus E*can be calculated using Eq. (1).

    where the ratio F/L is defined as the stiffness S, and g is known as the geometrical factor, which is calculated from the sample's dimensions. Both S and g are calculated using Eq. (2) and Eq. (3),respectively.

    where g is the geometrical factor for three-point bending specimen,and l, w, and b are the geometric dimensions of the sample identified in Fig.1.A schematic diagram of DMA test is shown in Fig.1.

    Fig.1. Schematic of functional principle of DMA.

    During DMA method, as temperature increases, the thermal expansion of specimens affects their geometry. The probe's pressure on specimens is 0.05 N in TMA method,whereas in the case of DMA method there is 1 N pressure on the upper clamp bar to ensure the assumed amplitude value can be 10 μm.The geometrical factor(g) can be corrected by the thermal expansion coefficient of the specimen during the heating process at the same heating rate and temperature range.

    Therefore, a method is required that is suitable for composite multi-layered flake gun propellant, which is easy to apply and calculate.The purpose of this paper is to describe such a method to correct for the effect of temperature on geometrical factor.Thermal expansion of the specimen in the thickness direction is considered,and the linear thermal expansion coefficient of the specimen is defined as α. When combined with the definition of thermal expansion coefficient,the thickness of the specimen at temperature T is expressed as Eq. (4).

    where b0(T0) is the original thickness of the specimen, and b(T) is the thickness of the specimen at temperature T.Eq.(3)can also be rewritten as Eq. (5).

    The complex modulus E*is expressed as Eq. (6).

    Therefore, the modulus can be corrected by thermal expansion coefficient.

    3. Experiment

    3.1. Materials

    All samples described herein were composite multi-layered flake gun propellant. The gun propellant had a layered structure.Both the inner and outer layers were composed of double-base propellant. The main component of the inner layered material was 1# double-base absorbent propellant (NC:68.2%, NG:31%,C2:0.8%, where NC is nitrocellulose, NG is nitroglycerin, C2 is centralite, the nitration level is 13.00%). The main component of the outer layered material was 2# double-base absorbent propellant(NC:88.26%, NG:10.04, C2:1.7%, the nitration level is 12.76%). The structure of the composite multi-layered flake gun propellant is shown in Fig.2,inner layer was coated by outer layer,and the inner layer was twice as thick as outer layer, the corresponding geometrical size of the samples is shown in Table 1. And the solubility parameters of nitrocellulose and alcohol/acetone mixed solvents (1:1) are similar, therefore, alcohol/acetone mixed solvents are used as solvent.

    Table 1 Geometrical size of the samples composite multi-layered flake gun propellant.

    Fig. 2. The structure of composite multi-layered flake gun propellant.

    3.2. Thermo-mechanical analyzer (TMA)

    The theory of CTE comprises the deformation of materials was determined as the temperature changes. In some temperature ranges, the linear CTE was constant. When the temperature changed, TMA measured the change in length of materials. The linear CTE of materials was calculated by measuring its deformation within a specific temperature range. The instrument's sensitivity was 0.1 μm.

    Furthermore, the composite multi-layered flake gun propellant belongs to double-based gun propellant, which turns soft above 300 K. In order to reduce the influence of probe on the thickness measurement, a flat probe 3 mm in diameter was used to reduce the pressure of probe on the specimens. TMA results showed that 0.007 MPa causes false negative expansion of the sample above 300 K;however,in DMA test,the specimens were subjected to the same pressure. ΔL values were obtained at various temperature from ΔL-T curves; hence, the thickness of the specimens were modified by ΔL-T curves in DMA test.

    A smooth part of the composite multi-layered flake gun propellant sample (double-base propellant) was selected. The sample was a flake shaped cuboid with the length,width and thickness of 0.5 cm, 0.5 cm, and 0.3-1.0 mm, respectively. The thermal expansion coefficient of the sample was measured using TMA-Q 400(TA instruments, USA) with 2 K/min heating rate under dry nitrogen atmosphere over the temperature range of 213-323 K. The linear CTE can be expressed using Eq. (7).

    where α, Δb, b0, and ΔT are the linear CTE, change in dimension,original thickness of the material, and corresponding temperature change, respectively. The linear CTE (α) is average linear thermal expansion coefficient for the sample in the temperature range of T2to T1.

    3.2.1. Dynamic mechanical analyzer (DMA)

    DMA test was performed using DMA1 analyzer(Mettler Toledo)in 3-point bending mode with the sample dimensions of 60 ± 0.05× 7 ± 0.05 ×b0(0.3-1.0) mm.A liquid nitrogen cooling accessory was used for the experiments. The non-isothermal measurements were conducted within the temperature range of 273-393 K at the heating rate of 2 K/min and loading frequency of 1 Hz. The amplitude of the test sample's deformation was established at ±10 μm. The isothermal measurements were carried out over the frequency range of 100-0.01 Hz with 10 data points for each order of magnitude. Moreover, the amplitude was 10 μm for the temperatures of 228 K, 263 K, and 293 K.

    4. Experimental results and analysis

    4.1. TMA test results

    TMA linear thermal expansion curve of specimens with thickness between 0.3 and 1.0 mm, scanned within the temperature range of 263-283 K,is shown in Fig.3.The obtained results showed that the curves were all linear.Between this temperature range,the linear thermal expansion coefficient lied within the range of 0.1800 × 10-4- 0.2081 × 10-4K-1, whereas the linear CTE increased with increasing sample thickness (see Table 2). The composite multi-layered flake gun propellant expansion during this stage was mainly caused by thermal expansion of nitrocellulose. The specimens were extruded under 5 MPa pressure in the forming process.However,solvent still remained in the specimens,as well as having creep recovery characteristics. In the forming process, as the sample was extruded by smaller export molds, the thickness direction of the sample experienced high pressure. The nitrocellulose was also subject to greater bondage stress. As the stress disappeared,nitrocellulose molecule rebounded greatly,and the solvent in the samples was completely volatilized after a long period of time. Hence, nitrocellulose molecule recovery occurred.During TMA test, as the temperature increased, the thermal expansion rebounding of nitrocellulose molecule in samples of smaller thicknesses was small. Therefore, the linear CTE increases with increasing thickness.

    Table 2 Linear thermal expansion coefficient of composite multi-layered flake gun propellant.

    Table 3 Transition temperature of composite multi-layered flake gun propellant.

    The ΔL-T curves of composite multi-layered flake gun propellant with different thickness values within the temperature range of 233-323 K is shown in Fig.4.There was an obviously turning point on each curve at about 300 K,which indicated that the sample had softened above this transition temperature.Therefore,the thermal expansion coefficient of composite multi-layered flake gun propellant could not be calculated by Eq. (4) at temperatures beyond 300 K. This was due to the sample being compressed because of probe pressure beyond 300 K. As a result, a false negative expansion trend above 300 K was observed. However, ΔL-T curves showed shape variable values of thickness for composite multi-layered flake gun propellant with temperatures above 300 K.

    Fig. 3. Relationship between ΔL and T of composite multi-layered flake gun propellant (263-283 K).

    Fig. 4. Relationship between ΔL and T of composite multi-layered flake gun propellant (233-323 K).

    According to the mechanism of the interaction of nitrocellulose and nitroglycerin, the essence of the interaction between the nitrocellulose macromolecule and nitroglycerin molecule was that nitroglycerin swelled or dissolved nitrocellulose. The gravitational force between nitrocellulose and nitroglycerin decreased with increasing temperature. Below 300 K, the nitrocellulose macromolecular chain froze. Low molecular energy was insufficient to overcome the molecular bond rotation of barrier potential,in which only small molecular units could move.Therefore,it was unable to achieve conformational shift. Above 300 K, most of the molecular chain started to move,and the molecular kinetic energy increased and overcame the molecular bond rotation of barrier potential. It was found that even some part of the molecular chain slipped.This phenomenon was realized with a decrease in stiffness. As a result,the samples turned soft and were not sufficiently stiff to withstand the probe pressure of 0.05 N. Meanwhile, the probe was partially embedded into the sample,and showed a false negative expansion above 300 K.Therefore,there was an obvious turning point on ΔL-T curves.

    The transition temperature decreased with increasing thickness within the range of 0.3-0.7 mm. However, the transition temperature of No.5, No.6, and No.7 was approximately 300 K due to the nitrocellulose molecule in thinner samples being rebounded greatly. Hence, more energy was required to make the molecular chain move. Therefore, the transition temperature decreased with increasing thickness. When the thickness of specimen was greater than 0.7 mm, the content of solvent in the specimens was comparable. Due to the addition of solvent, the swelling phenomenon occurred, which increased the distance of nitrocellulose molecule,and the thicker specimens exhibited similar free volume fraction with thinner specimens. Therefore, No.5, No.6, and No.7 had the same transition temperatures.

    4.1.1. DMA test results

    In Fig. 5, the reference temperatures for loss modulus and loss factor were 223.42 K and 233.59 K, respectively. It is well-known that, the glass transition temperature represents the temperature region of phase transition for materials. Therefore, glass transition temperature is not a strict point. This temperature region can be defined by peak temperature of loss modulus and peak temperature of loss factor as lower and upper limits.This temperature range is the usage temperature limit for the material.

    According to Eq.(6),DMA test curves were modified by thermal expansion coefficient. Fig. 6 shows the modified curve for storage modulus of No.2.As temperature increased, the thickness of specimen also increased above 213 K, and the geometrical factor decreased. Therefore, the modified storage modulus value was lower than the test value.

    Fig. 5. DMA test result for No.2.

    Fig. 6. Test curve and corresponding modified curve for storage modulus of No.2.

    Fig. 7. Loss modulus curves for different samples.

    Fig. 8. Storage modulus curves for different samples.

    The loss modulus and temperature test curves are shown in Fig. 7. It is worth noting that there were two obvious peaks in the loss modulus curves,which indicated the occurrence of two phase transition processes during the increase of composite multi-layered flake gun propellant within the temperature range of 193-373 K.The first peak temperature corresponds to the glass transition temperature of the energetic material, whereas the second peak temperature occurred due to high elastic transition.Table 3 shows the corresponding glass transition temperature of composite multilayered flake gun propellant with different thickness values,where the glass transition temperature decreased with increasing thickness of composite multi-layered flake gun propellant. This phenomenon may be caused by thermal expansion of specimens due to changes in geometrical dimensions with increasing temperature.The intermolecular distance increased with increasing temperature. According to the free volume theory, the energy required for nitrocellulose molecule movement was relatively low. Based upon the results presented in Table 1, the linear thermal expansion coefficient increased with increasing thickness. Therefore, as the thickness of specimens increased,the glass transition temperature decreased.

    Moreover,a smaller peak was observed in the loss modulus and temperature curve at approx.336.07 K,which was attributed to the secondary transition temperature. Fig. 8 shows that the storage modulus decreased with increasing temperature. As the temperature increased, the molecular motor activity increased, and the spacing between the molecules increased. Additionally, the stiffness of this energetic material and value of storage modulus decreased. Furthermore, the fluidity of the energetic material was enhanced,and the loss modulus became the dominant factor above 336.07 K.This resulted in a small peak being observed in the curves.

    As described in the theoretical background, it was confirmed that the linear thermal expansion coefficient could be used for the modification of DMA test results.The left curves in Fig.9 show the unmodified loss modulus curves and right curves the modified loss modulus curves, and the peak temperatures are shown in Table 4.The first peak temperature in the corresponding curves was comparable, which meant that the glass transition temperature was Due to the frequency sweep being performed at a certain temperature, the effect of thermal expansion was eliminated, hence, the curve and modified E’and temperature curve were mirror images.

    Fig. 9. Unmodified and modified loss modulus curves of different samples.

    Table 4 Glass transition temperature of composite multi-layered flake gun propellant.

    Fig.11 shows the mirror images for E′and temperature and E’and frequency curves,where the storage modulus was the same at 253 K (on left curve) and 1 Hz (on right curve). Under the conditions of 293 K and 0.1 Hz, or 213 K and 10 Hz, the storage moduli were equal.Hence,the storage modulus curve modified by thermal expansion coefficient was proven to be effective.

    5. Conclusions

    (1) A sample of composite multi-layered flake gun propellant has been tested using TMA and DMA methods. The linear thermal expansion coefficient has been determined and lays within the range of 0.1800 × 10-4K-1- 0.2081 × 10-4K-1with the temperature range of 263-283 K.The coefficient of thermal expansion and impact on the use of material is small.DMA test shows that the glass transition temperature of the composite multi-layered flake gun propellant is 222.94 K. This result was obtained by eliminating the effect of thermal expansion.

    Fig.11. Mirror images of E′ and temperature and E′ and frequency curves.

    All the modified glass transition temperatures are presented in Table 4. The modified values were lower than the test values and that the modified value was 222.94 K. This temperature was regarded as the lower limit temperature for the energetic material.

    4.1.2. Verification of modified results

    According to the time-temperature equivalence principle, the same mechanical relaxation phenomenon of polymer can be observed at higher temperature and shorter time (or higher frequency), or at lower temperature and longer time (or lower frequency), which means that E′and temperature curve and E′and frequency curve are mirror images. Based upon this theory, frequency sweep of this specimen was performed at a certain temperature, and E′and frequency curve was obtained. E′and frequency curve and E’ and temperature curve were found to be mirror images.

    The frequency sweep experiments were carried out at 293 K within the frequency range of 0.01-10 Hz. The amplitude of test sample's deformation was established at 10 μm.The obtained curve is shown in Fig.10.

    Fig.10. Modified loss modulus curves for different samples.

    Fig. 10 shows the frequency and E′curve for No.2 specimens.223.01-223.50 K, which conforms to the double-based propellant.

    (2) According to DMA test, geometric dimensions have a great impact on the storage and loss moduli. Furthermore, as the geometric dimensions increase, the loss modulus decreases.Therefore, the thermal expansion in the test process has a certain influence on the storage and loss modulus.According to the time-temperature equivalence principle, this method is an effective and reliable method to modify the storage and loss modulus.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    99国产精品一区二区蜜桃av| 可以在线观看的亚洲视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美人与善性xxx| 午夜福利在线观看免费完整高清在 | 中国美女看黄片| 色视频www国产| 免费一级毛片在线播放高清视频| 永久网站在线| 国产av不卡久久| videossex国产| 人妻久久中文字幕网| 91狼人影院| 97碰自拍视频| 免费看美女性在线毛片视频| 韩国av在线不卡| 窝窝影院91人妻| 夜夜看夜夜爽夜夜摸| 久久久成人免费电影| 色综合站精品国产| 日本撒尿小便嘘嘘汇集6| 男女那种视频在线观看| 身体一侧抽搐| 免费不卡的大黄色大毛片视频在线观看 | 日本a在线网址| 97超级碰碰碰精品色视频在线观看| 日日夜夜操网爽| 日本成人三级电影网站| 亚洲av美国av| 久99久视频精品免费| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 国产蜜桃级精品一区二区三区| 丝袜美腿在线中文| 内地一区二区视频在线| 国产视频一区二区在线看| 国产精品av视频在线免费观看| 69人妻影院| 欧美3d第一页| 欧美不卡视频在线免费观看| 国产免费av片在线观看野外av| 夜夜爽天天搞| 久久人妻av系列| 精品不卡国产一区二区三区| 精品一区二区免费观看| 午夜福利视频1000在线观看| 天堂网av新在线| 精品一区二区三区视频在线观看免费| 床上黄色一级片| 午夜精品在线福利| 少妇的逼水好多| 欧美+亚洲+日韩+国产| 久久久久久国产a免费观看| 能在线免费观看的黄片| 91在线观看av| www日本黄色视频网| 男女之事视频高清在线观看| 亚洲图色成人| 九九在线视频观看精品| 九色成人免费人妻av| 最近在线观看免费完整版| 精品一区二区三区av网在线观看| 波多野结衣高清无吗| 九九热线精品视视频播放| 少妇人妻精品综合一区二区 | 午夜激情欧美在线| 久久久久久久久久成人| 久久亚洲精品不卡| 春色校园在线视频观看| 国产精品一区二区免费欧美| 22中文网久久字幕| 国内毛片毛片毛片毛片毛片| 九九爱精品视频在线观看| 日本黄色视频三级网站网址| 亚洲av免费在线观看| 成年人黄色毛片网站| 亚洲在线观看片| 我的老师免费观看完整版| 国产亚洲91精品色在线| 久久久久性生活片| 国产毛片a区久久久久| 一区福利在线观看| 欧美黑人巨大hd| 别揉我奶头 嗯啊视频| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久久亚洲 | 欧美潮喷喷水| 国产伦一二天堂av在线观看| 白带黄色成豆腐渣| 不卡一级毛片| 自拍偷自拍亚洲精品老妇| 亚洲美女视频黄频| 九色成人免费人妻av| 久久6这里有精品| 男插女下体视频免费在线播放| 热99re8久久精品国产| 女的被弄到高潮叫床怎么办 | 国产精品免费一区二区三区在线| 国产精品爽爽va在线观看网站| 免费观看的影片在线观看| 精品午夜福利视频在线观看一区| 国产黄片美女视频| 国产男人的电影天堂91| 啦啦啦啦在线视频资源| 黄色日韩在线| 亚洲av美国av| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | 自拍偷自拍亚洲精品老妇| 嫩草影院新地址| 国产精品久久电影中文字幕| 国产视频一区二区在线看| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 亚洲,欧美,日韩| 精品久久久久久久久亚洲 | 在线观看av片永久免费下载| 国产伦人伦偷精品视频| 观看免费一级毛片| 最好的美女福利视频网| 桃红色精品国产亚洲av| 18+在线观看网站| 亚洲成人精品中文字幕电影| 日本 欧美在线| 久久人妻av系列| ponron亚洲| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 免费看光身美女| 精品一区二区免费观看| 国产精品嫩草影院av在线观看 | 亚洲自偷自拍三级| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品合色在线| 欧美精品国产亚洲| 成人亚洲精品av一区二区| 小说图片视频综合网站| 亚洲av成人av| 美女 人体艺术 gogo| 天天一区二区日本电影三级| 哪里可以看免费的av片| 老司机福利观看| 成人特级黄色片久久久久久久| 女生性感内裤真人,穿戴方法视频| 夜夜夜夜夜久久久久| 少妇高潮的动态图| 国产黄a三级三级三级人| 亚洲va在线va天堂va国产| 亚洲一区高清亚洲精品| 日日夜夜操网爽| 欧美日韩综合久久久久久 | 毛片女人毛片| 日韩,欧美,国产一区二区三区 | 中文资源天堂在线| 国产精品国产三级国产av玫瑰| 别揉我奶头~嗯~啊~动态视频| 色5月婷婷丁香| 草草在线视频免费看| 国产亚洲av嫩草精品影院| 中文字幕高清在线视频| 久久精品国产清高在天天线| 99精品久久久久人妻精品| 欧美+日韩+精品| 国产主播在线观看一区二区| 午夜爱爱视频在线播放| 波野结衣二区三区在线| 亚洲精品国产成人久久av| 精品午夜福利视频在线观看一区| 国产乱人视频| 欧美xxxx黑人xx丫x性爽| 精品国内亚洲2022精品成人| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 人妻制服诱惑在线中文字幕| 亚洲熟妇中文字幕五十中出| 91麻豆av在线| 日本免费一区二区三区高清不卡| 欧美性感艳星| av女优亚洲男人天堂| 99久久成人亚洲精品观看| 天天一区二区日本电影三级| 午夜福利高清视频| 天堂动漫精品| 国产精华一区二区三区| 亚洲第一电影网av| 国产毛片a区久久久久| 日韩欧美精品v在线| 极品教师在线免费播放| 国内精品久久久久久久电影| 国产av麻豆久久久久久久| 国产女主播在线喷水免费视频网站 | 亚洲欧美日韩高清在线视频| 亚洲中文字幕一区二区三区有码在线看| 他把我摸到了高潮在线观看| 非洲黑人性xxxx精品又粗又长| 免费无遮挡裸体视频| 国产国拍精品亚洲av在线观看| 日本黄色片子视频| 精品一区二区三区人妻视频| 亚洲黑人精品在线| 亚洲无线观看免费| 天堂动漫精品| 精品午夜福利视频在线观看一区| 九九爱精品视频在线观看| 日韩中文字幕欧美一区二区| 天美传媒精品一区二区| 可以在线观看毛片的网站| 97人妻精品一区二区三区麻豆| 午夜免费激情av| 尾随美女入室| 国产欧美日韩精品亚洲av| 蜜桃亚洲精品一区二区三区| 精品一区二区三区人妻视频| 最近最新免费中文字幕在线| 看黄色毛片网站| 草草在线视频免费看| 在线免费观看不下载黄p国产 | 简卡轻食公司| 在线天堂最新版资源| 最近视频中文字幕2019在线8| 内地一区二区视频在线| 国产三级在线视频| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区三区四区免费观看 | 九九热线精品视视频播放| 直男gayav资源| 亚洲av中文av极速乱 | 日本免费一区二区三区高清不卡| 一本久久中文字幕| 精品福利观看| 我要看日韩黄色一级片| 成人欧美大片| 91麻豆精品激情在线观看国产| 色视频www国产| 全区人妻精品视频| 国产高清有码在线观看视频| av专区在线播放| avwww免费| 国产成人av教育| 黄色一级大片看看| a级毛片a级免费在线| 国产又黄又爽又无遮挡在线| 亚洲精品日韩av片在线观看| 国产美女午夜福利| 国产一区二区三区视频了| 久久6这里有精品| 小蜜桃在线观看免费完整版高清| 国产精品一区二区三区四区久久| 一本精品99久久精品77| 精品国产三级普通话版| 啦啦啦啦在线视频资源| 日韩欧美国产一区二区入口| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看| 日本a在线网址| 天堂av国产一区二区熟女人妻| 国产一区二区在线观看日韩| 91久久精品国产一区二区三区| 国产精品综合久久久久久久免费| 亚州av有码| 一个人免费在线观看电影| 国产成人av教育| 91狼人影院| 两个人的视频大全免费| 免费高清视频大片| 乱系列少妇在线播放| 日日摸夜夜添夜夜添av毛片 | 一进一出好大好爽视频| 一个人观看的视频www高清免费观看| 亚洲男人的天堂狠狠| 在线国产一区二区在线| 久久欧美精品欧美久久欧美| 国产精品女同一区二区软件 | 中文字幕熟女人妻在线| 国产在线精品亚洲第一网站| 亚洲av.av天堂| 91精品国产九色| 精品人妻偷拍中文字幕| 国产精品嫩草影院av在线观看 | 制服丝袜大香蕉在线| 欧美潮喷喷水| 99在线人妻在线中文字幕| 尤物成人国产欧美一区二区三区| 少妇人妻一区二区三区视频| 亚洲一区高清亚洲精品| 国产精品福利在线免费观看| 欧美色欧美亚洲另类二区| 99久久久亚洲精品蜜臀av| 久久久久国内视频| 精品不卡国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 欧美高清成人免费视频www| 亚洲av成人av| 亚洲三级黄色毛片| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 亚州av有码| 直男gayav资源| www日本黄色视频网| 在线看三级毛片| ponron亚洲| 亚洲四区av| 小说图片视频综合网站| 在线观看午夜福利视频| 小说图片视频综合网站| 成人无遮挡网站| 人妻制服诱惑在线中文字幕| 日本熟妇午夜| 欧美日韩中文字幕国产精品一区二区三区| 国产综合懂色| 国产伦精品一区二区三区视频9| 女人被狂操c到高潮| 大型黄色视频在线免费观看| 一进一出抽搐动态| 欧美+亚洲+日韩+国产| 好男人在线观看高清免费视频| 又黄又爽又免费观看的视频| 日本五十路高清| 波野结衣二区三区在线| 免费av毛片视频| 国产淫片久久久久久久久| 小说图片视频综合网站| www.www免费av| 看免费成人av毛片| 亚洲aⅴ乱码一区二区在线播放| 亚洲无线观看免费| 在线观看美女被高潮喷水网站| 18+在线观看网站| 久久久久久国产a免费观看| 精品久久久噜噜| 久久久久久久久大av| 久久精品国产亚洲av天美| 久久精品久久久久久噜噜老黄 | 国产精品无大码| 很黄的视频免费| 欧美一区二区精品小视频在线| 毛片女人毛片| 国产免费一级a男人的天堂| 99久久精品国产国产毛片| 色综合婷婷激情| 国产高清激情床上av| 免费观看的影片在线观看| 亚洲国产日韩欧美精品在线观看| 香蕉av资源在线| 九九爱精品视频在线观看| 欧美bdsm另类| 一级黄片播放器| 九九在线视频观看精品| 亚洲国产欧美人成| 成人亚洲精品av一区二区| 日韩精品青青久久久久久| 18禁黄网站禁片午夜丰满| 小说图片视频综合网站| 国模一区二区三区四区视频| 国产老妇女一区| 中国美女看黄片| 国产黄片美女视频| 啪啪无遮挡十八禁网站| 国产精品国产高清国产av| 亚洲七黄色美女视频| 日本熟妇午夜| 国产精品日韩av在线免费观看| 人妻丰满熟妇av一区二区三区| 亚洲自偷自拍三级| 成人欧美大片| 中文字幕熟女人妻在线| 国产伦一二天堂av在线观看| 欧美成人性av电影在线观看| 一个人看的www免费观看视频| 91久久精品国产一区二区成人| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av第一区精品v没综合| 国产成人一区二区在线| 亚洲七黄色美女视频| 看免费成人av毛片| 久久天躁狠狠躁夜夜2o2o| 国产精品1区2区在线观看.| 最近最新免费中文字幕在线| 在线天堂最新版资源| 男插女下体视频免费在线播放| 欧美三级亚洲精品| 久99久视频精品免费| 大型黄色视频在线免费观看| 夜夜夜夜夜久久久久| 亚洲成人精品中文字幕电影| 中文字幕免费在线视频6| 3wmmmm亚洲av在线观看| 欧美国产日韩亚洲一区| 国产成年人精品一区二区| 亚洲人成网站高清观看| 精品久久国产蜜桃| 欧美日本视频| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 日本 欧美在线| 国产精品亚洲一级av第二区| 欧美色视频一区免费| 我的老师免费观看完整版| 午夜爱爱视频在线播放| 最好的美女福利视频网| 午夜福利成人在线免费观看| 久久精品国产清高在天天线| 啦啦啦韩国在线观看视频| 久久久久久国产a免费观看| 人人妻人人看人人澡| 俄罗斯特黄特色一大片| 国产精品一区二区三区四区免费观看 | 久久久久国内视频| 悠悠久久av| 三级毛片av免费| 日本与韩国留学比较| 蜜桃亚洲精品一区二区三区| 精品一区二区免费观看| 成人综合一区亚洲| 嫩草影院精品99| 亚洲,欧美,日韩| 亚洲欧美日韩卡通动漫| 亚洲成人免费电影在线观看| 国产真实乱freesex| 久久久午夜欧美精品| 99热这里只有是精品在线观看| 亚洲国产精品sss在线观看| 欧美另类亚洲清纯唯美| 高清日韩中文字幕在线| 成人一区二区视频在线观看| 极品教师在线免费播放| 看十八女毛片水多多多| 男人和女人高潮做爰伦理| 日本免费一区二区三区高清不卡| 干丝袜人妻中文字幕| 男女下面进入的视频免费午夜| av黄色大香蕉| 精品久久国产蜜桃| 久久精品国产亚洲av涩爱 | 免费观看人在逋| 久久久久国产精品人妻aⅴ院| 一区二区三区免费毛片| 1000部很黄的大片| 日本黄色片子视频| 久久精品国产亚洲av香蕉五月| 天堂影院成人在线观看| 亚洲内射少妇av| 在线观看舔阴道视频| av视频在线观看入口| 女同久久另类99精品国产91| 亚洲午夜理论影院| 欧美精品国产亚洲| 91狼人影院| 99久久无色码亚洲精品果冻| 国产一级毛片七仙女欲春2| 亚洲中文字幕日韩| 可以在线观看的亚洲视频| 亚洲乱码一区二区免费版| 成年女人毛片免费观看观看9| 草草在线视频免费看| 成人美女网站在线观看视频| 亚洲精品国产成人久久av| 不卡视频在线观看欧美| 欧美激情在线99| 在线天堂最新版资源| 可以在线观看的亚洲视频| 国产一区二区三区av在线 | 国产单亲对白刺激| 亚洲av.av天堂| 日韩一区二区视频免费看| 欧美性猛交╳xxx乱大交人| 伦理电影大哥的女人| 日本欧美国产在线视频| 最好的美女福利视频网| 久久精品国产鲁丝片午夜精品 | 国产精品99久久久久久久久| 亚洲精品影视一区二区三区av| 身体一侧抽搐| 免费看av在线观看网站| www.色视频.com| 国内精品久久久久精免费| 露出奶头的视频| 成人国产综合亚洲| 麻豆国产av国片精品| 午夜影院日韩av| 亚洲第一电影网av| 亚洲av日韩精品久久久久久密| 亚洲,欧美,日韩| 国模一区二区三区四区视频| 97碰自拍视频| 欧美另类亚洲清纯唯美| 国产国拍精品亚洲av在线观看| 伦精品一区二区三区| 欧美成人a在线观看| 12—13女人毛片做爰片一| 国产单亲对白刺激| 亚洲av免费在线观看| 国产亚洲精品综合一区在线观看| 色精品久久人妻99蜜桃| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 国产高清视频在线观看网站| 久久久久久久久久成人| 99久久精品国产国产毛片| 午夜免费成人在线视频| 黄色丝袜av网址大全| 色综合站精品国产| 国语自产精品视频在线第100页| 精品一区二区三区视频在线| 91av网一区二区| 男女那种视频在线观看| 我要搜黄色片| 在线观看免费视频日本深夜| 可以在线观看的亚洲视频| 国产一区二区在线观看日韩| 欧美最黄视频在线播放免费| 亚洲欧美日韩卡通动漫| 12—13女人毛片做爰片一| 中出人妻视频一区二区| 亚洲熟妇中文字幕五十中出| 伊人久久精品亚洲午夜| 久久午夜亚洲精品久久| 97人妻精品一区二区三区麻豆| 亚洲人成网站高清观看| 国产伦人伦偷精品视频| 亚洲一区二区三区色噜噜| 欧美中文日本在线观看视频| 亚洲一级一片aⅴ在线观看| 亚洲午夜理论影院| 俺也久久电影网| 人人妻人人澡欧美一区二区| 真实男女啪啪啪动态图| 国产精品女同一区二区软件 | 别揉我奶头~嗯~啊~动态视频| 成人国产一区最新在线观看| 在线免费观看不下载黄p国产 | 啦啦啦啦在线视频资源| 国产真实伦视频高清在线观看 | 无遮挡黄片免费观看| 91麻豆av在线| 久久久成人免费电影| 天堂网av新在线| 成人特级av手机在线观看| 亚洲无线观看免费| 欧美黑人巨大hd| 亚洲一区二区三区色噜噜| 国产美女午夜福利| 亚洲精品一区av在线观看| 成人亚洲精品av一区二区| 18禁黄网站禁片免费观看直播| 精品日产1卡2卡| 国产高清激情床上av| av在线老鸭窝| 九九热线精品视视频播放| 免费看光身美女| 亚洲精品久久国产高清桃花| 我的老师免费观看完整版| 51国产日韩欧美| 最好的美女福利视频网| 人妻久久中文字幕网| 久久久成人免费电影| 国产色爽女视频免费观看| 久久精品综合一区二区三区| 国产亚洲精品av在线| 亚洲aⅴ乱码一区二区在线播放| 天天一区二区日本电影三级| 日本 欧美在线| 国产极品精品免费视频能看的| 狂野欧美白嫩少妇大欣赏| 亚洲精品成人久久久久久| 日韩欧美国产在线观看| 女生性感内裤真人,穿戴方法视频| 欧美人与善性xxx| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 麻豆国产av国片精品| 午夜免费成人在线视频| 一级av片app| 精品国产三级普通话版| 尤物成人国产欧美一区二区三区| 亚洲av免费在线观看| 欧美日本视频| 亚洲精品国产成人久久av| 欧美最黄视频在线播放免费| 狂野欧美激情性xxxx在线观看| netflix在线观看网站| 美女xxoo啪啪120秒动态图| 真人做人爱边吃奶动态| 国产精品精品国产色婷婷| 网址你懂的国产日韩在线| 男人和女人高潮做爰伦理| 午夜福利在线在线| 精品午夜福利视频在线观看一区| 亚洲精品日韩av片在线观看| 级片在线观看| 床上黄色一级片| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 又粗又爽又猛毛片免费看| 女的被弄到高潮叫床怎么办 | 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲91精品色在线| 欧美+日韩+精品| 久久天躁狠狠躁夜夜2o2o| 成人毛片a级毛片在线播放| 91精品国产九色| 亚洲av免费在线观看| 日韩欧美在线乱码| а√天堂www在线а√下载| 国产高清不卡午夜福利| 乱系列少妇在线播放| 少妇人妻一区二区三区视频| 校园春色视频在线观看| 国产男人的电影天堂91| 18+在线观看网站| 色播亚洲综合网| 成年女人看的毛片在线观看| 很黄的视频免费| 男人舔女人下体高潮全视频| 精品国内亚洲2022精品成人|