Method for Quantitative Analysis of Cavity Zone of the Top of Concrete-Filled Steel Tube

20210310993 · 2021-10-07

    Inventors

    Cpc classification

    International classification

    Abstract

    The invention provides a method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube, comprising the following steps: By substitution of the determined inner radius of the steel tube, the thickness of the tube wall and the propagation speed of the ultrasonic wave in the steel tube, the propagation speed of the ultrasonic wave in the concrete, and the starting time of the first wave when the ultrasonic wave propagating between the top and the bottom of the concrete-filled steel tube into the calculation model of the cavity height of the top of the concrete-filled steel tube, to obtain the cavity height of the top of the concrete-filled steel tube; the calculation model of the cavity height is:

    [00001] t = ( x 1 + d 2 + 2 x 2 ( d + r ) v s + 4 r 2 + d 2 + 4 rd - 2 r h - 2 hd - d 2 + 2 x 2 ( d + r ) v c ) ;

    wherein, x.sub.1 and x.sub.2 are both calculation variables, and their values are:

    [00002] x 1 = d 2 + 2 h d + 2 r h h + d .Math. arcsin ( h + d d 2 + 2 h d + 2 r h ) x 2 = x 3 + x 3 2 - ( - 4 r 2 - d 2 + 2 r h - 4 r d + h d ) ( h 2 - 2 r h ) d 2 - 4 r 2 - d 2 + 2 r h - 4 r d + h d

    Wherein, x.sub.3=−4r.sup.3−rh.sup.2−rd.sup.2+4 hr.sup.2−4dr.sup.2+4drh−dh.sup.2. The method can be used to quantify the cavity height of the top of the concrete-filled steel tube, with small relative errors and high accuracy.

    Claims

    1. A calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube, comprising the following steps: Step 1: Determine the inner radius of Steel Tube (9), the thickness of the tube wall; determine the propagation speed of the ultrasonic wave in Steel Tube (9) and the propagation speed of the ultrasonic wave in Concrete (8); determine the starting time of the first wave when the ultrasonic wave propagating between the top and the bottom of the concrete-filled steel tube; Step 2: By substitution of the determined inner radius of Steel Tube (9), the thickness of the tube wall and the propagation speed of the ultrasonic wave in Steel Tube (9) and the propagation speed of the ultrasonic wave in Concrete (8), and the starting time of the first wave when the ultrasonic wave propagating between the top and the bottom of the concrete-filled steel tube obtained in Step 1 into the calculation model of the cavity height of the top of the concrete-filled steel tube, to obtain the cavity height of the top of the concrete-filled steel tube; the calculation model of the cavity height is t = ( x 1 + d 2 + 2 x 2 ( d + r ) v s + 4 r 2 + d 2 + 4 r d - 2 r h - 2 h d - d 2 + 2 x 2 ( d + r ) v c ) wherein, t is the starting time of the first wave, with the unit as seconds; r is the inner radius of Steel Tube (9), with the unit as meters; d is the wall thickness of Steel Tube (9), with the unit as meters; h is the cavity height of the top of the concrete-filled steel tube, with the unit as meters; v.sub.s and v.sub.c are the propagation speeds of the ultrasonic wave in Steel Tube (9) and in Concrete (8), with the unit as m/s.

    2. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein in Step 1, the method for measuring the starting time of the first wave when ultrasonic propagating between the top and the bottom of the concrete-filled steel tube is as follows: install ultrasonic Sensors (10) on the top and the bottom of the outer wall of Steel Tube (9) using ultrasonic detectors, measure the shortest time of the ultrasonic wave propagating between Sensor (10) at the top and Sensor (10) at the bottom of Steel Tube (9), that is, the starting time of the first wave.

    3. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein by substitution of the cavity height of the top of the concrete-filled steel tube obtained in Step 2 into the cavity arc length calculation “ ” model to obtain the cavity arc length, and the cavity arc length calculation model is α = 4 r .Math. arcsin ( 2 r h 2 r ) , wherein, α is the cavity arc length, with the unit as meters; l is the cavity chord length, with the unit as meters.

    4. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein by substitution of the cavity height of the top of the concrete-filled steel tube obtained in Step 2 into the cavity chord length calculation model to obtain the cavity chord length, and the cavity chord length calculation model is l=2√{square root over (2hr−h.sup.2)}, wherein α is the cavity arc length, with the unit as meters; l is the cavity chord length, with the unit as meters.

    5. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein by substitution of the cavity height of the top of the concrete-filled steel tube obtained in Step 2 into the cavity area calculation model to obtain the cavity area, and the cavity area calculation model is S = r 2 .Math. arcsin ( 2 r h - h 2 r ) - ( r - h ) 2 r h - h 2 Wherein S is the cavity area, with the unit as square meters.

    6. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein Step 2 is calculated via an iterative method.

    7. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein the method for measuring the propagation speed of ultrasonic waves in Steel Pipe (9), Sensors (10) are in good contact with Steel Pipe (9) through the ultrasonic couplant, and Sensors (10) are placed on opposite sides of the outer wall of Steel Pipe (9), N times measurements are performed, and the average value is taken as the ultrasonic propagation speed v.sub.s in Steel Pipe (9), wherein n≥3.

    8. The calculation method for quantitative analysis of cavity zone of the top of the concrete-filled steel tube according to claim 1, wherein the cavity height calculation model is calculated in the radian measure.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0028] FIG. 1 is a distribution diagram of the measuring points for measuring the propagation speed of ultrasonic waves in the steel tube.

    [0029] FIG. 2 is a distribution diagram of the measuring points for measuring the propagation speed of ultrasonic waves in the concrete.

    [0030] FIG. 3 is a schematic diagram of the correlations between the cavity height, the cavity arc length and the cavity chord length.

    [0031] FIG. 4 is a schematic diagram of the propagation path of the ultrasonic wave and the relevant parameters of the calculation model.

    [0032] FIG. 5 is an enlarged schematic diagram of the M block zone in FIG. 4.

    [0033] Reference numbers: 7—cavity zone; 8—concrete; 9—steel tube; 10—sensors.

    DETAILED DESCRIPTION OF THE INVENTION

    [0034] The invention is described in detail in combination with the attached drawings below.

    [0035] In order to make the purpose, technical solutions and advantages of the present invention to be clear, the present invention is further described in detail in combination with the attached drawings and embodiments below. It is understood that the specific embodiments described herein only aims to explain the invention but not to limit it.

    Example 1

    [0036] This example is a specific embodiment of the quantitative analysis method for the Cavity Zone 7 of the top of the concrete-filled steel tube, comprising the following steps:

    [0037] (1) Determine the Basic Parameters of the Cross Section of Steel Tube 9:

    [0038] The inner radius r and the tube wall thickness d of Steel Tube 9 are determined by measurement. The measured results are r=0.069 m, d=0.006 m.

    [0039] (2) Determine the Propagation Speeds of the Ultrasonic Wave in Steel Tube 9 and in Concrete 8:

    [0040] The propagation speed v.sub.s of the ultrasonic wave in Steel Tube 9 and the propagation speed v.sub.c of the ultrasonic wave in Concrete 8 are measured by the ultrasonic detectors respectively. When measuring the propagation speed of the ultrasonic wave in Steel Tube 9, Sensors 10 are in good contact with Steel Tube 9 through the ultrasonic couplant. Sensors 10 are placed on the two opposite sides of the outer wall of Steel Pipe 9 and 4 times of measurements are performed, and the average value is taken as the propagation speed v.sub.s of Steel Pipe 9, the distribution of measuring points {circle around (1)}, {circle around (3)}, {circle around (3)}, {circle around (4)} are shown in FIG. 1. When measuring the propagation speed of the ultrasonic wave in Concrete 8, Sensors 10 are in good contact with the test specimen of Concrete 8 through the ultrasonic couplant. The test specimen of Concrete 8 is a cube specimen of 150 mm×150 mm×150 mm. Choose five measuring points on the surface, and the average value of the measured values is the propagation speed v.sub.c of the ultrasonic wave in Concrete 8. The distribution of the measuring points {circle around (1′)}, {circle around (2′)}, {circle around (3′)}, {circle around (4′)}, {circle around (5′)} is shown in FIG. 2. The measured results are: v.sub.s=5702 m/s; v.sub.c=3843 m/s.

    [0041] (2) Establish the Calculation Model of the Starting Time of the First Wave Based on the Cavity Height of the Top of the Concrete-Filled Steel Tube:

    [0042] According to the inner radius of Steel Tube 9 and the tube wall thickness determined in Step (1) and the propagation speeds of the ultrasonic wave in Steel Tube 9 and in Concrete 8 determined in Step (2), a correlation model of the starting time of the first wave and the cavity height of the top of the concrete-filled steel tube is established as:

    [00007] t = ( x 1 + 0.00 6 2 + 0.15 x 2 5 7 0 2 + 0.020 7 3 6 - 0 . 1 5 h - 0.00 6 2 + 0.15 x 2 3 8 4 3 )

    [0043] Wherein t is the starting time of the first wave, with the unit as seconds (s), h is the cavity height of the top of the concrete-filled steel tube, with the unit as meters (m); x.sub.1 and x.sub.2 are both calculation variables, and their values are:

    [00008] x 1 = 0 . 0 0 0 0 3 6 + 0 . 1 5 h h + 0 . 0 0 6 .Math. arcsin ( h + 0 . 0 0 6 0.000 0 3 6 + 0 . 1 5 h ) x 2 = x 3 + x 3 2 - 0 . 0 0 0 0 3 6 ( - 0 . 0 2 0 7 3 6 + 0 . 1 4 4 h ) ( h 2 - 0 . 1 3 8 h ) - 0 . 0 2 0 7 3 6 + 0 . 1 4 4 h

    [0044] Wherein x.sub.3=−0.001430784−0.075h.sup.2+0.0207h.

    [0045] (4) Determine the Starting Time of the First Wave when the Ultrasonic Wave Propagating Between the Top and the Bottom of the Concrete-Filled Steel Tube:

    [0046] Install ultrasonic wave Sensors 10 on the top and the bottom of the concrete-filled steel tube, and use ultrasonic detectors to measure the shortest time that the ultrasonic wave propagates between Sensor 10 at the top and Sensor 10 at the bottom of the concrete-filled steel tube, that is, the starting time of the first wave t, and the measured results of three specimen groups are respectively: 4.27×10.sup.−5 s, 4.41×10.sup.−5 s, 4.51×10.sup.−5 s;

    [0047] (5) Calculate the Cavity Height of the Top of the Concrete-Filled Steel Tube:

    [0048] According to the calculation model of the starting time of the first wave based on the cavity height of the top of the concrete-filled steel tube established in Step (3), use the starting time of the first wave determined in Step (4), the cavity height of the top of the concrete-filled steel tube is calculated by iterative analysis, the cavity heights are: 0.0104 m, 0.0191 m, 0.0311 m, respectively.

    [0049] (6) Calculate the Cavity Arc Length, the Cavity Chord Length and the Cavity Area of the Concrete-Filled Steel Tube:

    [0050] According to the cavity height of the top of the concrete-filled steel tube calculated in Step (5), and the calculation model of the cavity arc length, the cavity chord length and the cavity area, calculate the cavity arc length, the cavity chord length and the cavity area of the top of the concrete-filled steel tube, taking h=0.0104 m as an example:

    [00009] α = 0.2 76 .Math. arcsin ( 0.001 4 3 5 2 0 . 1 3 8 ) = 0 . 0 768 ( m ) l = 2 0.001 4 3 5 2 - 0.000 1 0 8 1 6 = 0 . 0 729 ( m ) S = 0 . 0 04761 .Math. arcsin ( 0.001327 0 4 0.069 ) - 0 . 0 586 .Math. 0.001 3 2 7 0 4 = 0 . 0 005133 m 2 .

    [0051] Wherein α is the cavity arc length, with the unit as meters (m), l is the cavity chord length, with the unit as meters (m): S is the cavity area, with the unit as square meters (m.sup.2).

    [0052] The calculation results of three groups are listed in Table 1. While, the testing Steel Pipe 9 is split and the actual measurement results obtained by measuring the cross section of the concrete-filled steel tube are also listed in Table 1.

    TABLE-US-00001 TABLE 1 Comparison of calculated values and measured values of the cavity arc length and the cavity chord length By the models of By the models of the present invention Document 1 & Document 2 Measured Calculated Relative Calculated Relative Group value value error value error 1 Cavity height/m 0.0100 0.0104 4.0% 0.0075   25% Cavity arc 0.0752 0.0768 2.1% 0.0653   15% length/m Cavity chord 0.0716 0.0729 1.8% 0.0629   12% length/m Cavity area/m.sup.2 0.0004815 0.0005133 6.6% 0.0003164   34% 2 Cavity height/m 0.0200 0.0191 4.5% 0.0148   26% Cavity arc 0.1080 0.1052 2.6% 0.0921   15% length/m Cavity chord 0.0970 0.0953 1.8% 0.0854   12% length/m Cavity area/m.sup.2 0.001338 0.001252 6.4% 0.0008625   36% 3 Cavity height/m 0.0300 0.0311 3.6% 0.0238   21% Cavity arc 0.1339 0.1365 1.9% 0.1182   12% length/m Cavity chord 0.1138 0.1153 1.3% 0.1043  8.3% length/m Cavity area/m.sup.2 0.002399 0.002525 5.2% 0.001721   28%

    [0053] Through comparison, it is found that the relative errors between the calculated values and the actual measured values of the characteristic values of Cavity Zone 7 of the top of the concrete-filled steel tube analyzed based on the method of the present invention are all less than 7.0%, indicating that the method of the present invention is better than the existing model, and thus verifies the effectiveness and accuracy of the present invention.

    Example 2

    [0054] This example is a specific embodiment of the quantitative analysis method of the top of Cavity Zone 7 of concrete-filled steel tube, comprising the following steps:

    [0055] (1) Determine the Basic Parameters of the Cross Section of Steel Tube 9:

    [0056] The inner radius r and the tube wall thickness d of Steel Tube 9 are determined by measurement. The measured results are r=0.638 m, d=0.022 m.

    [0057] (2) Determine the Propagation Speeds of the Ultrasonic Wave in Steel Tube 9 and in Concrete 8:

    [0058] The propagation speed v.sub.s of the ultrasonic wave in Steel Tube 9 and the propagation speed v.sub.c of the ultrasonic wave in Concrete 8 are measured by the ultrasonic detectors respectively. When measuring the propagation speed of the ultrasonic wave in Steel Tube 9, Sensors 10 are in good contact with Steel Tube 9 through the ultrasonic couplant. Sensors 10 are placed on the two opposite sides of the outer wall of Steel Pipe 9 and 4 times of measurements are performed, and the average value is taken as the propagation speed v.sub.s of Steel Pipe 9, the distribution of measuring points {circle around (1)}, {circle around (2)}, {circle around (3)}, {circle around (4)} is shown in FIG. 1. When measuring the propagation speed of the ultrasonic wave in Concrete 8, Sensors 10 are in good contact with the test specimen of Concrete 8 through the ultrasonic couplant. The test specimen of Concrete 8 is a cube specimen of 150 mm×150 mm×150 mm. Choose five measuring points on the surface, and the average value of the measured values is the propagation speed v.sub.c of the ultrasonic wave in Concrete 8. The distribution of the measuring points {circle around (1′)}, {circle around (2′)}, {circle around (3′)}, {circle around (4′)}, {circle around (5′)} is shown in FIG. 2. The measured results are: v.sub.s=5735 m/s; v.sub.c=3895 m/s.

    [0059] (2) Establish the Calculation Model of the Starting Time of the First Wave Based on the Cavity Height of the Top of the Concrete-Filled Steel Tube:

    [0060] According to the inner radius of Steel Tube 9 and the tube wall thickness determined in Step (1) and the propagation speeds of the ultrasonic wave in Steel Tube 9 and in Concrete 8 determined in Step (2), a correlation model of the starting time of the first wave and the cavity height of the top of the concrete-filled steel tube is established as:

    [00010] t = ( x 1 + 0.02 2 2 + 1.32 x 2 5 7 3 5 + 1.684 8 - 1.32 h - 0.02 2 2 + 1.32 x 2 3 8 9 5 )

    [0061] Wherein t is the starting time of the first wave, with the unit as seconds (s), h is the cavity height of the top of the concrete-filled steel tube, with the unit as meters (m); x.sub.1 and x.sub.2 are both calculation variables, and their values are:

    [00011] x 1 = 0 . 0 0 0 4 8 4 + 1 . 3 2 h h + 0 . 0 2 2 .Math. arcsin ( h + 0 . 0 2 2 0.000 4 8 4 + 1 . 3 2 h ) x 2 = x 3 + x 3 2 - 0 . 0 0 0 4 8 4 ( - 1 . 6 8 4 8 + 1 . 2 9 8 h ) ( h 2 - 1 . 2 7 6 h ) - 1 . 6 8 4 8 + 1 . 2 9 8 h

    [0062] Wherein x.sub.3=−1.074905−0.66h.sup.2+1.68432h.

    [0063] (4) Determine the Starting Time of the First Wave when the Ultrasonic Wave Propagating Between the Top and the Bottom of the Concrete-Filled Steel Tube:

    [0064] Install ultrasonic wave Sensors 10 on the top and the bottom of the concrete-filled steel tube, and use ultrasonic detectors to measure the shortest time that the ultrasonic wave propagates between Sensor 10 at the top and Sensor 10 at the bottom of the concrete-filled steel tube, that is, the starting time of the first wave t, and the measured results of three specimen groups are respectively: 3.608×10.sup.−4 s, 3.649×10.sup.−4 s, 3.711×10.sup.−4 s;

    [0065] (5) Calculate the Cavity Height of the Top of the Concrete-Filled Steel Tube:

    [0066] According to the calculation model of the starting time of the first wave based on the cavity height of the top of the concrete-filled steel tube established in Step (3), use the starting time of the first wave determined in Step (4), the cavity height of the top of the concrete-filled steel tube is calculated by iterative analysis, the cavity heights are: 0.0258 m, 0.0349 m, 0.0523 m, respectively.

    [0067] (6) Calculate the Cavity Arc Length and the Cavity Chord Length of the Concrete-Filled Steel Tube:

    [0068] According to the cavity height of the top of the concrete-filled steel tube calculated in Step (5), and the calculation model of the cavity arc length, the cavity chord length and the cavity area, calculate the cavity arc length, the cavity chord length and the cavity area of the top of the concrete-filled steel tube, taking h=0.0523 m as an example:

    [00012] α = 2.5 52 .Math. arcsin ( 0.066 7 3 4 8 1 . 2 7 6 ) = 0 . 5 203 ( m ) l = 2 0.066 7 3 4 8 - 0.002 7 3 5 2 9 = 0 . 5 060 ( m ) S = 0 . 4 07044 .Math. arcsin ( 0.063 4 0.638 ) - 0 . 5 854 .Math. 0.063 4 = 0 . 0 1779 m 2

    [0069] Wherein α is the cavity arc length, with the unit as meters (m), l is the cavity chord length, with the unit as meters (m): S is the cavity area, with the unit as square meters (m.sup.2).

    [0070] The calculation results of three groups are listed in Table 2. While, the testing Steel Pipe 9 is drilled and the actual measurement results obtained by measuring the cross section of the concrete-filled steel tube are also listed in Table 2.

    TABLE-US-00002 TABLE 2 Comparison of calculated values and measured values of the cavity arc length and the cavity chord length By the models of By the models of the present invention Document 1 & Document 2 Measured Calculated Relative Calculated Relative Group value value error value error 1 Cavity height/m 0.0270 0.0258 4.4% 0.0203   25% Cavity arc 0.3725 0.3641 2.3% 0.3228   13% length/m Cavity chord 0.3673 0.3592 2.2% 0.3194   13% length/m Cavity area/m.sup.2 0.006640 0.006204 6.5% 0.004335   35% 2 Cavity height/m 0.0360 0.0349 3.1% 0.0284   21% Cavity arc 0.4307 0.4240 1.6% 0.3818   11% length/m Cavity chord 0.4226 0.4162 1.5% 0.3761   11% length/m Cavity area/m.sup.2 0. 01020 0.009734 4.6% 0.007160   30% 3 Cavity height/m 0.0510 0.0523 2.6% 0.0439   14% Cavity arc 0.5137 0.5203 1.3% 0.4762  7.3% length/m Cavity chord 0.4999 0.5060 1.2% 0.4652  6.9% length/m Cavity area/m.sup.2 0.01700 0.01779 4.6% 0.01371   19%

    [0071] Through comparison, it is found that the relative errors between the calculated values and the actual measured values of the characteristic values of Cavity Zone 7 of the top of the concrete-filled steel tube analyzed based on the method of the present invention are all less than 7.0%, indicating that the method of the present invention is better than the existing model, and thus verifies the effectiveness and accuracy of the present invention.

    Example 3

    [0072] This example will describe the establishing process of the simplified calculation model in detail with reference to FIGS. 4 and 5.

    [0073] As shown in FIG. 4, the propagation path of the ultrasonic wave in the concrete-filled steel tube in the cavity condition in the present invention is from Point D.fwdarw.A.fwdarw.E.fwdarw.F, wherein the paths through Steel Pipe 9 include D.fwdarw.A and E.fwdarw.F, the path through Concrete 8 is A.fwdarw.E.

    [0074] The path D.fwdarw.A is regarded as a circular arc in the present invention, which is recorded as the calculation variable x.sub.1. In order to find x.sub.1, the radius r′ and the chord length L.sub.DA of this arc need to be known, wherein:

    [0075] The radius of the arc of the path

    [00013] r = d 2 + 2 h r + 2 d h 2 ( h + d )

    [0076] The chord length


    L.sub.DA=√{square root over (d.sup.2+2hr+2dh)}

    [0077] Then calculate the variable

    [00014] x 1 = 2 r ' .Math. arcsin ( L D A 2 r ) = d 2 + 2 h r + 2 dh h + d .Math. arcsin ( h + d d 2 + 2 h r + 2 dh ) .

    [0078] Where r is the inner radius of Steel Tube 9 with the unit as meters (m); d is the thickness of the tube wall with the unit as meters (m); h is the cavity height of the top of the concrete-filled steel tube with the unit as meters (m);

    [0079] The path A.fwdarw.E.fwdarw.F is regarded as a straight line segment in the present invention. To find the lengths of the path A.fwdarw.E and E.fwdarw.F, it is necessary to cross Point E as the perpendicular line of OF and intersect at Point G, and record the length of GH as the calculation variable x.sub.2, as shown in FIG. 5.

    [0080] To find the length of x.sub.2, it needs to establish two equations simultaneously:

    [0081] From the relationship of the circle, we can obtain:


    L.sub.EG.sup.2+(r−x.sub.2).sup.2=r.sup.2;

    [0082] From the relationship of similar triangles, we can obtain:

    [00015] r 2 - ( r - h ) 2 L E G = 2 r - h + d x 2 + d

    [0083] Solve simultaneously to obtain:

    [00016] x 2 = x 3 + x 3 2 - ( - 4 r 2 - d 2 + 2 r h - 4 r d + h d ) ( h 2 - 2 r h ) d 2 - 4 r 2 - d 2 + 2 r h - 4 r d + h d

    [0084] Wherein


    x.sub.3=−4r.sup.3−rh.sup.2−rd.sup.2+4hr.sup.2−4dr.sup.2+4drh−dh.sup.2.

    [0085] So it can be obtained that: the length of AF


    L.sub.AF=√{square root over (4r.sup.2+d.sup.2+4dr−2rh−2hd)};


    the length of EF


    L.sub.EF=√{square root over (d.sup.2+2x.sub.2(d+r))}.

    [0086] The propagation time of the ultrasonic wave in the concrete-filled steel tube should be composed of two parts, namely the propagation time in the steel tube and the propagation time in the concrete. Combined with the path length calculation model, the calculation model of the starting time of the first wave based on the cavity height of the top of the concrete-filled steel tube of the present invention can be obtained, namely the exact model:

    [00017] t = ( x 1 + d 2 + 2 x 2 ( d + r ) v s + 4 r 2 + d 2 + 4 r d - 2 r h - 2 h d - d 2 + 2 x 2 ( d + r ) v c )

    [0087] The descriptions of the patents, patent applications and publications cited in the present invention are all incorporated into the present invention by reference. Any references cited should not be considered as allowing these references to be used as “prior art” in the present invention.

    [0088] The above descriptions are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modification, equivalent replacement and improvement made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.