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

    Long-term corona behaviour and performance enhancing mechanism of SiC/epoxy nanocomposite in SF6 gas environment

    2023-03-15 00:54:24JingruiWANG王靖瑞QingminLI李慶民YanfengGONG公衍峰QixinHOU侯啟鑫HengLIU劉衡JianWANG王健andHanwenREN任瀚文
    Plasma Science and Technology 2023年3期
    關(guān)鍵詞:王健

    Jingrui WANG (王靖瑞), Qingmin LI (李慶民), Yanfeng GONG (公衍峰),Qixin HOU (侯啟鑫), Heng LIU (劉衡), Jian WANG (王健) and Hanwen REN (任瀚文)

    State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources,North China Electric Power University, Beijing 102206, People’s Republic of China

    Abstract Surface coating technology is an effective way to solve the interface insulation problem of DC GIS/GIL basin insulators, but the performance of the coating will change greatly, and the insulation strength will be completely lost, after long-term use in the extreme conditions of corona erosion.In this research, the multi-needle-plate electrode platform was constructed to explore the long-term use performance of SiC-doped nanocomposite exposed to corona discharge in SF6 gas.Samples with a high SiC content have advantages in maintaining physical and chemical properties such as elemental composition, erosion depth, surface roughness and mass loss.The nanocomposite doped with 6 wt.%SiC has prominent surface insulation strength after long term exposure to corona, and the others are close to losing, or have completely lost,their insulating properties.Furthermore, the degradation mechanism of physicochemical properties of composite exposed to corona discharge was investigated with the proposed ReaxFF MD model of energetic particles from SF6 decomposition bombarding the epoxy surface.The reaction process of SF particles and F particles with the cross-linked epoxy resin, and the SiC nanoparticles providing shelter to the surrounding polymer and mitigating their suffering direct bombardment, have been established.The damage propagation depth, mass loss and surface roughness change of nanocomposite material bombarded by SF6 decomposition products is reproduced in this simulation.Finally, the deterioration mechanism of insulation properties for the SiC-doped composite was elucidated with DFT analysis.The band gap of the molecule containing S drops directly from the initial 7.785 eV to 1.875 eV,which causes the deterioration of surface electric properties.

    Keywords:corona discharge in SF6,SiC doped nanocomposite,epoxy resin,ReaxFF MD,DFT

    1.Introduction

    High voltage direct current (HVDC) technology is an important support for renewable energy consumption and the stable operation of new power systems [1-4].The DC gas insulated transmission line (DC GIL) and DC gas insulated switchgear (DC GIS), as alternatives for traditional equipment, are the most important technologies in future power systems.The gas-solid insulation system still faces two severe challenges [5, 6].On the one hand, the charge accumulation on the gas-solid interface under DC voltage leads to serious electric field distortion.On the other hand,metal particles and micro-nano dust produced for multiple reasons will be driven by a DC electric field and adsorbed on the spacer surface.The interaction of these two factors will more easily cause corona discharge and partial discharge,and even insulation failure and severe ablation on the surface [7].

    Recent research shows that surface coating technology is one of the most practical methods to alleviate these surface insulation problems[8].Various inorganic nanofillers,such as TiO2, BN, montmorillonite, etc, were applied to modify the coating composites to regulate surface charge [9].Compared with them, silicon carbide (SiC), which has a nonlinear conductivity property, could make the conductivity of epoxy resin composites adapt to the applied electric field [10, 11].Scholars have designed nano or micron SiC/epoxy coatings with nonlinear conductivity for the DC insulator, which can control the electric field distribution pattern along the surface,inhibit dust adsorption and improve the surface insulation strength[12-14].Therefore,the SiC-doped composite coating is regarded as the most promising method to enhance the surface insulation of DC basin insulators.

    The basic electrical properties of the SiC coating surface at the early stage of service in SF6gas have been studied.However,there are few studies on the performance change of SiC composite coatings under harsh conditions for long periods.Due to the inevitable existence of corona discharge and partial discharge inside the equipment,SF6would decompose into low fluoride sulfide, fluorine atoms, and so on [7, 15].These particles generated by discharge are not only highly reactive but also in a high temperature state of more than 8000 K, and they would rapidly impact the coating surface with the equivalent energy of 0.65-1.03 eV [16, 17].Under such harsh conditions,a SiC/epoxy coating will be destroyed by chemically active gas by-products, and the microstructure and chemical composition will be greatly changed, and even some unpredictable insulation failures will occur.Although some studies have pointed out that anti-corona materials modified by nano-SiC are used in fusion reactor coils[18],the operating conditions are totally different from the SF6environment.The existing research on the influence of corona discharge on the treated surface,such as a directly fluorinated surface, mainly use corona discharge model experiments to directly characterize the influence of discharge on epoxy surface deterioration and the possible functions of the treatment [7, 19].However, it is still difficult to comprehensively understand the influence mechanism of SF6decomposition products on the surface dissociation of epoxy polymer, so a targeted performance improvement strategy cannot be formed.

    Some representative simulation models have been developed to analyze the effects of particle bombardment on another material interface.Rahmani and van Duin [20, 21]studied the degradation mechanism of polymers, such as epoxy and polyimide,under the impact of atomic oxygen in a low earth orbit based on reactive force field molecular dynamic (ReaxFF MD) simulation.Bai et al [22] carried out MD simulation of the phase transformation of fused silica under nanoparticle impact in high power laser systems.Zhao et al [23] elucidated the microscopic mechanism of the interaction between OH radicals in plasma with biofilms using MD simulation methodology.Based on these studies,the model was further developed and applied to discharge fault simulation.Wang et al[24]transformed the dissociation of bisphenol F epoxy resin under partial discharge into the crack process induced by the impact of high-energy particles.Huang et al[25]have applied this method to study polyimide damage caused by DC corona in the air.The ReaxFF MD method provides an effective way to study the decomposition mechanism of SiC-doped composite under the high energetic particle’s bombardment.

    The prominent corona resistance characteristics and longterm service stability are the key to supporting the application of SiC-doped nanocoating.In this research, SiC/epoxy nanocomposite was prepared.The platform for a polymer exposed to corona discharge under electrothermal coupling was established, and the surface damage characteristics of nano-SiC/epoxy composite under corona erosion in SF6gas were realized in the experimental study.Furthermore, the change rules of microstructure, elemental composition, surface roughness, mass loss, and surface insulation strength of SiC/epoxy composite exposed to corona for different times were analyzed.Finally, atomic-scale simulation analysis was carried out based on ReaxFF MD to reveal the mechanism of corona erosion-induced changes in the physicochemical properties of the composite surface, and DFT calculations were performed to reveal the microscopic mechanism of longterm corona-induced surface insulation degradation and insulation failure.

    2.Sample preparation and key performance test

    2.1.Preparation of SiC/epoxy nanocomposite

    Diglycidyl ether of Bisphenol A (DGEBA) epoxy resin,methyl tetrahydrophthalic anhydride (MTHPA) and the accelerant, 2,4,6-dimethylaminomethylphenol (DMP-30)were used to prepare the polymer specimens.The nano-SiC with a purity of 99.9%and a particle size of 40 nm is used as the dopant.The preparation of pure epoxy resin polymer and SiC/EP nanocomposite is carried out according to the mature craftsmanship [26], the thickness of the obtained samples is about 1 mm.Nano-SiC doped composites (abbreviated SEC)containing fillers in five different proportions,2 wt.%,4 wt.%,6 wt.%, 8 wt.% and 10 wt.%, were prepared.The SiC/epoxy composite with a mass fraction of 2 wt.% is taken as an example, which is abbreviated as SEC2 in the following description.

    2.2.Experimental platform for composite exposed to corona discharge

    In order to further study the influence of long-term corona discharge on physicochemical properties and the insulation performance of the composite surface, the experimental platform for polymer exposed to corona discharge under electrothermal coupling was designed and developed in this research,and its main structure diagram is shown in figure 1.The platform is mainly composed of a high voltage power supply, bushing, a seal chamber, an electrode system, an oil heating system, and a gas supply and disposal system.The multi-needle-plate electrode system can simulate the tip discharge fault on the high voltage conductor near the insulator in GIS/GIL,and it can be ensured that each sample is facing three needle electrodes and the distance between the sample and the tip of the needle electrode is 30 mm.For each needle electrode, the radius of curvature of the needle tip is about 100 μm.When the applied voltage is +45 kV, the schematic diagram of electric field distribution between the electrodes is shown in the cloud diagram of figure 1.The electric field strength at the needle tip could be calculated using Mason’s equation [27] as equation (1) and the maximum value of its magnitude could reach 127 kV mm-1, which is much larger than the onset field of positive corona calculated by Nitta’s model [28].

    Figure 1.Schematic diagram of experimental platform for polymer exposed to corona discharge under electrothermal coupling.

    where U is applied voltage;d is distance from the plate to the needle tip; R is the curvature radius at the tip.

    The corona discharge generated by the needle electrode directly sputtered onto the epoxy composite samples’ surface below the electrode.The maximum allowed pressure of the seal test system could reach 0.6 MPa.The temperature inside the chamber can be controlled between 50°C and 120°C by the external oil circulation heating system to simulate the working conditions of the gas insulated equipment, and test experiments were carried out at 90 °C.Due to the toxic components of decomposition products, the sealing property of the chamber should be guaranteed reliably,and the tail gas disposal system can remove harmful substances and recover SF6gas.The experiments were carried out in 0.2 MPa SF6gas for 10 h, 20 h and 30 h respectively.After that, the key physical and chemical properties of the samples were compared and analysed.

    2.3.Surface insulation strength and surface potential decay test platform

    Since the SiC/epoxy composite is a coating material applied to the DC insulator’s surface, the changes in surface insulation strength and surface potential decay after long-term exposure are the keys to evaluating the feasibility of its application.

    The surface flashover test platform shown in figure 2(a)was applied to evaluate the critical surface insulation strength after corona erosion.The flashover test was conducted with finger electrodes in the seal chamber with the SF6gas pressure of 0.1 MPa.The radius of curvature of the hemisphere of the two finger electrodes is 8 mm, and the distance between the two electrodes is 5 mm.Figure 2(b) gives the charge injection and dissipation test system.The sample is placed on the movable platform, and the electrode system consists of a needle, a metal net and a grounded plate.During the experiment, the sample was charged under the corona with the DC voltage of +10 kV for 10 min, and then the power supply was turned off and the sample was moved below the probe to perform measurements of the surface potential change.In each experiment, a data acquisition card was used to record the change of the surface potential within 1800 s.The test experiments were conducted at an ambient temperature from 20°C to 25°C and an air humidity of about 20%.

    Figure 2.Schematic diagram of experimental platform for key electric properties.(a)Surface flashover test platform,(b)charge injection and dissipation test system.

    3.Variation of key physicochemical and electric properties after long-term service in corona discharge

    3.1.Microstructure and elemental composition on eroded area

    The surface morphology and element composition of the exposed area of the composite sample obtained by the scanning electron microscope(SEM)at their initial state and after being exposed to corona for 30 h are shown in figure 3.The SEM of the surface and cross section of pure epoxy resin are given in figures 3(a) and (b), and figures 3(f) and (g) show these states of SEC10, which are all smooth and clean.As shown in figure 3(c),after being exposed to corona for a long time, the pure epoxy surface was obviously damaged, and many defects like cracked bubbles appeared.However,taking SEC10 as an example, after the same exposure time, the sample surface shown in figure 3(h) is still relatively smooth with only local imprint.Some important characteristics related to erosion behavior can also be found from the cross sections of the sample.After 30 h exposure, obvious stratification phenomena appeared on the cross sections of both samples.Figures 3(d)and(i)show the specific characteristics,and the surface material thickness of the pure epoxy resin sample and SEC10 is about 24.8 μm and 7.7 μm, respectively.The energy dispersive x-ray spectroscopy (EDX) of these two samples’ surface was further analyzed.Comparing figures 3(e) and (j), it can be found that the characteristic peaks of S and F elements appear on both surfaces after 30 h exposure under corona.This indicates that there is physical and chemical reaction between the exposed surface molecules and the decomposed S- or F-containing groups under the erosion of corona discharge, which stay on the composite surface,penetrate to the deep area and promote the formation of the surface erosion layer.Compared with pure epoxy, the peak of S on the SEC10 surface is lower and the peak of F is similar.

    Figure 3.SEM and EDX of epoxy composite initially and after being exposed to corona discharge for 30 h.

    To further analyze the different effects of SiC content,the surface layer thickness was defined as the corona erosion depth.The depth of six erosion sites is selected and averaged for each sample, and the average values with different components exposed to corona for 30 h are obtained.The average corona erosion depth of all samples after 30 h exposure is shown in figure 4.The average corona erosion depth of pure epoxy is the largest,which is 23.1 μm.The depth of composite decreases gradually with the increase of SiC mass fraction,and the erosion depth of SEC10 is the smallest at 9.8 μm.The average values of SEC6,SEC8 and SEC10 are between 10 μm and 15 μm, and the values are relatively close.

    From the analysis of SEM and EDX,it can be speculated that SiC dopants can effectively resist the behavior of SF6decomposition products eroding more seriously.

    Figure 4.The average corona erosion depth of pure epoxy and composites with different components after 30 h exposure.

    3.2.Surface roughness and mass loss

    Corona erosion will cause changes in the microscopic topography and elemental composition of the sample surface at the microscale level, and it will be directly manifested as the change of surface roughness at the macroscale level.The surface roughness tester was used to measure the surface roughness of these six kinds of composite with different components, and the value was also the average of the six measurement results, as shown in figure 5(a).The doped nanoparticles had no significant effect on the surface roughness of the composites in the initial state.When prepared,the surfaces of these materials are smooth, with a roughness in the range of 0.03-0.05 μm.While the surface roughness of all kinds of samples increases gradually with the increase of exposure time under corona.After exposed to corona for 10 h,the surface roughness of the samples of each component is similar, all around 0.2 μm.It is worth noting that after exposure for 20 h or 30 h,the roughness of the samples doped with SiC nanoparticles is smaller than that of the pure epoxy samples, and the advantages of nano-doping are gradually reflected.Among the sample exposed to corona for 30 h, the roughness value of SEC10 is the smallest.

    Figure 5.(a)Surface roughness and (b) mass change of epoxy resin composites exposed to corona discharge for different durations.

    The precision electronic balance was used to accurately weigh the mass of these samples, and the obtained sample mass changes are shown in figure 5(b).All mass changes given are the difference between the mass gained after the fixed exposure time and the initial mass.Through comparison and analysis, it can be found that after 10 h exposure, the mass of the six kinds of samples all increased slightly with little difference, all around 0.001 g.After 20 h corona exposure, the mass changes of the composite samples with different components began to show significant differences.The mass of the pure epoxy sample was similar to that of the 10 h exposure one,remaining at around 0.001 g.The mass increase of SEC2 and SEC4 is about 0.0030 g, and the increments of the other samples with higher SiC doping content reach 0.005-0.006 g, which is about twice that of the former.After 30 h corona exposure,the mass of the pure epoxy sample and SEC2 sample decreased compared with that of 30 h exposure.The former was even 0.005 g less than the initial state,and the latter was 0.0005 g less than the 20 h sample.The mass gain of SEC4 at 30 h was the same as that at 20 h exposure,remaining at 0.0020 g.The mass of other SiC-containing samples still maintains an increasing trend after exposure for 30 h.It is not difficult to see that with the increase in exposure time, only the mass of pure EP and SEC2 increased first and then decreased, while the quality of other samples kept increasing.

    From the analysis of the above mass changes,we observe that when the doping content of the SiC nanoparticles is low,SF6decomposition products with high reactivity are further converted into gaseous derivatives after reacting with epoxy resin molecules.When the SiC content was higher than 4 wt.%,the nanoparticles effectively suppressed the intensity of the chemical reaction on the surface and reduced the mass loss caused by the decomposition of the epoxy resin into gaseous products.With reference to the aforementioned EDX,it can be inferred that the heavy groups containing an S or F atom formed by SF6decomposition stay on the surface of the corona erosion zone, resulting in mass increase.

    3.3.Surface insulation strength and surface potential decay

    The DC surface flashover voltage of the composite sample obtained through experiments is shown in figure 6(a),and the surface flashover phenomenon of SEC6 after being exposed for 30 h and photographed by a high-speed camera is shown in figure 6(b).In the initial state, the DC surface flashover voltage of the composite showed a trend of first increasing and then decreasing with increasing SiC doping content.When the doping content is 6 wt.%,the DC flashover voltage is the highest, which is about 16.9% higher than that of the pure sample.The flashover voltage after corona erosion is decreased to different degrees, while the change trend of the flashover voltage of the samples with the doping ratio of SiC is consistent with that of the pure one, which is not affected by exposure time.After being exposed to the corona for 10 h,the flashover voltage of pure epoxy resin decreased by 11.9%compared with the initial state.The flashover voltage decrease percentage of SEC2 and SEC4 is similar to that of pure epoxy,and the drop of SEC6 is slightly smaller, about 9.6%.The flashover voltage of SEC8 and SEC10 has the smallest drop, which is close numerically to the original state.When exposed for 20 h,the surface flashover voltage of pure epoxy resin drops significantly, to about 75.8% lower than that of the original one.At the same time,the flashover voltage of the SiC-doped polymer drops by 50%to 60%,while the flashover voltage of the composite is still higher than that of the pure sample.In particular, the flashover voltage is 6.3 kV and 7.9 kV higher than those of SEC6 and SEC8, respectively.When the corona erosion time reaches 30 h, the insulating ability of pure epoxy resin and SEC2 has completely failed.However,the other composites still retain a certain insulation strength, and the flashover voltage of SEC6 is the highest among them.It can be determined that long-term exposure to corona discharge will lead to great deterioration of the composite surface insulation performance.

    Figure 6.DC flashover of epoxy resin after being exposed to corona.(a) Flashover voltage variation.(b) Surface flashover phenomenon.

    The pure epoxy sample and SEC6 sample are selected in the initial state and after 30 h exposure for surface potential decay analysis, and the surface potential normalized by the initial value is given in figure 7.For the samples at initial state,the potential decay on the SEC6 surface is much quicker than that for the pure epoxy,since the SiC dopants introduce a large number of shallow traps and reduce the density and energy level of the deep traps.Due to the 30 h exposure, the pure epoxy completely lost its surface insulation and could not obtain the potential decay behavior on its surface.After corona erosion, the potential decay rate of SEC6 is greatly increased.The surface electric conductivities of these four kinds of samples obtained by the three-electrode test system are given in table 1.It can be seen that long-term corona erosion will greatly decrease the electrical conductivity of the composite surface.The conductivity of pure epoxy resin increased from the initial value, 6.56 × 10-15S m-1, to 5.39 × 10-9S m-1after exposure, and this abnormal increase is the main reason for the surface insulation failure.The increase in the conductivity of SEC6 directly leads to the rise in the potential decay rate.

    Figure 7.Normalized surface potential decay.

    4.Degradation and performance improvement mechanism of SiC doped nanocomposite

    ReaxFF MD and DFT are important theoretical methods for the analysis of the microstructure evolution of the composite.The former can show the dynamic behavior of the physical and chemical reaction process and explain the performance change after being exposed to long-term corona, while the latter can provide theoretical explanation support for the electrical properties change.

    Since the solid insulation in the gas-insulated apparatus is immersed in SF6gas, the S/F/O reactive force field has always been crucial to studying the dynamic characteristics of its long-term corona behaviour using the ReaxFF MD method.Although S, F and O elements have appeared in the existing force fields, the connection between these three is empty and these force fields cannot be used to study the molecular dynamics behavior of epoxy composite/SF6systems.Based on the C/H/O/S/Si/F force field suitable for the SF6-EP-SiC system which could be obtained from our previous study through the proposed synergetic optimization methodology [29, 30], the model of the highly energetic particles from SF6decomposition bombarding the epoxy surface was constructed to simulate the corona erosion process in SF6gas, and it is given in figure 8(a).In this model,DGEBA and MTHPA are applied as the cross-linker,and the polymerization degree n of the epoxy resin monomer is set to 0, as shown in figures 8(b) and (c).The SiC nanoparticle model shown in figure 8(d) with diameter of 1.5 nm was constructed and optimized.The chemical molecular structures of the material involved are shown in figure 8(e).Besides the pure epoxy model, three SiC/epoxy composite models were also constructed,and their SiC contents were 3 wt.%,6 wt.%,and 9 wt.%,respectively.The pure epoxy model is denoted as Model P,and the other three models are denoted as Model S1,Model S2 and Model S3 in the above order.The crosslinked model was placed at the bottom of the simulation box with a z-axis length of 160 ?, and the z-direction is aperiodic.The most representative particles, SF and F, among the SF6decomposition products [15] are generated at a random position at a fixed height from the composite surface and shot vertically towards the polymer along the z-axis.A particle will inject every 2000 simulation steps from the beginning to the end of the simulation,and simulation steps are carried out with a micro-canonical ensemble (NVE).These energetic particles are all assigned energy similar to the actual situation,0.75 eV [8].An elastic wall was added at the bottom of the box to avoid ‘drift’ of the entire model.All ReaxFF MD calculations were performed with the software AMS.

    Figure 8.Molecular structure of (a) bombardment model, (b) DGEBA, (c) MTHPA, (d) SiC nanoparticle with diameter of 1.5 nm and(e) various types of atoms involved in the study.

    The meta-GGA calculation method with the 2-zeta basis group, M06-2X/6-31 G**level [31, 32], is employed to investigate and optimize the geometries of the simplified crosslinked epoxy resin and typical reaction products containing an S atom or F atom.The keyword,int=ultrafine,is added to improve the calculation accuracy.The crosslinked polymer molecule contains one DGEBA molecule and one MTHPA molecule, from which all the other reactants are derived.The DFT analysis is conducted with ORCA 4.2 software package [33] and the data processing is performed with Multiwfn [34] and VMD [35].

    4.1.Mechanism of change in physicochemical properties

    Taking Model P and Model S3 as examples, their ReaxFF MD simulation results were analyzed, and the states of the epoxy composite system after being impacted by different particles were intercepted, as shown in figure 9.It can be found that under the 50 ps continuous bombardment of F or SF,the surface molecules of pure epoxy become loose,while the molecular structure of the surface layer is still relatively tight after doping with SiC nanoparticles.Compared with the pure epoxy resin system,the gas side small molecule products from the SiC doped epoxy composite system are obviously less, indicating that the damage degree of the solid side molecular chain is smaller.At 100 ps, the decomposition degree of the solid side in the reaction system increased significantly compared with that at 50 ps,and more molecular fragments drifted into the gas side.However,the model doped with SiC is better than pure epoxy resin in both the number of small molecules and the degree of looseness.After the bombardment of the same time,the structural changes caused by SF particles are much more obvious than those caused by F particles.There are not only looser structures but also more small molecule products.

    It cannot be ignored that the chemical reaction process of incident particles leading to surface polymer chain erosion is basically the same regardless of whether or not SiC nanoparticles are doped.The SF particles are easily combined with the O atoms on the hydroxyl and ester groups in the epoxy resin molecule, thereby destroying the stability of the epoxy resin molecule.One of the most probable reaction pathways for S atoms to stay on the long chain of epoxy resins is shown in figure 10(a),and HF is also generated at the same time.The F particles will take away the H atoms on the C-H structure to generate HF, and at the same time, other F atoms will complement and form C-F bonds.The typical process of F atoms staying on the polymer surface is shown in figure 10(b).In addition, other intermolecular forces between nanoparticles and these highly energetic particles also provide conditions for the residence.It is through the physical and chemical reactions between the SF6decomposition products and the composite surface substances that the surface features are changed,as shown in figure 3.Not only will the density of the surface material change significantly, resulting in a delamination phenomenon,but also S and F elements will appear on the exposed sample surface.These representative physicochemical reaction processes will occur regardless of whether SiC is doped or not.

    The improvement mechanism of SiC nanoparticle doping on the corona resistance of epoxy resin is mainly manifested in the following aspects.First, SiC nanoparticles have strong physical and chemical stability, so the nanoparticles exposed to the surface of composites can effectively block the direct bombardment of epoxy resin molecules by incident highenergy particles without reacting with them.Then, nanoparticles can physically or chemically adsorb incident highenergy particles, thereby reducing the chance of reactive particles reacting with the crosslinked epoxy resin system.In addition,through covalent bonds or intermolecular forces,the nanoparticles are closely combined with the surrounding epoxy resin molecules, and the nanoparticles use their own stable structure to enhance the stability of the surrounding epoxy resin molecules,making the epoxy resin molecules less prone to random movement into the unprotected area,thereby reducing the possibility of epoxy resin molecules being attacked by high-energy particles generated by corona discharge.Figure 11 shows the atomic temperature distribution in model P and model S3 suffering the bombardment process at 80 ps.It can be seen that under the continuous bombardment of high-energy particles, the local atomic temperature can reach about 2000 K.A gradient distribution of temperature occurs along the z-axis across the solid side.When the composite is doped with nanoparticles, it will block and adsorb the incident high-energy particles.During this process,the kinetic energy of the incident particles is directly transferred to the nanoparticles themselves, so the temperature of the nanoparticles is higher than that of the surrounding molecule, as shown in figure 11(b).Considering that the thermal conductivity of SiC nanoparticles is tens of times higher than that of epoxy resin matrix,the energy obtained by nanoparticles from blocking and adsorbing incident particles is easier to disperse rapidly inside the particles and to the molecules around the nanoparticles.Localized overheating is thus avoided, thereby mitigating thermal decomposition of solid materials due to temperature buildup.Therefore, the corona erosion suffered by the composite is greatly alleviated.

    Table 1.Surface electric conductivity of pure epoxy and SiC/epoxy composite under the electric field of 1 kV mm-1.

    To fully grasp the influence of nanoparticles on the corona resistance of crosslinked epoxy, it is necessary to carry out a more in-depth quantitative analysis.These molecules with fewer than six carbon or sulfur atoms formed a compound and gas molecules are defined as the mass loss part, and the curve of the residual mass over time is depicted in figure 12.

    Under the bombardment of F particles, the normalized residual mass of model P and model S3 decreases with a slow trend, and the mass of them was about 0.975 at the end.The residual mass of these two models is slightly higher than the original value in the time range of 2-15 ps and 4-38 ps,respectively.The maximum value of the former is only 1.003,while the maximum value of the latter is only close to 1.005.This is closely related to the chemical reaction process in the corresponding time.At the beginning,the F particles are close to the surface,and the main process causing the mass increase is physical adsorption.As the simulation continues, the energy of the system gradually increases, and the reaction between the composite and F particles begins to accelerate,resulting in a decrease in the subsequent residual mass.It is worth noting that the final mass of model P and model S3 is very close, so it can be inferred that the erosion damage caused by F particles is relatively light, and the SiC nanoparticles in Model S3 have not even played a protective role.

    The residual mass of model P and model S3 bombarded by SF particles showed a significant trend of first increasing and then decreasing, and the final mass difference between them was obvious.The residual mass of Model P has a maximum value of 1.033, which appears around 51 ps; the maximum value of Model S3 is around 1.024 at 65 ps.Subsequently,the residual mass of both systems turns down.The decline rate of mass for the model P was significantly faster than that for the model S3.

    Figure 9.Representative snapshots of(a)model P and(b)model S3 under the impact with F particles, and snapshots of (c) model P and(d) model S3 under the impact with SF particles.From left to right,the simulation times are 50 ps and 100 ps in sequence.

    Figure 10.Typical reaction process of (a) SF particle and (b) F particle with cross-linked epoxy resin.

    Figure 11.Molecular structure of (a) model P and (b) model S3 under the impact of SF particles at 80 ps.Red balls represent atomic temperature above 2000 K and blue balls represent low temperature below 353.15 K.

    At the end of the simulation, the mass of model P is 0.912,and that of model S3 remains at only 0.987,which was even less than that of model S3 under the F particle’s bombardment.This is attributed to the relatively large relative molecular mass of S atoms in SF, and the initial reaction of SF particles on the composite surface is not very violent.Through physical adsorption, chemical adsorption and other processes, it becomes part of the composite, causing the overall quality of the two models to exceed the initial state.On the other hand,the SF particles have strong reactivity and undergo complex and violent chemical reactions with the crosslinked epoxy chain, resulting in a rapid decrease in the mass of model P.The presence of SiC nanoparticles plays an important protective role, so the subsequent mass loss of model S3 is lower.

    For experiments and tests, SF6is decomposed under corona discharge, and the epoxy resin sample is bombarded by both of F particles and SF particles, so F and S elements will be found in the EDX of figure 3.Since a large number of SF particles are bound to the composite surface by physical adsorption, they may then be separated from the SiC doped composites surface by other processes.Therefore,the S and F elements on the surface of SEC10 are less than those of pure epoxy resin.It is the effect of SF particles that leads to the situation in figure 6 that the mass of pure epoxy resin first increases and then decreases, and the mass of SiC doped samples increases with the increase of SiC content.

    To better clarify the damage degree of epoxy resin systems exposed to highly reactive particles,the average damage propagation depth (DPD) was calculated in this work, which was determined on the basis of a comparison between the normalized mass density profiles of the initial and final states of the systems, along the z-axis, parallel to the particle bombardment direction.The reference value for DPD is set to 0.90 in this research.The F-induced DPD change shown in figure 13 is an average calculated from the last 1000 steps of the simulation.DPD induced by other particles is shown in table 2.For the same composite model,the DPD value caused by SF particles is larger than that caused by F particles,which is inseparable from the higher chemical reactivity of SF particles and S atoms.Under the bombardment of the two kinds of particles, the DPD values of the model all showed a decreasing trend with the increase of the proportion of SiC.The minimum DPD values all appear in the model with SiC mass fraction of 9wt.%, which are 11.4 ? under the bombardment of S particles and 13.8 ? under the bombardment of SF particles.It is the penetration of these particles into the composite that causes the density change, which in turn creates the delamination phenomenon in figure 3.When the SiC content is increased, the DPD of energetic particles can be reduced and the delamination phenomenon can be alleviated.Then, model P and model S3 were processed to remove all gas products and only the molecular chains connected to the matrix were retained.The surface states of Model P and Model S3 were drawn as shown in figure 14.It can be seen that SiC nanoparticle effectively reduces the surface roughness of the composite after high energetic particle erosion and bombardment.The bombardment model constructed in this MD study effectively reflects the role of SiC nanoparticle in resisting corona erosion by comparing with the experimental measurement data given in figure 5(a).

    The simulation system data after 100 ps of bombardment by SF particles is processed, and all small molecule products containing C number less than 5 are removed.A python script is used to draw the surface microscopic morphology according to the atomic position distribution of solid substance,as shown in figure 14.The surface roughness values of model P,model S1,model S2 and model S3 were obtained by calculation, which are 5.43 ?, 4.55 ?, 3.75 ? and 3.42 ?,respectively.As the content of SiC nanoparticles increases,the surface roughness after bombardment decreases gradually,which is the same trend as the experimental values after 30 h exposure in figure 5(a).

    Figure 12.Normalized residual mass of nanocomposite model.

    Figure 13.Normalized density along the z-axis under impact with F.

    Table 2.DPD values of different composite models after 100 ps.

    4.2.Deterioration mechanism of electric properties

    It can be seen from the measured value of surface roughness in figure 5(a) that the surface roughness of the composite increases with increasing corona erosion time.The current research on the influence of roughness on the insulation performance of epoxy resin composites shows that, in the range of micron scale, a proper increase in roughness could improve the DC flashover voltage [36-38].The value of composite surface roughness after corona erosion should meet the applicable range of this law.However, the flashover voltage along the surface decreases rapidly with the increase of corona treatment time in figure 6.This indicates that the roughness is not the key to influence the composite along surface insulation performance after corona erosion.On the other hand,combining the experimental tests in figure 3 with the reaction mechanism analysis in figure 10, it can be speculated that the change of the elemental composition in the corona-exposed and eroded region is an important reason for the surface insulation failure.

    Since the HOMO-LUMO gap of the oligomer is closely related to the band gap of an infinitely extended polymer,the band gap can be extrapolated from the HOMO-LUMO gap calculated by a quantum chemical program.The molecule models,DGEBA/MTHPA,F-DGEBA/F-MTHPA,S-DGEBA/S-MTHPA, shown in figure 15, are respectively constructed based on the typical reaction products of F particle and SF particle with crosslinked epoxy resin molecules in figure 10.

    The orbital cloud map and energy level orbital submap of the above three single-chain molecules are obtained by calculation, as shown in figure 16.The cloud images on the top and bottom are LUMO and HOMO, respectively.

    The LUMO value of F-DGEBA/F-MTHPA is not affected by a small amount of F element substitution, and its value is almost the same as that of DGEBA/MTHPA.There is a slight difference in the HOMO values of the two,and the value of F-DGEBA/F-MTHPA is slightly smaller than that of DGEBA/MTHPA.The energy gaps of DGEBA/MTHPA and F-DGEBA/F-MTHPA are 7.785 eV and 7.844 eV,respectively,and the latter is 0.059 eV larger than the former.The increase of the band gap indicates the effects of F element substitution, elevating the energy required by the electrons and raising the barrier of electron migration between different lattices.Meanwhile,for DGEBA/MTHPA,there are obvious electron and hole traps.The F element makes the electron trap of F-DGEBA/F-MTHPA deeper and the hole trap slightly shallower.It is shown that the F element can change the effect of the potential decay process of the polymer and improve the insulation performance to a certain extent.

    For S-DGEBA/S-MTHPA, the S element significantly reduces the LUMO and slightly increases the HOMO, which makes the gap drop to only 1.875 eV, which is no longer an insulating medium.It can be concluded that after the excess SF particle reacts with the crosslinked polymer,the S element resides on the composite surface exposed to the corona,and it is the key factor causing the surface insulation failure and significant increase in conductivity.The longer the exposure time to corona, the more S element is accumulated on the composite surface, which will degrade the surface insulation more seriously.In addition, for the initial state and the samples exposed to corona for different amounts of time, the flashover voltage along the surface has a maximum value when the mass fraction of SiC is 6 wt.%, which is consistent with the nonlinear conductivity characteristics of SiC and its regulation of deep and shallow traps in composites.The SiC nanoparticles maintain the electrical performance along the surface by reducing the accumulation of S elements on the surface exposed to long-term corona.

    Figure 14.Surface topography after 100 ps bombardment by SF particles.

    Figure 15.Optimized representative molecular models:(a)DGEBA/MTHPA, (b) F-DGEBA/F-MTHPA, (c) S-DGEBA/S-MTHPA.

    5.Conclusions

    In this study, through experimental research, the long-term service performance of SiC-modified nanocomposite materials under the extreme condition of corona erosion has been established.Moreover, ReaxFF MD and DFT were used to reveal the deterioration mechanism of the composite material by corona erosion and the mechanism of corona erosion resistance of nanocomposite enhanced by SiC nanoparticles from the atomic scale.The detailed conclusions are as follows:

    (1) The variation law of key physicochemical properties and insulation strength after long-term service in corona discharge has been established.A multi-needle-plate electrode platform was used to carry out long-term corona erosion experiments on the SiC-doped nanocomposite.With increasing exposure time, various key properties and indicators were deteriorated, and S and F elements were found on the exposed surface.The erosion depth and surface roughness of the sample doped with 10 wt.% SiC remained optimal.Regarding insulating properties, the DC flashover voltages of composites doped with 6 wt.% and 8 wt.% SiC are 6.3 kV and 7.9 kV after 30 h of erosion,respectively, while the others are close to losing, or have completely lost, their insulating properties.

    Figure 16.Energy level distribution and molecular orbitals of epoxy resin chain after being exposed to bombardment.

    (2) The degradation mechanism of the physicochemical properties of composite exposed to corona discharge was elucidated using ReaxFF MD simulation methodology.Based on the developed reactive force field containing C, H, O, S, Si, and F elements, a model suitable for simulating the corona discharge erosion of solid insulation by high-energy particles bombarding the surface is constructed.Typical chemical reactions of F particles and SF particles with cross-linked epoxy resins leading to the decomposition of composite were found, which explained the change in the elemental composition obtained by EDX.The corona resistance of the nanocomposite was improved by SiC enhancing the interaction with the adjacent epoxy resin long chain and exposing itself to the direct bombardment of highenergy particles.The data on the dynamic damage characteristics of the interface is analyzed by the python script developed.The normalized residual mass and damage propagation depth (DPD) are defined to evaluate the interface features, and a 3D surface topography of the solid part is constructed to compare with the macroscopic experimental results.It can be confirmed that SF particles are more reactive and more destructive, and this is the primary reason for the changes in various physicochemical properties of the nanocomposite surface.

    (3) The deterioration mechanism of the insulation properties for the SiC-doped nanocomposite was elucidated with DFT from the atomic-scale perspective.Energy level orbital analysis was performed on typical molecular structures obtained from ReaxFF MD analysis.The band gap of molecules containing F was raised and the insulating properties were enhanced.For the molecule containing S,the band gap drops directly from the initial 7.785 eV to 1.875 eV.Thus,the continuous accumulation of S-containing derivatives on the composite surface after corona erosion leads to the decline and failure of the surface insulation.

    (4) Comparing the simulation results to experimental results shows that the corona erosion model based on microscale simulation effectively reflects some of the characteristics in the experiment.The experimental and theoretical analysis methods proposed in this study have important industrial applications and reference value in designing, developing and evaluating the performance of novel corona resistant insulating materials in gas compressed equipment.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China (Nos.51737005, 51929701, 52177147 and 52127812).

    猜你喜歡
    王健
    王健
    當代作家(2023年5期)2023-07-10 22:13:48
    美聯(lián)儲的艱難選擇:穩(wěn)通脹還是穩(wěn)金融市場
    鉗工機械操作的質(zhì)量控制策略
    Exact solution of an integrable quantum spin chain with competing interactions?
    疫情期間,甲狀腺結(jié)節(jié)患者該怎么辦
    抗癌之窗(2020年1期)2020-05-16 01:37:26
    每天只練五分鐘
    “王健扇藝展”
    你像春神
    初冬
    “平移”檢測題
    亚洲成人久久爱视频| 精品久久久久久久人妻蜜臀av| 欧美成人免费av一区二区三区| 晚上一个人看的免费电影| 女的被弄到高潮叫床怎么办| 日本在线视频免费播放| 久久久a久久爽久久v久久| 成人综合一区亚洲| 看黄色毛片网站| 欧美变态另类bdsm刘玥| 麻豆成人av视频| 特级一级黄色大片| 白带黄色成豆腐渣| 亚洲电影在线观看av| 国产极品精品免费视频能看的| 亚洲中文字幕一区二区三区有码在线看| 91aial.com中文字幕在线观看| 男女边吃奶边做爰视频| 亚洲欧美中文字幕日韩二区| 一级毛片aaaaaa免费看小| 美女脱内裤让男人舔精品视频 | 一边亲一边摸免费视频| av卡一久久| 国产成人a∨麻豆精品| 精品无人区乱码1区二区| 亚洲av成人精品一区久久| 一本精品99久久精品77| 欧美激情国产日韩精品一区| 国产精品一区二区性色av| 日韩一本色道免费dvd| 97超视频在线观看视频| 一本一本综合久久| 男女那种视频在线观看| 国产亚洲精品久久久久久毛片| 久久精品夜夜夜夜夜久久蜜豆| 美女xxoo啪啪120秒动态图| 午夜爱爱视频在线播放| 天堂影院成人在线观看| 全区人妻精品视频| 丰满的人妻完整版| 国内精品宾馆在线| 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 国产亚洲5aaaaa淫片| 国产一区二区在线av高清观看| 午夜福利视频1000在线观看| 别揉我奶头 嗯啊视频| 久久精品国产亚洲av涩爱 | 女的被弄到高潮叫床怎么办| 国内精品一区二区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 91在线精品国自产拍蜜月| 色视频www国产| 午夜视频国产福利| 午夜福利在线在线| 日本免费a在线| 久久久久性生活片| 亚洲在久久综合| 日韩精品青青久久久久久| 中文欧美无线码| 精品国内亚洲2022精品成人| 婷婷色综合大香蕉| 中文欧美无线码| 哪个播放器可以免费观看大片| 国产午夜精品久久久久久一区二区三区| 麻豆久久精品国产亚洲av| 国产精品综合久久久久久久免费| 免费看a级黄色片| 午夜a级毛片| 国产精品人妻久久久影院| 舔av片在线| 免费搜索国产男女视频| 麻豆成人午夜福利视频| 久久精品久久久久久久性| 色尼玛亚洲综合影院| 成人美女网站在线观看视频| 97超碰精品成人国产| 国产精品.久久久| 一级毛片aaaaaa免费看小| 91久久精品国产一区二区成人| 波多野结衣高清作品| 国产成人精品婷婷| 免费看a级黄色片| 又爽又黄无遮挡网站| 国产精品美女特级片免费视频播放器| 亚洲婷婷狠狠爱综合网| 亚洲欧美中文字幕日韩二区| 亚洲成a人片在线一区二区| 午夜福利在线观看吧| 在线观看午夜福利视频| 久久午夜福利片| 99久久无色码亚洲精品果冻| 亚洲五月天丁香| www.色视频.com| 高清日韩中文字幕在线| 亚洲熟妇中文字幕五十中出| 亚洲av免费在线观看| 中国国产av一级| 日韩av在线大香蕉| 国产片特级美女逼逼视频| 一级av片app| 成人一区二区视频在线观看| 欧美激情国产日韩精品一区| 国产极品精品免费视频能看的| 村上凉子中文字幕在线| 久久久成人免费电影| 亚洲最大成人av| 久久午夜福利片| 成年免费大片在线观看| av在线老鸭窝| h日本视频在线播放| 久久这里只有精品中国| 一本久久中文字幕| 欧美高清成人免费视频www| 亚洲精品乱码久久久久久按摩| 综合色av麻豆| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 午夜老司机福利剧场| 别揉我奶头 嗯啊视频| 日本成人三级电影网站| 五月玫瑰六月丁香| 久久精品国产亚洲网站| 国语自产精品视频在线第100页| 国产精品av视频在线免费观看| 欧美激情在线99| 国产91av在线免费观看| 亚洲人成网站在线播放欧美日韩| 97在线视频观看| 青春草国产在线视频 | 久久热精品热| 超碰av人人做人人爽久久| 晚上一个人看的免费电影| kizo精华| 12—13女人毛片做爰片一| 男女边吃奶边做爰视频| 中文字幕免费在线视频6| 一级黄色大片毛片| 91麻豆精品激情在线观看国产| 国产亚洲av嫩草精品影院| 校园春色视频在线观看| 悠悠久久av| 九九热线精品视视频播放| 国产中年淑女户外野战色| 尾随美女入室| 日日撸夜夜添| 免费在线观看成人毛片| 夫妻性生交免费视频一级片| 男女边吃奶边做爰视频| 免费看日本二区| 国产精品一区二区三区四区免费观看| 国产精品久久久久久精品电影| 给我免费播放毛片高清在线观看| 91精品一卡2卡3卡4卡| 少妇人妻一区二区三区视频| 午夜精品在线福利| 精品少妇黑人巨大在线播放 | 两个人视频免费观看高清| 亚洲国产精品国产精品| 亚洲激情五月婷婷啪啪| 99久久无色码亚洲精品果冻| 国产在线男女| 热99在线观看视频| 成人av在线播放网站| 国产午夜精品论理片| 欧美bdsm另类| 精品人妻熟女av久视频| 欧美激情久久久久久爽电影| 亚洲av男天堂| 国产亚洲欧美98| 免费av不卡在线播放| 国产av不卡久久| 国产亚洲欧美98| 一个人看的www免费观看视频| av又黄又爽大尺度在线免费看 | 免费看日本二区| 美女脱内裤让男人舔精品视频 | 十八禁国产超污无遮挡网站| 91午夜精品亚洲一区二区三区| 一级毛片久久久久久久久女| 又黄又爽又刺激的免费视频.| 男女啪啪激烈高潮av片| av卡一久久| 亚洲三级黄色毛片| 精品久久久久久久久亚洲| 午夜老司机福利剧场| 成人午夜精彩视频在线观看| 精品一区二区免费观看| 亚洲一区高清亚洲精品| 国产一级毛片七仙女欲春2| or卡值多少钱| 国产在线精品亚洲第一网站| av国产免费在线观看| 成人无遮挡网站| 亚洲丝袜综合中文字幕| 草草在线视频免费看| 中国美女看黄片| 亚洲丝袜综合中文字幕| 精品久久久久久久久亚洲| 成人无遮挡网站| 亚洲人成网站在线观看播放| 国产精品,欧美在线| 日韩 亚洲 欧美在线| 日本熟妇午夜| 长腿黑丝高跟| 夜夜夜夜夜久久久久| 亚洲人成网站在线播放欧美日韩| 国产黄片视频在线免费观看| 国产精品野战在线观看| 村上凉子中文字幕在线| 麻豆国产av国片精品| 不卡视频在线观看欧美| 一区福利在线观看| 国产人妻一区二区三区在| 久久精品国产亚洲av涩爱 | 天堂av国产一区二区熟女人妻| 精品人妻视频免费看| 国产精品久久视频播放| 日韩一区二区视频免费看| 自拍偷自拍亚洲精品老妇| av卡一久久| 波多野结衣巨乳人妻| 少妇被粗大猛烈的视频| 国产视频内射| 插逼视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成人精品一,二区 | 亚洲人成网站高清观看| 久久久精品大字幕| 97人妻精品一区二区三区麻豆| 性色avwww在线观看| 尾随美女入室| 国产成年人精品一区二区| 青春草国产在线视频 | 国产成人福利小说| 亚洲美女搞黄在线观看| 精品人妻视频免费看| 精品欧美国产一区二区三| 伦精品一区二区三区| 久久中文看片网| 中文字幕精品亚洲无线码一区| 少妇熟女aⅴ在线视频| av又黄又爽大尺度在线免费看 | 99久久成人亚洲精品观看| 2022亚洲国产成人精品| 亚洲欧美中文字幕日韩二区| 夜夜看夜夜爽夜夜摸| 天美传媒精品一区二区| 午夜视频国产福利| 又粗又硬又长又爽又黄的视频 | 亚洲真实伦在线观看| 国产精品久久视频播放| 国内精品一区二区在线观看| 国产精品日韩av在线免费观看| 亚洲不卡免费看| 久久久久久九九精品二区国产| 久久6这里有精品| 欧美一级a爱片免费观看看| 国产亚洲5aaaaa淫片| 给我免费播放毛片高清在线观看| 网址你懂的国产日韩在线| 九九热线精品视视频播放| 国产大屁股一区二区在线视频| 亚洲激情五月婷婷啪啪| 中国美女看黄片| 亚洲性久久影院| 免费不卡的大黄色大毛片视频在线观看 | 你懂的网址亚洲精品在线观看 | 晚上一个人看的免费电影| ponron亚洲| 一个人观看的视频www高清免费观看| 看黄色毛片网站| 日韩三级伦理在线观看| 精品99又大又爽又粗少妇毛片| 中文字幕av成人在线电影| 成人无遮挡网站| 日韩国内少妇激情av| 久久久久久久午夜电影| 免费人成在线观看视频色| 国产蜜桃级精品一区二区三区| 啦啦啦啦在线视频资源| 日本五十路高清| 麻豆国产av国片精品| 国产高清不卡午夜福利| 毛片一级片免费看久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 成人特级av手机在线观看| 男女那种视频在线观看| 少妇丰满av| 免费在线观看成人毛片| 又粗又爽又猛毛片免费看| 男人的好看免费观看在线视频| 国内精品一区二区在线观看| 校园人妻丝袜中文字幕| 亚洲成av人片在线播放无| 亚洲av免费在线观看| 亚洲人成网站在线播| 亚洲精品久久久久久婷婷小说 | 乱系列少妇在线播放| 亚洲成人av在线免费| 欧美性感艳星| 中文资源天堂在线| 舔av片在线| 国产精品电影一区二区三区| 黑人高潮一二区| 一卡2卡三卡四卡精品乱码亚洲| 日日啪夜夜撸| 少妇人妻精品综合一区二区 | 18禁在线无遮挡免费观看视频| 国产免费男女视频| 黄色视频,在线免费观看| 成人一区二区视频在线观看| 久久精品影院6| 国产乱人偷精品视频| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| 黄片无遮挡物在线观看| 日本撒尿小便嘘嘘汇集6| av在线蜜桃| 中国美女看黄片| 大香蕉久久网| 一级二级三级毛片免费看| 啦啦啦韩国在线观看视频| 91麻豆精品激情在线观看国产| 成年av动漫网址| 只有这里有精品99| 秋霞在线观看毛片| 日韩一本色道免费dvd| 亚洲欧美日韩高清专用| 亚洲欧美精品自产自拍| 亚洲最大成人手机在线| 亚洲欧美精品自产自拍| 日韩欧美精品免费久久| 九草在线视频观看| 日日干狠狠操夜夜爽| 国产伦在线观看视频一区| 国国产精品蜜臀av免费| 午夜视频国产福利| 欧美性猛交黑人性爽| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| a级一级毛片免费在线观看| 国产男人的电影天堂91| 国产成人精品一,二区 | 国产探花极品一区二区| 亚洲av一区综合| 日韩成人伦理影院| 黄色视频,在线免费观看| 赤兔流量卡办理| 欧美精品国产亚洲| 欧美一级a爱片免费观看看| 在线免费观看的www视频| 色尼玛亚洲综合影院| 国内揄拍国产精品人妻在线| 精品久久久久久成人av| 一区福利在线观看| 亚洲国产精品国产精品| 丰满的人妻完整版| 99久久人妻综合| 久久这里只有精品中国| 精品久久久久久久久久久久久| 国产不卡一卡二| 晚上一个人看的免费电影| 亚洲国产精品sss在线观看| 欧美丝袜亚洲另类| 亚洲精品乱码久久久久久按摩| 成人性生交大片免费视频hd| 欧美+亚洲+日韩+国产| 一进一出抽搐动态| 亚洲精品久久久久久婷婷小说 | 中文字幕制服av| av免费观看日本| 91精品国产九色| 免费观看a级毛片全部| 桃色一区二区三区在线观看| 亚洲欧美精品专区久久| 亚洲一区高清亚洲精品| 欧美日韩精品成人综合77777| 久久久国产成人精品二区| 蜜臀久久99精品久久宅男| 人妻系列 视频| 波多野结衣高清无吗| 少妇人妻精品综合一区二区 | 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 一进一出抽搐gif免费好疼| 好男人在线观看高清免费视频| 麻豆久久精品国产亚洲av| 99在线视频只有这里精品首页| 成人毛片60女人毛片免费| 久久久色成人| 91午夜精品亚洲一区二区三区| 国产美女午夜福利| 日日啪夜夜撸| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一及| 免费观看精品视频网站| 一级毛片电影观看 | 少妇高潮的动态图| 日韩欧美精品免费久久| 精品一区二区三区人妻视频| 亚洲人与动物交配视频| 中文欧美无线码| 嘟嘟电影网在线观看| 少妇人妻一区二区三区视频| 亚洲久久久久久中文字幕| 久久婷婷人人爽人人干人人爱| 亚洲国产精品久久男人天堂| 男女下面进入的视频免费午夜| 中文字幕制服av| 九九久久精品国产亚洲av麻豆| av黄色大香蕉| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 久久久精品欧美日韩精品| kizo精华| 岛国在线免费视频观看| 国产老妇女一区| 午夜免费激情av| 色哟哟哟哟哟哟| 国产黄色视频一区二区在线观看 | 欧美性猛交黑人性爽| 日韩成人伦理影院| 欧美zozozo另类| 白带黄色成豆腐渣| 美女cb高潮喷水在线观看| 天堂中文最新版在线下载 | 变态另类成人亚洲欧美熟女| 99久国产av精品国产电影| 91久久精品电影网| 99热6这里只有精品| 国产精品无大码| 如何舔出高潮| 久久99热这里只有精品18| 级片在线观看| 色哟哟·www| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清| 亚洲国产日韩欧美精品在线观看| av视频在线观看入口| 夫妻性生交免费视频一级片| 老司机影院成人| 亚洲精品亚洲一区二区| 国产精品,欧美在线| 97在线视频观看| 99热这里只有是精品50| 久久精品国产清高在天天线| 亚洲久久久久久中文字幕| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 日本免费一区二区三区高清不卡| 麻豆国产av国片精品| 日韩欧美 国产精品| 成人性生交大片免费视频hd| 一进一出抽搐gif免费好疼| 日韩欧美精品免费久久| 免费av不卡在线播放| 99久久无色码亚洲精品果冻| 网址你懂的国产日韩在线| 91久久精品国产一区二区成人| 国产成人freesex在线| 亚洲欧洲日产国产| 久久这里有精品视频免费| 国产一区二区在线av高清观看| 国产成人a区在线观看| 免费av毛片视频| 亚洲av不卡在线观看| 久久热精品热| 在现免费观看毛片| 熟女电影av网| 六月丁香七月| 免费看光身美女| 超碰av人人做人人爽久久| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| 又黄又爽又刺激的免费视频.| 国产成人91sexporn| 精品日产1卡2卡| 特级一级黄色大片| av国产免费在线观看| 热99re8久久精品国产| 日韩一区二区视频免费看| 日本黄色视频三级网站网址| 99久久精品国产国产毛片| 亚洲成人av在线免费| 日韩高清综合在线| av天堂中文字幕网| 亚洲五月天丁香| 亚洲欧美日韩东京热| or卡值多少钱| 国产色婷婷99| 免费电影在线观看免费观看| 一级av片app| 成人毛片a级毛片在线播放| 午夜老司机福利剧场| 久久鲁丝午夜福利片| 九九热线精品视视频播放| 欧美潮喷喷水| 亚洲国产欧美在线一区| 国产精品久久电影中文字幕| 国产极品精品免费视频能看的| 日韩av在线大香蕉| 国产伦理片在线播放av一区 | 亚洲精品粉嫩美女一区| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄 | 一区二区三区免费毛片| 国产精品电影一区二区三区| 精品少妇黑人巨大在线播放 | 全区人妻精品视频| 嫩草影院入口| 欧美不卡视频在线免费观看| 激情 狠狠 欧美| 在线国产一区二区在线| 国国产精品蜜臀av免费| 成人一区二区视频在线观看| www日本黄色视频网| 免费不卡的大黄色大毛片视频在线观看 | 成人一区二区视频在线观看| 久久人人爽人人爽人人片va| 亚洲第一电影网av| 色综合亚洲欧美另类图片| 在线播放国产精品三级| 精品一区二区免费观看| 在线播放国产精品三级| 亚洲av男天堂| 久久久精品大字幕| 国产av不卡久久| 看黄色毛片网站| 精品免费久久久久久久清纯| 国产毛片a区久久久久| 一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 久久99精品国语久久久| 国产人妻一区二区三区在| 1000部很黄的大片| 99久国产av精品国产电影| 99久久九九国产精品国产免费| 最近中文字幕高清免费大全6| 亚洲经典国产精华液单| 99riav亚洲国产免费| 天堂中文最新版在线下载 | 欧美激情在线99| 内射极品少妇av片p| 精品久久国产蜜桃| 亚洲av中文字字幕乱码综合| 有码 亚洲区| 国产成人aa在线观看| 在线天堂最新版资源| 午夜老司机福利剧场| 成人特级黄色片久久久久久久| 在线观看午夜福利视频| 国产极品精品免费视频能看的| 亚洲欧美中文字幕日韩二区| 国内少妇人妻偷人精品xxx网站| 久久99热6这里只有精品| 亚洲欧美成人综合另类久久久 | 日韩欧美国产在线观看| 九色成人免费人妻av| 人人妻人人看人人澡| 18禁在线播放成人免费| 亚洲在久久综合| 综合色丁香网| 看免费成人av毛片| 一个人观看的视频www高清免费观看| 欧美不卡视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 日本黄色视频三级网站网址| 久久久国产成人免费| 亚洲精品日韩在线中文字幕 | 自拍偷自拍亚洲精品老妇| 亚洲一区高清亚洲精品| eeuss影院久久| 日本在线视频免费播放| 亚洲av免费在线观看| 免费一级毛片在线播放高清视频| 综合色av麻豆| 亚洲欧美日韩东京热| 97热精品久久久久久| 色哟哟哟哟哟哟| 精品无人区乱码1区二区| 97人妻精品一区二区三区麻豆| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 午夜精品一区二区三区免费看| 欧美一区二区国产精品久久精品| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 天堂中文最新版在线下载 | 亚洲精品456在线播放app| 26uuu在线亚洲综合色| 国产色爽女视频免费观看| 亚洲av第一区精品v没综合| 成人毛片a级毛片在线播放| 亚洲第一电影网av| 在线观看美女被高潮喷水网站| 亚洲成人久久性| 国产精品永久免费网站| av在线天堂中文字幕| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av天美| 偷拍熟女少妇极品色| 一区二区三区高清视频在线| 亚洲国产精品国产精品| 日本免费a在线| 中文字幕av成人在线电影| 国产又黄又爽又无遮挡在线| 两个人视频免费观看高清| 亚洲精品日韩在线中文字幕 | 欧美性感艳星| 国产高潮美女av|