Abstract

The paper describes an investigation into the comprehensive sensitivities of Sino high strength hydraulic concrete to the addition levels of fly-ash powder and silica fume. Four fly-ash powder addition schemes (100, 120, 140, and 160 kg/m3) and four silica fume addition schemes (12.5, 17.5, 22.5, and 27.5 kg/m3) were considered, respectively. The curing ages’ effects were incorporated in this study. The sensitive feedbacks, including creep development, temperature rise procedure, and strain-stress constitution, have been analyzed. The simulation and test results indicate that the Sino high strength hydraulic concrete, prepared by 140 kg/m3 fly-ash powder and 22.5 kg/m3 silica fume, gives the superb strength, durability, and creep stability.

1. Introduction

Sino hydraulic construction triggered up the heavy application of high strength and fully graded concrete that has been the predominant material for these monstrous structures during the past 30 years [1, 2]. Meanwhile, the mechanical-physical characteristics of the massive high strength hydraulic concrete (HSHC) controlled primarily these hydraulic structures’ working behaviors [3, 4].

However, there is little work to study the comprehensive properties of Sino HSHC material of which the expected amount in the future is still above 6 hundred million m3 in China [5, 6].

The main philosophy for the potential-exploiting of HSHC is to optimize the mixture’s constituents creatively in order to increase its strength and durability and reduce its creep and hydration heat generation (that can cause dreadful damage, fracture, etc.) as much as possible.

Some work paid close attention to the constituents’ effects on the working behavior of the normal concrete. The aggregates could be treated by chemical process to strengthen its connection with the mortar binder [7]; the water curing activity for the aggregates has also been employed in most Sino concrete construction including the work in this paper. Some supplement aggregates, such as furnace slag, were utilized in concrete material that showed distinct cost-advantage [8]; however, the vehement abrasion in the hydraulic environment may be a challenge against it.

Particularly, the study of the special concrete has been one hotspot of material science advances. Some specific additions in the concrete material were adopted and their mechanical-physical effects were also studied. The thermal conductivity could be the key factor for the nuclear reactor concrete [9]; the application of magnetite aggregate was helpful for the thermal conductivity innovation. The pumice material was used to substitute the rock aggregates and the cost could be cut down effectively [10]; according to the relative study, the pumice aggregate concrete could get the maximally compressive strength up to 30 MPa that met the requirements of the conventional concrete construction.

In terms of HSHC study, it is the crucial methodology to qualify the functions of its constituents so as to exploit optimally its potential capacity.

The internal addition method was employed in this paper to study the comprehensive sensitivities of Sino HSHC to the addition levels of silica fume and fly-ash powder. Namely, the cement powder was replaced by these agents directly and variably.

In the present study, the well-grained silica fume and fly-ash powder (the level of their particles’ sizes is 1 × 10−2 mm) were utilized with different addition schemes by mass (12.5, 17.5, 22.5, and 27.5 kg/m3 for silica fume; 100, 120, 140, and 160 kg/m3 for fly-ash powder). Particularly, the curing ages were considered in detail in the thermodynamics and strain-stress constitutive model.

The thermodynamic characteristics and strain-stress constitution were discussed in this study with the diverse addition schemes. Particularly, several updated models on the temperature rise procedure and the thermodynamic creep development as well as the stress path variation were established for HSHC whose properties can be then investigated more appropriately than by the conventional ones.

The second part of this paper stated the constituents’ preparation and the samples’ production as well as the experimental equipment and programs; the third part expressed the methodologies on the comprehensive tests that included the laboratory and the in situ ones; the fourth and fifth parts introduced the primary tests’ results of thermodynamic characteristics and strain-stress constitutive model; the sixth part discussed in detail the systematical sensitivities of HSHC. With the diverse addition schemes of silica fume and fly-ash powder, the multimechanical-physical characteristics including the temperature rise procedure, the thermodynamic creep development, and the strength evolution of HSHC were investigated.

2. Experimental Programs

2.1. Material Constituents of Sino HSHC

We produced the target concrete material for our study in the State Key Laboratory of Hydroscience and Hydraulic Engineering of Tsinghua University. The original place for all the parental materials is the construction site in Southern-Western China where the gravel, sand, cement, and silica fume were collected. We designed the target concrete material’s standard constituents percentages based on their mass proportions (Figure 1). The expected compressive strength of the target concrete material is above 40 MPa.

The information on the parental materials was offered based on their mass proportions.

Gravel and sand, as the aggregate constituents, were taken from the Jinshajiang River catchment (Figure 2). Gravel with the size range of 5–20 mm was excavated from the basalt layer of Jinshajiang mountain. River sand with the relative compaction value of 0.35 was collected from the valley bed [11, 12].

The cement producer in the study is XUANWEI Company in Southern-Western China. The cement powder (Figure 3) with grade 425 of Sino mainland standard is moderate-heat one [13, 14].

The mixing water with pH of 7.0 is Northern China pure water and can be applied instantly for lives.

The fly-ash powder producer is HAICHUANG Company in Nantong city of Jiangsu province. Its average grain diameter is 0.045 mm (Figure 4). The content percentage of SO3 is <3%. As the crucial constituent of Sino high strength hydraulic concrete, it can help diminish effectively the hydration heat generation, in particular, during the early stage of massive concrete pouring procedure [15]. Our work will establish the qualified relation between the sound amount of fly-ash powder and the thermodynamic characteristics of Sino HSHC.

The silica fume producer is LANGTIAN Company of Sichuan province in Southern-Western China. Its average grain diameter is 0.074 mm (Figure 5). Its proper application, in spite of the small amount, is the critical study of massive and high strength hydraulic concrete construction. This paper will study how the silica fume can strengthen the Sino HSHC material [16].

The water reducing agent is the twins agent for silica fume whose application will raise the water requirement and weaken the workability of Sino HSHC material. Hence, the application of water reducing agent can help decrease the water supplement [17]. The water reducing agent in this study is tripolycyanamide series and produced by XIANGBANG Company in Jiangsu province of China.

The heavy seeping pressure is another mortal threat against the Sino HSHC material. Air entraining agent can help reduce the original voids of concrete material and densify the massive concrete block. The air entraining agent in this paper is colophony series and produced by QICHENG Company in Shanghai.

2.2. Tests’ Standards and Equipment

The constituents proportion design and concrete samples’ modeling procedure obey mainly the Sino standards and refer partly to the American codes [13, 1820].

According to the standard constituents percentages, the plastic concrete paste was produced in the mixer of which the type was JG244-2009 (maximal power 3 kW; capacity 60 L) produced by Jingwei Company of Hebei province in China. The standard mixing time was 30 minutes. The homogeneous paste then was poured into the normal moulds (Figure 6; standard dimension, i.e., length × width × height: 100 mm × 100 mm × 300 mm) that were put on the vibroplatform for vibrating compaction (vibrating period 15 minutes; vibration frequency 50 ± 3 Hz; vibration amplitude 0.35 mm).

These concrete samples were maintained in the curing chamber where the temperature and percentage relative humidity (RH) were 20 ± 5°C and 95%, respectively. The curing ages including 7 days and 28 days were different for these tests in this paper.

2.3. Thermodynamic Creep Tests

Figure 7 shows the standard concrete sample and creep loading test. The creep loading test instrument for Sino HSHC material was TXB-50 system produced by GANGYUAN Company in Tianjin. TXB-50 system (Figure 7) is spring instrument (maximal range 20 mm; maximal tensile force 90 MPa).

2.4. Temperature Rise Procedure Tests

The system for temperature rise procedure test is automatic equipment JDC-2 that includes the data collector, sensors, and close circuit. The data of the whole temperature rise procedure can be collected continuously by them (Figure 8).

2.5. Stress Characteristics Tests

Multiaxial compressive loading tests were realized with American Instron-8508 multiaxial loading system by which we simulated physically the whole multiaxial loading procedures of the target material (Figure 9). The maximal compressive force of Instron-8508 is 1 MN. The maximal strokes for vertical and horizontal axis are 900 mm and 750 mm, respectively. The tensiometers (the marked area in Figure 9) are the key supplementary units in these tests to help get the complete stress-strain constitutive curves.

2.6. Systematic Curing Activities

Temperature rise procedure tests were implemented in the construction site. In order to realize the standard curing conditions, the systematic curing activities (Figure 10) were applied, including regular sprinkling; insulation blanket, board, and painting installation; water cooling pipe nets layout; full-time temperature monitoring.

3. Proposed Methodologies

3.1. Thermodynamic Creep Tests

As for the concrete samples with age of 7 days or 28 days, the eventual creep loading duration increment is 200 days that is the generally accepted standard of the concrete creep loading test in Sino mainland. The creep loading tests were applied in the specific chamber where the temperature was kept at 20 ± 5°C and RH was above 75%. The creep loading ages include 3, 7, 10, 15, 30, 45, 60, 90, 120, 150, and 200 days. The standard concrete samples were loaded with one incremental force (i.e., compensating force) at the end of every age and the strain induced by the incremental force was 10−5.

3.2. Temperature Rise Procedure Tests

The temperature rise procedure tests were implemented in one arch dam construction site of Sichuan province in Southern-Western China. The local average temperature during the test was 15°C and RH was above 85%. The tests began as soon as the concrete paste was set and the total age of the test was 90 days. The temperature sensors were eternally embedded in the HSHC paste. The full-time engineers, after the tests, will monitor the temperature evolution during the whole working period of the hydraulic structure.

3.3. Stress Characteristics Tests

Multiaxial compressive loading tests include uniaxial and biaxial compressive loading cases (Figure 9(b) was the uniaxial case; Figure 9(c) was the biaxial case). The multiaxial compressive loading tests adopted strain controlling technique. The strain velocity was 2 × 10−4/minute for the uniaxial compressive loading case. The vertical and horizontal axial strain velocities for the biaxial case were 1 × 10−4/minute and 5 × 10−5/minute, respectively. The temperature and RH during the tests were kept at 15 ± 5°C and 50%. The standard concrete samples were loaded after 28 days with standard curing activities.

3.4. Sensitivity Tests

Based on the standard constituents percentages of Sino HSHC material, how will the addition levels of fly-ash powder and silica fume affect HSHC’s working behavior and comprehensive characteristics including thermodynamics, temperature rise procedure, and strain-stress constitutions? We will study these comprehensive sensitivities on Sino high strength hydraulic concrete material.

4. Thermodynamic Characteristics on High Strength Hydraulic Concrete

The record of gravity dam with maximal concrete consumption is kept by Three Gorges gravity dam that has devoured 28,000 thousand m3 concrete material whose construction period lasted for 18 years. In the next decade, there will come more than 20 monstrous dams to be launched in Sino mainland and their construction period may occupy half of this century. During just this lengthy construction period, the massive concrete block will generate the ghastly amount of hydration heat that can split easily these hydraulic monsters [21]. Particularly, these colossal structures may be destroyed absolutely if the temperature cracks cannot be controlled properly [22, 23]. Therefore, the quantification on Sino HSHC material’s thermodynamic characteristics is the key shield against the hazardous thermodamage development.

The achievement of effective simulation models is the crucial approach for the quantification on Sino HSHC material’s thermodynamic characteristics. The models should help engineers to ascertain the thermodynamic creep development and dynamic temperature rise procedure.

In terms of the thermodynamic creep models, they should be able to compute the recoverable creep and unrecoverable creep, respectively. Most cases of Sino hydraulic construction show that both the recoverable creep and unrecoverable creep accumulation are the lurking peril for temperature cracks development of massive hydraulic structures.

Figure 11 shows the distribution of temperature creep cracks of some typical Sino super dams that include arch dams and gravity ones with height over 150 m. Chinese engineers and researchers adopted the compound exponential model (i.e., the referred model) in the past to realize the thermodynamic creep calculation: where represents the thermodynamic creep values; denotes the static age from the curing start point; refers to the total age that includes the static age and loading period; , , , , , , , and are the fitting parameters to calculate the recoverable creep of high strength hydraulic concrete material; and are the fitting parameters to calculate the unrecoverable creep; is the loading duration increment.

The philosophy of thermodynamic creep calculation is to quantify the recoverable creep and unrecoverable creep, respectively. The hydraulic application experiences told us that both of them could create temperature creep cracks, but they played different roles in different stages.

Based on the standard constituents percentages, the creep loading test and numerical simulation results are shown in Tables 1 and 2.

It can be deduced that the recoverable creep will develop in the early stage of concrete material construction. Although this period is short one, the time accumulation gradient of recoverable creep development is very wild and reaches 0.21 × 10−6/MPa/day. Many mortal temperature cracks may be generated at the stage. On the contrary, the unrecoverable creep development is a long period, but its time accumulation gradient is only 0.009 × 10−6/MPa/day. Surely, the lower time accumulation gradient of unrecoverable creep development does not mean the lower probability. There are still many temperature cracks cases caused by it according to Figure 11.

This paper established a concise model to simulate numerically the hybrid creep development of Sino high strength hydraulic concrete as the following equation: where is the creep coordinate functions; the model can simulate the recoverable or unrecoverable creep development when or , respectively:

in (2) is defined as a loading duration function:

Based on the standard constituents percentages and (2), the comparisons of the creep loading test data and numerical simulation results are showed in Tables 3 and 4. Particularly, the longer the curing time is, the lower the parameters’ values of are.

According to this model, the curing activities will have no effects theoretically over the creep development when the curing age exceeds loading duration increment of recoverable creep. Then, will be a constant one. Therefore, our study is very important because it can optimize basically the constitutive constituents of material by which the creep development of Sino HSHC can be controlled rationally.

Meanwhile, the comparisons between the numerical simulation results of (1) and (2) indicate clearly that the model of this paper has higher accuracy than the referred model (Tables 14 and Figure 12).

The hydration heat discharge of massive concrete material is the intrinsic trigger for temperature creep cracks of hydraulic structures. Hence, the proper calculation of hydration heat discharge is helpful for the precise assessment of hydraulic structures’ thermodynamic characteristics. The temperature rise model was adopted in this paper to qualify the dynamic procedure of massive concrete material’s hydration heat discharge.

Here are four conventional temperature rise models [24], namely, power-exponential composite model (see (5)), hyperbolic model (see (6)), single-exponential model (see (7)), and double-exponential model (see (8)). One has where is the temperature rise variable; is the ultimate temperature rise value; is the early stage temperature rise value; is the mature temperature rise value; , , , and are the test parameters. As for Sino high strength hydraulic concrete material whose compressive strength is above 40 MPa, , , , , , , and .

However, the hydration heat discharge of massive concrete material has close relation with the constitutive constituents of material. Hence, the temperature rise model should incorporate the material’s physical effects. At least, the proper model should consider the concrete material’s sensitivity to its constituents. We established one updated hyperbolic model that will be introduced in the discussion.

5. Stress Characteristics on High Strength Hydraulic Concrete

Most cases of hydraulic concrete structures failure during the operating period were caused by diverse stress overloading conditions [25, 26]. Particularly, the stress states of the massive hydraulic structures are so complicated that the attempt to establish one generalized model will often fail. We will study the targeted stress characteristics of Sino HSHC material and express its stress development behaviors.

Stress-strain constitutive model plays the primary role in translating concrete material’s mechanical behavior [27]. In China, the generally accepted stress-strain constitutive model for concrete material is the bisection nonlinear model which can be expressed as where and are the stress and strain variable, respectively; is the peak strain under compressive loading case; , , and are the concrete material parameters that indicate the continuousness of the bisection nonlinear model at ; and are the shape parameters that indicate the geometrical characteristics of the nonlinear strain-stress curves; is the nominal modulus.

In terms of standard constituents percentages of Sino HSHC material, the simulation values of these constitutive parameters are shown in Table 5.

6. Discussion

The study on comprehensive characteristics sensitivity is valuable for ascertaining the effects of constituents’ application in Sino high strength hydraulic concrete material. The comprehensively sensitive study in this paper was based on internal addition method. Namely, the cement powder was replaced by these agents directly and variably [2830].

Fly-ash powder and silica fume in this paper were the targeted agents to be researched for their effects. Based on the standard constituents percentages, four fly-ash powder addition schemes (100, 120, 140, and 160 kg/m3) and four silica fume addition schemes (12.5, 17.5, 22.5, and 27.5 kg/m3) were considered, respectively. Then, we studied the sensitive feedback of stress and thermodynamic characteristics of Sino high strength hydraulic concrete.

6.1. Temperature Rise Procedure Sensitivity

We monitored the temperature rise procedures of the targeted concrete material with different fly-ash powder addition levels and the data were compared with the four conventional models in Figure 13. According to our study, the fly-ash powder is the indispensable constituent of Sino HSHC. Its application can effectively reduce the hydration heat discharge in the early stage of the concrete pouring procedure.

Figure 13 shows that the values computed by these models converge in the mature stage and the obvious discrepancy just lives in the early stage of hydration heat development. Meantime, the temperature rise procedures of Sino HSHC material show obvious discreteness from the conventional models’ calculation. Hence, the updated hyperbolic model was chosen to express properly the effects of fly-ash powder addition levels:where test parameters and are the functions of fly-ash powder addition level and can be computed bywhere , , , , and .

The theoretical temperature rise procedures based on (10) can be expressed by Figure 14.

In (10), the value of will rise with the increasement of fly-ash powder addition level ; on the contrary, the value of will fall clearly (in Figure 15). Therefore, it is believed that the updated hyperbolic model can incorporate the thermodynamic effect of fly-ash powder addition level over the temperature rise procedures of Sino HSHC.

6.2. Thermodynamic Creep and Stress Characteristics Sensitivities

Figure 16 shows the sensitive curves of creep development of the targeted materials with different ages (7 days and 28 days).

The persistent curing activities are also helpful for the reduction of creep value level. The longer curing age will facilitate the lower level of both recoverable and unrecoverable creep developments. The ribbon of creep development on the silica fume addition levels (Figure 16(b)) is narrower than that on fly-ash powder addition levels (Figure 16(a)). Therefore, the thermodynamic creep of Sino HSHC material is not sensitive to the addition level of silica fume. Under the standard constituents percentages, the cumulative gradient of thermodynamic creep at 7 days age is 0.07 × 10−6 MPa/day that meets the requirements of the international concrete construction code [31]:(a1)Fly-ash powder case under uniaxial compressive loading condition, sample’s age: 28 days.(a2)Silica fume case under uniaxial compressive loading condition, sample’s age: 28 days.(b1)Fly-ash powder case under biaxial compressive loading condition, for vertical axis, sample’s age: 28 days.(b2)Fly-ash powder case under biaxial compressive loading condition, for horizontal axis, sample’s age: 28 days.(b3)Silica fume case under biaxial compressive loading condition, for vertical axis, sample’s age: 28 days.(b4)Silica fume case under biaxial compressive loading condition, for horizontal axis, sample’s age: 28 days.

Figure 17 shows the complete stress-strain constitutive curves under multiaxial compressive loading tests. The mechanical behaviors of the targeted concrete material comply with nonlinear constitutive characteristic under multiaxial compressive loading conditions. Furthermore, the nonlinear constitutive characteristics of Sino HSHC show highly sensitive feedback to the addition levels of the targeted agents, that is, fly-ash powder and silica fume.

Under the standard constituents percentages condition, the values of uniaxial ultimate strength () and vertical axial ultimate strength () keep steadily in the range of 55–65 MPa. However, the horizontal axial ultimate strength () fluctuates clearly with the diverse replacement of the targeted agents: the application of fly-ash powder will degrade the horizontal axial ultimate strength down to 35 MPa; the application of silica fume will promote the horizontal axial ultimate strength up to 70 MPa. The causes for this phenomenon include the active material with higher strength than Ca(OH)2 crystal, namely, 2CaO·SiO2·H2O generated from the silica fume stuffs densely and continually the interfaces of the cements hydrates and the aggregates particles by which the concrete material gets super homogeneity and can bear safely the complex loading case; higher replacement level of fly-ash powder fabricates larger voids between the cements hydrates and the fly-ash particles that will devour more hydrates whose propagation may consume much time and cause particularly higher heterogeneity and lower strength at the early stage.

The degenerative level of ultimate strength rises sharp when the silica fume addition level decreases under uniaxial compressive loading condition: under the case of 27.5 kg/m3, the uniaxial ultimate strength will rocket up to 80 MPa; however, under the case of 12.5 kg/m3, it will reduce down to 45 MPa.

The vertical axial ultimate strength of biaxial compressive loading cases is sensitive to the fly-ash powder addition levels. Its value under the case of 160 kg/m3 is lower than that under the case of 100 kg/m3 by 20 MPa.

The mechanical strain under multiaxial also shows the obvious sensitivities to the addition levels of fly-ash powder and silica fume.

As a whole, the peak strain levels with the silica fume application are higher than that with the fly-ash powder application (Figure 18). The phenomenon means that the application of silica fume helps Sino HSHC develop the continuous deformation with higher level before it is destroyed seriously. The super quality is realized by its homogeneity due to the active material generation. Under the biaxial compressive loading condition, the value of    of the targeted concrete attains above 0.36 (Figure 18(b4)). The value of , even under the uniaxial compressive loading condition that is a more disadvantageous stress state, can reach up to 0.002 (Figure 18(a2)). Therefore, Sino HSHC material with the silica fume application can undergo the rigorous stress states.

The stress path characteristics are the inspiring tools for understanding the quality of Sino HSHC material under complex stress states. The cognition in its stress path development helps scientists and engineers apply properly the HSHC material in complicated circumstances (Figure 19).

The normalized stress space and its coordinate system can be expressed as . The further coordinate transformation can be expressed by (12). The updated stress path can be showed under the transformed coordinate system as in Figure 19:where and are the intermediate coordinate system with a rotation transformation from ,  ; is the rotation angle for normal stress space; and are the eventual coordinate system with a translation transformation from and ; is the translation distance from the intermediate coordinate system and to the eventual coordinate system and ; , , , and are the relative values on stress, namely, nondimensional parameters calculated by (12).

With the same internal addition method, the sensitivity of stress path characteristics of the HSHC material was studied on fly-ash powder and silica fume addition. The results show that the higher the content of silica fume is, the longer the total length of the stress path is; the higher the content of fly-ash is, the shorter the total length of the stress path is. Therefore, the two main agents, that is, fly-ash powder and silica fume, play different roles in the multimechanical-physical characteristics of the high strength hydraulic concrete material. The longer stress paths can form larger areas in stress space and the larger stress space areas will offer higher strength and resistance. Furthermore, the larger stress space areas can generate super durability under the complicate stress states. The results in Figure 19 consist with that in Figure 18.

Although the application of fly-ash powder will weaken the strength of the concrete material, its proper addition level can restrain effectively the thermodynamic creep development (Figures 16 and 20).

The compressive yield strength of the concrete material will decline with the increase of fly-ash powder addition level. By contrast, its creep values will run lower, too. The compressive yield strength of the concrete material will grow with the increase of silica fume addition level. However, the increasing effect tends to converge. Meanwhile, its creep values will reach minimum (17.65 × 10−6/MPa, 7 days age; 15 × 10−6/MPa, 28 days age) with the increase of silica fume addition level (Figure 20).

Due to the double application of fly-ash powder and silica fume, there must be somewhat replacement for the cement amount in Sino HSHC material. However, the replacement influence can be attenuated effectively with the curing age continuing (Figures 20(b) and 20(d)).

Moreover, the compressive yield strengths of Sino HSHC show explicit sensitivities to the curing ages. Under the case of 160 kg/m3 for fly-ash powder addition, the value of compressive yield strength can reach only 17 MPa at the age of 7 days. However, the value can grow up to 29 MPa at the age of 28 days. The cumulative gradient is 0.52 MPa/day. By contrast, under the case of 27.5 kg/m3 for silica fume addition, the cumulative gradient of compressive yield strength is 0.47 MPa/day.

7. Conclusions

Sino hydraulic construction shows both unbeatable ambition and unsurpassable complexity. This paper tried to introduce comprehensively the key aspects of Sino HSHC material. The following conclusions are drawn:(1)Four fly-ash powder addition schemes and four silica fume addition schemes were designed. Based on them, the concrete material’s thermodynamic creep sensitivities to the targeted agents’ addition levels were studied: the completely thermodynamic creep developments at different ages were tested and compared; the creep developments show obvious sensitivities to the fly-ash powder addition levels; in order to qualify the creep development, one new thermodynamic creep model was established and its simulation results show higher accuracy than the old one.(2)The temperature rise procedures of the targeted concrete were monitored on the spot; the measured data were compared with the conventional temperature rise model and the results show explicit discreteness from each other; the updated hyperbolic model was established and its parameters’ functions on the fly-ash powder addition levels were offered.(3)Based on the multiaxial loading tests, the stress-strain constitutive curves of the targeted concrete were completed; the sensitivity of the stress path characteristics was studied; the results prove that the proper addition of silica fume in the targeted concrete can facilitate the strength development. Particularly, the strength sensitivities to the curing ages were studied. The application of fly-ash powder can promote the strength development at the early stage.

Competing Interests

The author declares that there are no competing interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 51109118) and Zhejiang Provincial Natural Science Foundation of China (Grant no. LY14E090001).