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

    Experimental study of hydrodynamic performance of full-scale horizontal axis tidal current turbine*

    2017-03-09 09:09:42FengmeiJing荊豐梅WeijiaMa馬偉佳LiangZhang張亮ShuqiWang王樹齊XiaohangWang王曉航
    關(guān)鍵詞:張亮

    Feng-mei Jing (荊豐梅), Wei-jia Ma (馬偉佳), Liang Zhang (張亮), Shu-qi Wang (王樹齊), Xiao-hang Wang (王曉航)

    1. College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China, E-mail: jingfengmei@hrbeu.edu.cn

    2. School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China

    3. Harbin Electric Machinery Company Limited, Harbin 150040, China

    (Received February 15, 2015, Revised November 11, 2015)

    Experimental study of hydrodynamic performance of full-scale horizontal axis tidal current turbine*

    Feng-mei Jing (荊豐梅)1, Wei-jia Ma (馬偉佳)1, Liang Zhang (張亮)1, Shu-qi Wang (王樹齊)2, Xiao-hang Wang (王曉航)3

    1. College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China, E-mail: jingfengmei@hrbeu.edu.cn

    2. School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China

    3. Harbin Electric Machinery Company Limited, Harbin 150040, China

    (Received February 15, 2015, Revised November 11, 2015)

    In this paper, experiments of both the model turbine (1 kW) and the full scale (10 kW) turbine are carried out in a towing tank and a basin, respectively, and the test of the full scale turbine on the sea is conducted. By comparison between the model turbine (D=0.7 m)and the full scale turbine (D=2.0 m), it is shown that the maximum power coefficient increases with the increase of the diameter of the turbine. The test results on the sea are used to study the hydrodynamic performances of the horizontal axis turbine, and provide a basis for the design. Experimental results can validate the accuracy of the numerical simulation results.

    Tidal current energy, horizontal axis turbine, model test, full scale test, sea trial

    Introduction

    The oceans cover more than 70% of the earth’s surface and offer huge energy resources for the production of electricity and/or fresh water[1-3]. The ocean energy is available in various forms, including tidal wave, tidal current, thermal energy, salinity gradients and biomass[4]. The tidal current energy technology shares many similar fundamental characteristics with the wind energy. Like a wind turbine, the tidal current turbine converts the kinetic energy of the fluid into the electrical energy. However, due to the difference ofthe working environment, there are some distinct phenomena of the tidal current turbine from the wind turbine, including cavitations, free surface, and waveinduced velocity gradient, etc.[5]. For example, compared with the wind turbine, the tidal current turbine has an advantage that those reversing gravity loads, such as the buoyancy, can be used to balance the weight, which is very important for the horizontal-axis tidal current turbine. Given these considerations, it is necessary to study the characteristics of the tidal current turbine, which is of interest for designers.

    The investigation of the tidal current energy is a classic topic in recent decades and many fundamental issues were studied by both numerical and experimental approaches. In the design of the horizontal axis turbine, many parameters should be considered, including loads, depth and cavitation occurrence, so the experiment is an effective means to explore the characteristics of the horizontal axis turbine. This paper presents only the experiment work on the horizontalaxis tidal current turbines, and the experiment work on the vertical-axis tidal current turbines were described in the literature[6,7].

    A number of experiments were conducted to cha-racterize the turbine performance. Bahaj et al.[8,9]conducted experiments under conditions of different airfoils, pitches and yaws, blade cavitations and twinrotor interactions, to validate blade element momentum models. Wang et al.[10]conducted experimental and numerical study under different rotor preset angles. Maganga et al.[11]conducted experiments in a closedloop water tunnel, examining the effects of the inflow turbulence intensity, varying from 8% to 25%, the yaw angle ranging from ?10oto 20o, and the depth beneath the free surface, including observations at 0.94D, 1.37D and 2.04Din a channel of 3.14Ddeep. Galloway et al.[12]carried out small-scale mesh disk rotor simulators to study the structure of the wake, specifically with regard to the depth. Tests were conducted in a flume of 0.3 m deep with disk of 0.1 m in diameter. The depth, measured from the center of the disk, was set to 0.33D,0.50D,0.66Dand0.75Dand the velocity deficit was examined from a distance of 20D. The velocity deficit was found to be higher for shallow settings, and at the deepest settings the disk was considered as a submerged obstruction. Barltrop et al.[13]conducted experiments in a towing-tank using a three-bladed turbine with the diameter of 0.4 m and featuring an S814 airfoil cross-section.

    To the authors’ knowledge, there are few experimental studies of the full scale turbine, which serves as a challenge. The full scale turbine experiment is very expensive while the sites that can be used to test the turbine performance are hard to find; because the loads on a full scale turbine are larger than those on a model turbine, the requirements for the test equipment are very severe, which imposes another difficulty for the full scale turbine experiment. However, the full scale turbine test data are important for the study of the horizontal axis tidal current turbine. In this paper, the model turbine experiment is conducted in a towing tank and the full scale turbine is examined in a deep basin, moreover, the full scale turbine is tested in the real sea. The relationship between the model turbine and the full scale turbine is studied by comparing the power coefficient. Experimental results could be used as the basis for the design of the horizontal axis turbine, and also could be used to validate the accuracy of numerical simulation results.

    As we all know, in an experiment, if the Froude number is made to satisfy the requirements, it is very difficult to make the Reynolds numbers also satisfy the requirements[14,15], so the scale effect is unavoidable in the experiment. In order to reduce the experiment error, the larger the model, the better the test result will be. However, the large model will increase the cost and make the requirements for the test equipment more complex. The large model is better to be tested in a large basin where the blockage effect can be reduced. Therefore, this paper makes a comparison of results between the 1 kW model turbine (D= 0.7 m)and the 10 kW full scale turbine (D=2.0 m), which would be significantly important for analysis of the scale effect on the horizontal axis tidal current turbine.

    1. Turbine experiments

    The turbine used in the experiment is a twobladed horizontal axis turbine, where the diameter of the model turbine is 0.7 m and the diameter of the full scale turbine is 2.0 m. The model experiment is carried out in the towing tank of Harbin Engineering University (HEU), the full scale turbine test is carried out in a basin of HEU, and the test of the full scale turbine on the sea is carried out in the Daishan Island of China.

    Table 1 Particulars of turbine blades

    1.1Description of blades

    To improve the utilization efficiency of tidal energy, the blade element momentum (BEM) based Wilson model is applied to the design of the horizontal tidal turbine[16]. The blades of the turbine are developed by Institute of Ocean Renewable Energy System of HEU, based on the simplified wind theory and the Glauert vortex theory[17,18].The blades aredeveloped based on the profile S809, and blade radius (R), the chord (c), pitch distribution (?)and thickness(σ)are shown in Table 1. The blades of the model turbine are made of copper(D=0.7 m)while the blades of the full scale turbine are made of stainless steel(D=2.0 m).

    1.2Turbine experiments

    1.2.1 Experimental methods

    The output power of the tidal current turbine gives electric powers, which include the bearing and seal losses[19,20]. The electric power (Pd)is expressed as

    whereu(t)is the instantaneous voltage generated by generator,i(t)is the instantaneous current generated by generator,tis the cumulative sampling time.

    The signals of the voltage and the current of the electric power are discretized.

    where Δtis the interval sampling time,Nis the sampling frequency,ukis the transient voltage value sampled,ikis the transient current value sampled.

    1.2.2 Model turbine equipment and test rig

    Experiments of the model turbine of diameter of 0.7 m are carried out in the towing tank at Harbin Engineering University. It is 108.0 m long, 7.0 m wide and 3.5 m deep with a maximum towing speed of 6 m/s. The tank is equipped with a piston type wave maker at the upstream end, and can produce regular and irregular waves with the maximum amplitude of 0.4 m and a period range of 0.4 s-14.0 s. At the downstream end there is a wave absorbing beach.

    In the test, the DHDAS data acquisition system is adopted. This system can accurately collect signals of test data, and can ensure the signals synchronized in time. There are two channels in the DHDAS system to display the voltage and the current, respectively.

    The output signals of the generator are AC signals, which are transformed to DC signals via rectifiers and are considered as the input of the control system. The electrical power is absorbed with rheostats via a rotary encoder mounted behind the generator, which also allows the regulation of the rotor speed. The power signals are transferred to PC via the data collector and then the voltage and the current can be read from PC. The test diagram is shown in Fig.1. Figure 2 is the chart of the experiment process.

    Fig.1 Schematic diagram of experiment test

    Fig.2 Chart of experiment process

    In order to test the model turbine, an experimental test rig that could be bolted to the carriage of the towing tank is designed. The horizontal axis turbine is fixed to an elliptic column, which is connected to a circular column by a flange, and the circular column is fixed to the girder of the carrier. The test rig is shown in Fig.3 .The axis of the turbine is propelled to the flow direction, the rotor is attached to a main shaft, When the rotor rotates, a DC generator is driven through the main shaft, thus the kinetic energy of the fluid is converted to the electric energy. The model copper blades are milled to their shapes and in order to improve sealing, the nacelle is cast from aluminum alloy, with a hole inside to house the generator. For the model turbine (D=0.7 m), the rated current speed is 2.0 m/s and the rated rotation speed of the generator is 350 rpm, and the rated output power of the generator is 1 kW. The objectives of the model turbine experiment include:

    (1) Measuring the output voltage and current to obtain the power performance of the model turbine, when the flow speed is 1.0 m/s, 1.5 m/s and 2.0 m/s, and the installation pitch angle is0o,?2.5o,2.5oand 5o.

    Fig.3 (Color online) Photos of experiment model (1 kW)

    (2) Measuring the output voltage and current to obtain the power performance of the model turbine, when the turbine is installed in the opposite direction, that is, the flow goes from the cabin to the blades.

    (3) To obtain the power characteristics of the model turbine, while the turbine is working under wave and current.

    In this paper, only the first point is discussed, and the others would be discussed in the next paper.

    1.3Full scale turbine experiment and test rig

    The experiment of the full scale turbine with diameter of 2 m is conducted in the basin of HEU. The basin is 50 m long, 30 m wide and 10 m deep. The basin is also equipped with a piston type wave maker. The rated current speed of the full scale turbine is 2.0 m/s, and the rated rotation speed of the generator is 90 rpm/s and the output power of the generator is 10 kW. And the parameters of the full scale turbine are shown in Table 2. The objectives of the full scale turbine experiment are:

    (1) Based on the test of the hydrodynamic performance, to obtain the power characteristics of the turbine.

    (2) Based on the test of the hydrodynamic performance, to obtain the load characteristics of the turbine.

    (3) Based on the full scale turbine experiment, to validate the accuracy of the numerical simulation results.

    (4) To solve some unpredictable issues with regard to the design and manufacture process, which would be significant for the generator working in the real sea.

    (5) Compared with the model turbine, to provide information about the scaling effect.

    In this paper, all objectives mentioned above are discussed except for point two, which would be discussed in the next paper.

    Table 2 Parameters of full scale turbine

    The rig layout is shown in Fig.4. The generator (1#) is connected to an elliptic column (5#) by a flange (2#), the elliptic column is connected to the circular column (12#), and then they are connected to the trough girders (13#) and are welted in the girder of the carrier (22#) by the trough girders (21#). The foam layout is shown in Fig.4. The length of the foam is 2.0 m, and its diameter is 0.4 m. A picture of the full scale turbine experiments carried out in basin is given in Fig.5.

    Fig.4 The rig layout of full scale turbine

    Fig.5 (Color online) Photo of the full scale turbine experiment in basin

    1.4Full scale turbine test in real sea

    The full scale turbine is tested under the real sea condition, in the Xiaomentou channel, Daishan island, China. The maximum flow speed is about 2.5 m/s, the minimum flow speed is about 1.2 m/s, and the average flow speed is about 1.7 m/s. The highest water level is about 7.5 m, while the lowest water level is about 3.6 m/s. The seabed is covered with stone. All environment parameters mentioned above were measured by the Second Institute of Oceanography, SOA.

    The tidal current power station is fixed on the seabed, named “Haiming I” and developed by HEU. The generator is supported by a triangular basement and a hexagon framework, and the generator can rotate around the upper and lower points of the hexagon framework, as shown in Fig.6. The advantage of this device is that the generator can rotate following the incoming flow in the range of ±180o. The full scale turbine working in the real sea is shown in Fig.7. The objectives of the test are:

    (1) To validate the structural strength of the tidal current power station under the real sea condition.

    (2) To validate the characteristic performance of the turbine, by the measurement of the flow speed and the power output.

    (3) To validate the control system, to see whether it realizes the optimal power controlling or not.

    (4) To validate the reliability of the seal components.

    In this paper, all objectives mentioned above are discussed except for point three, which will be discussed in the next paper.

    Fig.6 Full scale turbine layout

    Fig.7 (Color online) Photo of full scale turbine working in the real sea

    The test of the full scale turbine in the real sea includes two parts: one is to measure the flow speed using Aquadopp profiler 2 MHz with a standard head, which is 200 m depth rating, 9 MB of memory, with pressure sensors (0 m-500 m), compasses and tiltautomatic detectors up and down, 10 m polyurethane cables of 8 pin connector, the other is to measure the power output using software by recording the voltage and the current. The “Haiming” tidal current power station operates from November 2011, and it is shownthat the structural strength and the seal components of the turbine can satisfy the related requirements.

    2. Results and discussions

    2.1Model test results

    The power coefficient of the turbine is defined as

    The thrust coefficient of the turbine is defined as

    wherePis the power output of the turbine,Dis the diameter of the turbine,Tis the thrust force of the turbine,Uis the free stream velocity, which is equivalent to the speed of carriage,ρis the fluid density.

    The tip speed ratio of the turbine is given as

    whereωis the angular velocity of the turbine.

    The hydrodynamic performances of the turbine are described by usingCp-λcurves andP-ωcurves, which indicate the turbine’s ability of capturing the tidal current energy. Due to the limitation of the experiment model, there is not enough space to install the torque meter in the cabin, therefore, only power characteristic curves are obtained. The output power of the model turbine is measured at different flow speeds.

    The incoming flow speed influences the Renoyld’s number of the blade, thus the hydrodynamic performances of the turbine would be affected by the incoming flow speed. During the experiment, the power output of the turbine is measured for different tip speed ratios and pitch angles. The tip speed ratio is dimensionalised by the incoming flow speed and the rotational speed of the turbine, controlled by changing the speed of the towing carrier and the rotational speed of the turbine in practice. TheCp-λcurves at different incoming flow speeds (1.0 m/s, 1.5 m/s and 2.0 m/s) are shown in Fig.8. Fig.8 shows that when the speed ratio is about 5, the maximum value ofCpis reached. The rotational speed of the turbine increases with the increase of the incoming flow speed for a constant tip speed ratio, and the trends of curves at different incoming flow speeds are similar. However, with the increase of the incoming flow speed, the maximum value ofCpdecreases gradually. There might be two reasons to explain the results mentioned above. One is that with the increase of the incoming flow speed, the free surface attack is severed, which affects the working condition of the turbine, the other one is that with the increase of the rotating speed of the turbine, cavitations might occur on the surface of the blade. The pitch angles tested in experiment are 0o,?2.5o,2.5oand 5o. The experiment results are shown in Fig.9. For the design turbine, when the pitch angle is0o, the power performance of the turbine is optimal. Therefore, the full scale turbine is a fixed pitch turbine, with the pitch angle of0o.

    Fig.8 Power coefficients with tip speed ratio at different flow speeds

    Fig.9 Power coefficients with tip speed ratio at different pitch angles

    2.2Full scale turbine test results in basin

    The full scale turbine experiment is carried out in a basin and the incoming flow speeds are 1.5 m/s, 1.8 m/s and 2.0 m/s. The axial force of the turbine increases with the increase of the tip speed ratio, and at the high range of the tip speed ratio, the axial force becomes too large for the towing carrier to support, therefore, only a part of experiment data in the low range of the speed ratio is obtained. The test results are shown in Fig.10. When the incoming flow speed is 1.5 m/s and 1.8 m/s, the output power of the turbine is2.0 kW and 3.6 kW, and the maximum power coefficients are 38.0% and 38.8%, respectively. When the current speed is 2.0 m/s, the loads become too great, resulting in a strong vibration of the support structure, so only data for one point are obtained (whenn= 70 rpm,P=4.9 kW), and the power coefficient is about 39% at this point.

    Fig.10 Power performance of full scale turbine

    Based on the full scale turbine test in the basin, the power and load characteristics are checked, and some unpredictable problems are found during testing. When the generator is working in a high rotation speed, the output voltage is too large to make the control system operate normally, which is another reason why the experimental results are not available in the high rotation speed range, and this issue can be solved by increasing the rated voltage of the control system. During the full scale turbine testing, the starting torque of the generator is tested, and the results show that the static magnetic torque is very large, to make the turbine a little difficult to start, and this issue should be solved during the design process of the generator. The full scale turbine experiment is the basis for the turbine working in the real sea. The comparison between the full scale turbine and the model turbine is important for the study of the scale effect, to provide useful information for the design and operation of larger scale tidal current power stations.

    Fig.11 Comparison between model turbine (1 kW) and full scale turbine (10 kW)

    For the tidal current turbine, the offset of the curves due to the scale effect should not be ignored, and it is important to determine the accuracy of the scale effect. The comparison between the model turbine results and the full scale turbine results is shown in Fig.11. The maximum power coefficient (Cpmax) of the full scale turbine (D=2.0 m)is higher than that of the model turbine(D=0.7 m), by about 6%. And the speed ratio with respect to the maximum power coefficient(Cpmax)in the full scale is a little smaller than that of the model turbine. With the increase of the diameter of the turbines, the power coefficient of the turbine increases, and the speed ratio responding to the maximum power coefficient decreases, and the reason might be the decrease of the friction loss. The comparison between the model turbine and the full scale turbine could provide the results of the scale effect, which would be the basis to study the scale effect.

    Fig.12 Current speed varied with time

    Fig.13 Probability density of current speed

    2.3Full scale turbine test results in sea trial

    The real sea test of the full scale turbine is carried out in Daishan, China. The distribution of the flow speed in the sea depends on the season and the lunar calendar date. The flow speed is the largest in September of the year; on the first day and the fifteen day of the lunar calendar, when one sees the largest current speed. Figure 12 shows the characteristics ofthe tidal current in 1 d (November 7, 2011). It is obvious that the flood (negative value) or the ebb (positive value) lasts for 6 h and 12 min, and the curve is similar to a sine curve. Figure 13 is the probability distribution of the current speed. In a day, the probability of the current speed in the range of 1.6 m/s-1.8 m/s is 35%, which takes the largest part.

    Figure 14 is the test results for the full scale turbine in the real sea. The output power of the generator always varies with the incoming flow speed at any time. Figure 15 is the comparison between the real test results and the design values. It is obvious that the test results are in good agreement with the design values. It is estimated that when the incoming flow speed is 2.3 m/s, the output power is nearly 10 kW, and the power coefficient of the generator is about 36.3%.

    Fig.14 (Color online) Power performance with current speed (October 31)

    Fig.15 Comparison between test results and calculation results

    3. Conclusions

    Experiments of the model turbine and the full scale turbine are carried out. The rigs are designed, manufactured and installed successfully for measuring the performance characteristics of the model tidal current turbine in the towing tank and the full scale current turbine in the basin and the real sea, respectively.

    The power coefficient is related with the incoming flow speed and the pitch angle of the blade. The maximum power coefficient decreases with the increase of the incoming flow speed, for the designed model turbine, when the pitch angle of the blade is 0o, the power performance is the best.

    A comparison between the model turbine and the full scale turbine shows that the power coefficients increase with the increase of the sizes of the turbines. For the designed turbine, the maximum power coefficient of the full scale turbine is larger than that of the model turbine, by about 6%, and the optimal characteristics of the turbine are achieved at a lower speed ratio with the increase of the scale. The results shed some light on the scaling effect for the further research. The comparison between the real sea test results and the design values of the full scale turbine shows that the test results are in good agreement with the design values, so the designed full scale turbine provides a good reference for the further study.

    [1] Bahaj A. S., Molland A. F., Chaplin J. R. Power and thrust measurements of marine current turbines under various hydrodynamic flow conditions in a cavitation tunnel and a towing tank [J].Renew Energy, 2007, 32(3): 407-426.

    [2] Bahaj A. S., Batten W. M. J., McCann G. Experimental verifications of numerical predictions for the hydrodynamic performance of a horizontal axis marine current turbines [J].Renew Energy, 2007, 32(15): 2479-2490.

    [3] Batten W. M. J., Bahaj A. S., Molland A. F. Experimentally validated numerical method for the hydrodynamic design of horizontal axis tidal turbines [J].Ocean Engineering, 2007, 34(7): 1013-1020.

    [4] Charles W. F., Roger C. Electrical power generation from ocean currents in the straits of florida: Some environment considerations [J].Renewable and Sustainable Energy Reviews, 2009, 13(9): 2597-2604.

    [5] Lust E. E., Luznik L., Flack K. A. et al. The influence of surface gravity waves on marine current turbine performance [J].International Journal of Marine Energy, 2013, 3-4: 27-40.

    [6] Jing F., Sheng Q., Zhang L. Experimental research on tidal current vertical axis turbine with variable-pitch blades [J].Ocean Engineering, 2014, 88: 228-241.

    [7] Zhang X. W., Wang S. Q., Wang F. et al. The hydrodynamic characteristics of free variable-pitch axis tidal turbine [J].Journal of Hydrodynamics, 2012, 24(6): 834-839.

    [8] Batten W. M. J., Bahaj A. S., Molland A. F. et al. The prediction of the hydrodynamic performance of marine current turbines [J].Renewable Energy, 2008, 33(5): 1085-1096.

    [9] Batten W. M. J., Bahaj A. S., Molland A. F. et al. Hydrodynamics of marine current turbines [J].Renewable Energy, 2006, 31(2): 249-256.

    [10] Wang X. H., Xiao G., Zhang L. Hydrodynamic research for horizontal axis turbine at different rotor preset angles [J].Acta Energiae Solaris Sinca, 2015, 36(5): 1253-1258(in Chinese).

    [11] Maganga F., Germain G., King J. et al. Experimental characteristion of flow effects on marine current turbine behaviour and on its wake properties [J].IET Renewable Power Generation, 2010, 4(6): 498-509.

    [12] Galloway P. W., Myers L. E., Bahaj A. S. Studies of scale turbine in close proximity to waves [C].3rd InternationalConference on Ocean Energy. Bilbao, Spain, 2010.

    [13] Barltrop N., VAryani K., Grant A. et al. Investigation into wave-current interaction in marine current turbine [J].Proceedings of the Institution of Mechanical Engineering, Part A: Journal of Power Energy, 2007, 221(2): 233-242.

    [14] Zhang L., Li Y. B. Fluid mechanics [M]. Harbin Engineering University Press, 2001, 179-188(in Chinese).

    [15] Liu Y. Z. Theory of ship wave making resistance [M]. National Defence Industy Press, 2003(in Chinese).

    [16] Sheng Q. H., Zhao D. Y., Zhang L. A design and numerical simulation of horizontal tidal turbine [J]. Journal of Harbin Engineering University, 2014, 35(4): 389-394(in Chinese).

    [17] Bao D., Liu Z. Z., Tian R. Designing on the blades of 1 kW bousehold wind turbine and experimental study on its performance [J].Renewable Energy Resources, 2008, 26(5): 90-95.

    [18] Lu H. Experimental study on aerodynamic performance of horizontal axis wind turbine and design [D]. Master Thesis, Nanjing, China: Nanjing University of Aeronutics and Astronautics, 2009(in Chinese).

    [19] Li Z. C. Numerical simulation and experimental study on hydrodynamic characteristic of vertical axis tidal turbine [D]. Doctoral Thesis, Harbin, China: Harbin Engineering University, 2012(in Chinese).

    [20] Li L. X., Liao H. S., Liu D. et al. Experimental investigation of the optimization of stilling basin with shallow water cushion used for Froude number energy dissipation [J].Journal of Hydrodynamics, 2015, 27(4): 522-529.

    * Project supported by the National Natural Science Foundation of China (Grant Nos. 51409057, 51579055 ), the Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20132304110009) and the Natural Science Foundation of Heilongjiang Province (Grant No. E2015048).

    Biography:Feng-mei Jing (1982-), Female, Ph. D., Associate Professor

    Corresponging author:Wei-jia Ma, E-mail: maweijia@hrbeu.edu.cn

    猜你喜歡
    張亮
    “失蹤”的母親
    公路橋梁靜載試驗檢測技術(shù)探討
    Spin-orbit stable dirac nodal line in monolayer B6O
    No.5 張亮和楊國福誰的“人設(shè)”更討喜?
    彩虹
    離家的孩子
    東方劍(2017年6期)2017-11-13 20:51:23
    順風(fēng)車
    賣房
    捕獲、跟蹤、瞄準(zhǔn)系統(tǒng)中光斑探測相機(jī)的定位精度
    新晉男神之張亮
    東方電影(2014年2期)2014-04-29 00:44:03
    久久精品国产99精品国产亚洲性色 | 人人妻人人澡人人看| 亚洲成国产人片在线观看| ponron亚洲| 大型黄色视频在线免费观看| 亚洲第一青青草原| 久久久久久久精品吃奶| 国产麻豆69| 久久人妻av系列| 亚洲一区中文字幕在线| 国产精品综合久久久久久久免费 | 亚洲一区高清亚洲精品| 中文字幕高清在线视频| 一区二区三区激情视频| 亚洲激情在线av| 亚洲av熟女| bbb黄色大片| 久久国产精品人妻蜜桃| 少妇 在线观看| 日韩 欧美 亚洲 中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美激情综合另类| 看黄色毛片网站| 老司机靠b影院| 精品国产美女av久久久久小说| 久久精品国产亚洲av高清一级| 日本vs欧美在线观看视频| 久久青草综合色| 亚洲熟女毛片儿| 亚洲精品久久成人aⅴ小说| 满18在线观看网站| 欧美黄色淫秽网站| 国产伦一二天堂av在线观看| 黑人巨大精品欧美一区二区蜜桃| 天堂动漫精品| 成年人黄色毛片网站| 国产人伦9x9x在线观看| 亚洲欧美日韩另类电影网站| 他把我摸到了高潮在线观看| 亚洲国产精品成人综合色| 最新在线观看一区二区三区| 搡老熟女国产l中国老女人| www.自偷自拍.com| 视频区欧美日本亚洲| 午夜精品在线福利| 国产成人av激情在线播放| 亚洲一区中文字幕在线| 黄色女人牲交| 麻豆av在线久日| 国产精品久久电影中文字幕| 精品福利观看| 一进一出好大好爽视频| 桃红色精品国产亚洲av| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 亚洲av美国av| 99在线视频只有这里精品首页| 成人18禁高潮啪啪吃奶动态图| 少妇 在线观看| 他把我摸到了高潮在线观看| 午夜免费激情av| 国产精品久久电影中文字幕| 国产视频一区二区在线看| av视频免费观看在线观看| 国产高清有码在线观看视频 | 色综合站精品国产| 日本一区二区免费在线视频| 一级毛片高清免费大全| 亚洲伊人色综图| 久久久久久国产a免费观看| 别揉我奶头~嗯~啊~动态视频| 91老司机精品| 亚洲国产欧美日韩在线播放| 亚洲久久久国产精品| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区字幕在线| 人人妻人人澡欧美一区二区 | 麻豆一二三区av精品| 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 夜夜爽天天搞| 亚洲欧洲精品一区二区精品久久久| 久热爱精品视频在线9| 婷婷六月久久综合丁香| 成年人黄色毛片网站| 免费在线观看日本一区| 97碰自拍视频| 久久中文字幕一级| 欧美中文综合在线视频| 侵犯人妻中文字幕一二三四区| 色综合亚洲欧美另类图片| 亚洲国产精品成人综合色| 法律面前人人平等表现在哪些方面| 亚洲午夜精品一区,二区,三区| 国产成人系列免费观看| 亚洲一码二码三码区别大吗| 制服人妻中文乱码| 国产一区二区三区视频了| 他把我摸到了高潮在线观看| 国产欧美日韩一区二区精品| 波多野结衣巨乳人妻| 亚洲精品一区av在线观看| 俄罗斯特黄特色一大片| 人妻久久中文字幕网| 一二三四在线观看免费中文在| 久久天堂一区二区三区四区| 男女下面插进去视频免费观看| 国产精品 欧美亚洲| 九色亚洲精品在线播放| 久久中文看片网| 欧美色视频一区免费| 国产成人免费无遮挡视频| 日韩欧美国产一区二区入口| 一个人免费在线观看的高清视频| 日韩一卡2卡3卡4卡2021年| 精品久久久精品久久久| 国产成人系列免费观看| 日韩大码丰满熟妇| 日日夜夜操网爽| 自线自在国产av| 天天添夜夜摸| 欧美日韩亚洲国产一区二区在线观看| 成人特级黄色片久久久久久久| 国产精品一区二区三区四区久久 | 欧美丝袜亚洲另类 | 黄频高清免费视频| 成人18禁在线播放| 999久久久国产精品视频| 亚洲无线在线观看| 久久草成人影院| 久久精品亚洲精品国产色婷小说| 国产精品久久久人人做人人爽| www.www免费av| 国产精品电影一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av高清一级| 在线播放国产精品三级| 黄网站色视频无遮挡免费观看| 天堂动漫精品| 很黄的视频免费| 国产亚洲精品一区二区www| 精品久久久久久,| 69av精品久久久久久| 国产aⅴ精品一区二区三区波| 国产一级毛片七仙女欲春2 | 国内精品久久久久久久电影| 天天一区二区日本电影三级 | 亚洲中文日韩欧美视频| 丁香欧美五月| 香蕉国产在线看| 国产1区2区3区精品| 男人舔女人的私密视频| 亚洲七黄色美女视频| 国产精品永久免费网站| 少妇粗大呻吟视频| 99国产综合亚洲精品| 欧美中文综合在线视频| 午夜福利高清视频| 少妇裸体淫交视频免费看高清 | 午夜免费激情av| 99国产精品99久久久久| 国产在线观看jvid| tocl精华| 国产免费男女视频| 亚洲久久久国产精品| 亚洲欧美激情综合另类| 亚洲五月色婷婷综合| 99在线人妻在线中文字幕| 日本vs欧美在线观看视频| 免费无遮挡裸体视频| 欧美激情 高清一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 超碰成人久久| 18禁国产床啪视频网站| 黄色片一级片一级黄色片| 欧美乱色亚洲激情| 精品乱码久久久久久99久播| 非洲黑人性xxxx精品又粗又长| 欧美性长视频在线观看| 亚洲久久久国产精品| 女人精品久久久久毛片| 亚洲三区欧美一区| 久久久国产精品麻豆| 大码成人一级视频| 亚洲精品中文字幕在线视频| 高清黄色对白视频在线免费看| 久久影院123| 自线自在国产av| 免费搜索国产男女视频| 午夜日韩欧美国产| 日本 av在线| tocl精华| 一边摸一边抽搐一进一出视频| 国产精品免费一区二区三区在线| 操美女的视频在线观看| 中文亚洲av片在线观看爽| 久久久久久久午夜电影| 日本三级黄在线观看| 黄色 视频免费看| 久久中文字幕一级| 中文字幕久久专区| 亚洲国产看品久久| 欧美日韩福利视频一区二区| 国产精品98久久久久久宅男小说| 久久精品人人爽人人爽视色| 黄色丝袜av网址大全| 精品电影一区二区在线| 欧美在线一区亚洲| 99国产精品99久久久久| 国产精品,欧美在线| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 大陆偷拍与自拍| 欧美大码av| 国产精品影院久久| 欧美成狂野欧美在线观看| 欧美乱色亚洲激情| 91精品三级在线观看| 首页视频小说图片口味搜索| 老汉色∧v一级毛片| 亚洲av第一区精品v没综合| 国产成人av激情在线播放| 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 国产一级毛片七仙女欲春2 | 午夜久久久久精精品| 久久人人97超碰香蕉20202| 欧美成人性av电影在线观看| 91精品国产国语对白视频| 丁香欧美五月| 国产成人啪精品午夜网站| 女人被躁到高潮嗷嗷叫费观| av超薄肉色丝袜交足视频| 亚洲在线自拍视频| 美女免费视频网站| av福利片在线| 久热爱精品视频在线9| 国产av精品麻豆| 高清在线国产一区| 久久国产精品男人的天堂亚洲| 国产av又大| 久久精品亚洲熟妇少妇任你| av有码第一页| 色尼玛亚洲综合影院| 乱人伦中国视频| 丝袜人妻中文字幕| 国产精品乱码一区二三区的特点 | 国产成+人综合+亚洲专区| 18禁国产床啪视频网站| 国产精品亚洲av一区麻豆| 亚洲精品在线美女| 免费观看人在逋| 欧美一级a爱片免费观看看 | 搡老岳熟女国产| 真人一进一出gif抽搐免费| 91九色精品人成在线观看| 18禁国产床啪视频网站| www.精华液| 女人高潮潮喷娇喘18禁视频| 嫁个100分男人电影在线观看| 久久精品亚洲精品国产色婷小说| 午夜福利免费观看在线| 欧美一级a爱片免费观看看 | avwww免费| 99久久国产精品久久久| 日韩中文字幕欧美一区二区| 亚洲国产精品成人综合色| 黑人操中国人逼视频| 宅男免费午夜| 丰满人妻熟妇乱又伦精品不卡| 亚洲自拍偷在线| 波多野结衣一区麻豆| 曰老女人黄片| 香蕉丝袜av| 免费在线观看日本一区| 88av欧美| 在线播放国产精品三级| 亚洲第一av免费看| 国产国语露脸激情在线看| 在线观看免费视频日本深夜| 欧美在线一区亚洲| 女人爽到高潮嗷嗷叫在线视频| 99精品欧美一区二区三区四区| 精品日产1卡2卡| 久久人人97超碰香蕉20202| 国产亚洲av嫩草精品影院| 欧美日韩瑟瑟在线播放| 国产精品自产拍在线观看55亚洲| 国产激情欧美一区二区| 国产av一区二区精品久久| 女人精品久久久久毛片| 国产高清视频在线播放一区| 一级,二级,三级黄色视频| 欧美国产精品va在线观看不卡| 亚洲中文av在线| 黄色丝袜av网址大全| 亚洲国产精品999在线| 青草久久国产| 又紧又爽又黄一区二区| 免费搜索国产男女视频| 亚洲午夜精品一区,二区,三区| 叶爱在线成人免费视频播放| 久久久久久久久中文| 长腿黑丝高跟| 久久热在线av| 他把我摸到了高潮在线观看| 美女扒开内裤让男人捅视频| 国产在线观看jvid| 在线av久久热| 久久久久亚洲av毛片大全| 男人的好看免费观看在线视频 | 熟妇人妻久久中文字幕3abv| 亚洲无线在线观看| 看黄色毛片网站| 91精品三级在线观看| 国产人伦9x9x在线观看| 丰满的人妻完整版| 亚洲第一欧美日韩一区二区三区| 精品人妻1区二区| 久9热在线精品视频| 黄色成人免费大全| 中亚洲国语对白在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产午夜福利久久久久久| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| 亚洲七黄色美女视频| 国产熟女xx| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 免费少妇av软件| a级毛片在线看网站| 国产精品野战在线观看| 99国产精品免费福利视频| 无人区码免费观看不卡| 亚洲五月色婷婷综合| a级毛片在线看网站| 久久久久国产精品人妻aⅴ院| 嫩草影院精品99| 亚洲中文日韩欧美视频| 精品久久久久久久毛片微露脸| 日日爽夜夜爽网站| 欧美日韩中文字幕国产精品一区二区三区 | 日本黄色视频三级网站网址| 久久精品成人免费网站| 黄色成人免费大全| www.www免费av| 国产精品98久久久久久宅男小说| 国产精品久久视频播放| 国产精品亚洲一级av第二区| 亚洲五月婷婷丁香| 久久婷婷成人综合色麻豆| 国产精华一区二区三区| 亚洲精品中文字幕在线视频| av福利片在线| 嫩草影视91久久| 99在线视频只有这里精品首页| 不卡一级毛片| 精品欧美国产一区二区三| 国产视频一区二区在线看| 亚洲精品久久国产高清桃花| 一级黄色大片毛片| 桃红色精品国产亚洲av| 亚洲国产欧美日韩在线播放| 国产极品粉嫩免费观看在线| 亚洲全国av大片| 别揉我奶头~嗯~啊~动态视频| 老汉色∧v一级毛片| 亚洲成av人片免费观看| 亚洲成人久久性| 如日韩欧美国产精品一区二区三区| 一边摸一边做爽爽视频免费| 又大又爽又粗| 亚洲三区欧美一区| 精品欧美国产一区二区三| 国产99白浆流出| tocl精华| 好男人电影高清在线观看| 精品熟女少妇八av免费久了| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 日本五十路高清| 精品一区二区三区四区五区乱码| 国产亚洲精品第一综合不卡| 亚洲欧美精品综合久久99| 国产精品久久久av美女十八| 法律面前人人平等表现在哪些方面| 国产高清激情床上av| 亚洲国产精品久久男人天堂| 精品乱码久久久久久99久播| 脱女人内裤的视频| 女性生殖器流出的白浆| 在线国产一区二区在线| 欧美日韩亚洲综合一区二区三区_| 搡老熟女国产l中国老女人| 亚洲九九香蕉| 曰老女人黄片| 亚洲人成电影观看| 色哟哟哟哟哟哟| 久久久国产成人免费| 精品电影一区二区在线| 操美女的视频在线观看| 午夜福利欧美成人| 国产高清videossex| 精品国产美女av久久久久小说| 中文字幕最新亚洲高清| 在线免费观看的www视频| 久久人妻av系列| 久久亚洲真实| 在线观看66精品国产| 一本久久中文字幕| 亚洲精品在线美女| 国产麻豆成人av免费视频| 久久久久亚洲av毛片大全| 国产亚洲欧美98| 无遮挡黄片免费观看| 亚洲欧美精品综合久久99| 国产精品香港三级国产av潘金莲| 成人18禁高潮啪啪吃奶动态图| 日韩欧美三级三区| 成人欧美大片| 国产精品日韩av在线免费观看 | 两个人视频免费观看高清| 熟妇人妻久久中文字幕3abv| 久久久久久久精品吃奶| 亚洲 欧美 日韩 在线 免费| 欧美人与性动交α欧美精品济南到| 黄色 视频免费看| 午夜老司机福利片| 午夜福利18| 成年女人毛片免费观看观看9| 午夜福利欧美成人| 国产精品永久免费网站| 老汉色av国产亚洲站长工具| 窝窝影院91人妻| 亚洲中文av在线| 亚洲一区二区三区不卡视频| 国产伦人伦偷精品视频| 岛国在线观看网站| 国产成人精品久久二区二区91| 在线十欧美十亚洲十日本专区| 免费人成视频x8x8入口观看| 在线观看免费视频日本深夜| 国产欧美日韩综合在线一区二区| 村上凉子中文字幕在线| 在线播放国产精品三级| 国产黄a三级三级三级人| 亚洲中文日韩欧美视频| 国产成人精品久久二区二区91| 69精品国产乱码久久久| 亚洲熟女毛片儿| a在线观看视频网站| 国产亚洲av高清不卡| 波多野结衣一区麻豆| 此物有八面人人有两片| 首页视频小说图片口味搜索| avwww免费| 久久久久久久久免费视频了| 变态另类丝袜制服| 在线观看舔阴道视频| 亚洲一码二码三码区别大吗| 亚洲 欧美 日韩 在线 免费| 久久久精品国产亚洲av高清涩受| 97人妻精品一区二区三区麻豆 | 高潮久久久久久久久久久不卡| 黄色视频不卡| 欧美国产精品va在线观看不卡| 一级片免费观看大全| 午夜两性在线视频| 久久这里只有精品19| 国产又爽黄色视频| 久久国产精品人妻蜜桃| 免费高清在线观看日韩| 久久久久国产精品人妻aⅴ院| 午夜福利在线观看吧| 黄色 视频免费看| 久久久精品国产亚洲av高清涩受| 亚洲第一av免费看| 欧美久久黑人一区二区| 欧美激情极品国产一区二区三区| 欧美乱色亚洲激情| 日本 欧美在线| 激情视频va一区二区三区| 国产极品粉嫩免费观看在线| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 国产亚洲欧美在线一区二区| 熟妇人妻久久中文字幕3abv| 大香蕉久久成人网| 悠悠久久av| 国产成人一区二区三区免费视频网站| 亚洲第一av免费看| 高清在线国产一区| 国产精品日韩av在线免费观看 | 精品欧美国产一区二区三| 乱人伦中国视频| 亚洲国产中文字幕在线视频| 亚洲,欧美精品.| 亚洲av片天天在线观看| 亚洲精华国产精华精| 欧美亚洲日本最大视频资源| 中文字幕色久视频| 丰满的人妻完整版| 精品高清国产在线一区| 日本 欧美在线| 一进一出抽搐gif免费好疼| 久久久久久亚洲精品国产蜜桃av| 国产三级在线视频| 天堂影院成人在线观看| 国产午夜精品久久久久久| 久久精品国产综合久久久| 亚洲黑人精品在线| 999精品在线视频| av片东京热男人的天堂| 18禁国产床啪视频网站| 一级毛片高清免费大全| 夜夜爽天天搞| 久久人人97超碰香蕉20202| 亚洲第一欧美日韩一区二区三区| 亚洲久久久国产精品| 9191精品国产免费久久| 久久久久久久久久久久大奶| 91字幕亚洲| 久久久国产成人精品二区| 午夜福利一区二区在线看| 脱女人内裤的视频| av视频在线观看入口| 亚洲成人久久性| 欧美中文日本在线观看视频| 亚洲第一青青草原| 99久久99久久久精品蜜桃| 中国美女看黄片| 久久欧美精品欧美久久欧美| 国产av一区在线观看免费| 女性生殖器流出的白浆| 成人国产综合亚洲| 免费在线观看黄色视频的| 一级a爱片免费观看的视频| 最新美女视频免费是黄的| 欧美日韩亚洲综合一区二区三区_| 国产精品98久久久久久宅男小说| 精品国产一区二区久久| 精品人妻在线不人妻| 丝袜美足系列| 亚洲国产欧美一区二区综合| 日韩精品青青久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲久久久国产精品| 亚洲最大成人中文| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av在线| 国产精品一区二区在线不卡| 国产精品av久久久久免费| 精品国产亚洲在线| 十分钟在线观看高清视频www| 国产91精品成人一区二区三区| 国产99久久九九免费精品| 丝袜美腿诱惑在线| 亚洲 欧美 日韩 在线 免费| 欧美中文日本在线观看视频| 99国产精品免费福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲第一青青草原| 午夜福利影视在线免费观看| 大香蕉久久成人网| 日本欧美视频一区| 一级毛片女人18水好多| 久久精品国产亚洲av高清一级| 午夜视频精品福利| 亚洲avbb在线观看| 国产精品爽爽va在线观看网站 | 成人国产一区最新在线观看| 精品日产1卡2卡| 女性生殖器流出的白浆| 侵犯人妻中文字幕一二三四区| 日本撒尿小便嘘嘘汇集6| 99精品欧美一区二区三区四区| 熟妇人妻久久中文字幕3abv| 在线观看日韩欧美| 国产精品香港三级国产av潘金莲| 午夜日韩欧美国产| 大香蕉久久成人网| 精品国产一区二区三区四区第35| 啦啦啦韩国在线观看视频| 亚洲av熟女| 嫩草影视91久久| xxx96com| 国产成人精品无人区| 日韩欧美国产一区二区入口| 黑人欧美特级aaaaaa片| 最近最新中文字幕大全免费视频| 欧美+亚洲+日韩+国产| 精品第一国产精品| 国产高清有码在线观看视频 | 午夜福利影视在线免费观看| 男人操女人黄网站| 国产高清激情床上av| 女人爽到高潮嗷嗷叫在线视频| 波多野结衣一区麻豆| 国产主播在线观看一区二区| 精品久久久精品久久久| 好男人在线观看高清免费视频 | 久久天躁狠狠躁夜夜2o2o| 最好的美女福利视频网| 免费观看精品视频网站| 亚洲五月天丁香| 久久国产乱子伦精品免费另类| 亚洲精品av麻豆狂野| 又大又爽又粗| 精品国产一区二区三区四区第35| 亚洲中文日韩欧美视频|