Abstract:
A method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading is disclosed, which is directed to the cyclic test of loading and unloading. With the assumption that only the deviator stress generates damage to the sample, a test method of the hydrostatic pressure unloading is invented in order to calculate mechanical parameters in different stages of stress and strain. Nine mechanical parameters can be calculated in connection with hexahedral pores connecting samples in the true triaxial test. Six mechanical parameters can be calculated for non-pores connecting samples. Nine mechanical parameters can be calculated in connection with hexahedral pores connecting samples in the traditional triaxial test. Six mechanical parameters can be calculated for non-pores connecting samples. The specific expressions and test methods are provided.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application claims priority from Chinese Patent Application No. 2015100861698, filed on Feb. 17, 2015, in the State Intellectual Property Office of China, the content of which are hereby incorporated by reference in their entirety for all purposes. 
       BACKGROUND OF THE INVENTION 
       [0002]    1. Field of the Invention 
         [0003]    The present invention relates to the technical field of civil engineering and experimental measurement of geological materials, in particular to a method of measurement of stress and strain whole process material parameter by using of hydrostatic pressure unloading. 
         [0004]    2. Description of the Related Art 
         [0005]    The measurement of the rules of variation of geological material whole process parameters is a scientific problem which has not been well solved so far. Its theory and method of measurement is still not sufficiently good. In addition, during the cycle of loading and unloading of the deviatoric stress, the geological material parameters are determined by the approximation to the linear segment of loading according to the test. This often results in large deviation, and even becomes difficult to determine. 
       SUMMARY OF THE INVENTION 
       [0006]    The triaxial test measurement of material parameters has a long history. Furthermore, determining material parameters of rock mass and soil mass in the process of cyclic loading test also has a long history. However, since the current measurement neglects the fact that the deformation change from the hydrostatic pressure occurs, the corresponding material parameters are difficult to be determined when the applied deviatoric stress is greater than the proportional limit stress under the cyclic loading and unloading case. 
         [0007]    The objective of the present invention is to overcome the existing deficiencies in the current techniques. A method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading based on the triaxial test is provided. 
         [0008]    The method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading of the present invention includes the following steps: 
         [0009]    (1.1) applying a hydrostatic pressure first in a true triaxial test, assuming σ 11 =σ 11   H  , σ 22 =σ 22   H , σ 33 =σ 33   H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0 ε jj   H   , i, j ∈ (1,3)   (1)
 
         [0000]    where C iijj   0  is an initial stiffness matrix. 
         [0010]    (1.2) applying a deviator stress q, q=σ 11 +σ 11   H −σ 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (2)
 
         [0000]      σ ii   H   =C   iijj   b ε jj   , i ∈ (2,3),  j ∈(1,3)   (3)
 
         [0011]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space. 
         [0012]    For the material parameters in stress state exceeding yield limit stress space, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 1111   b , C 1122   b , C 1133   b  is calculated using the equation (2). C 2222   b , C 2233   b , C 3333   b  is calculated using the equation (3). As known from the symmetry of the stiffness matrix, C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b , C 2211   b , C 3311   b , C 3322   b  are checked at the same time, that is, calculating the six material parameters and checking the three material parameters. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                     
                 
                  
                 
                   C 
                   
                     11 
                      
                     
                         
                     
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                       
                   
                    
                   
                     C 
                     
                       22 
                        
                       
                           
                       
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                       
                   
                    
                   
                     C 
                     
                       33 
                        
                       
                           
                       
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                   
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                     
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0013]    (1.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (2) and the equation (3) are expressed as 
         [0000]      σ 11 +σ 11   H −α 11   P=C   11jj   b ε jj    (4)
 
         [0000]      σ ii   H −α ii   P=C    iijj   b ε jj   , i ∈ (2,3),  j ∈(1,3)   (5)
 
         [0014]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ ii , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b  (Refer to ab in the FIGURE). The strain ε ii  is then sprung back from ε ii   a  to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =−Δε ii . The equation of increments for the equation (4) and the equation (5) are: 
         [0000]      α ii   ΔP=C   iijj   b Δε jj (α ii   , i ∈(1,3))   (6)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0015]    Three Biot coefficients α 11 , α 22 , α 33  can be calculated from the equation (6). 
         [0016]    The method of measurement of stress and strain whole process material parameter by using method for unloading hydrostatic pressure of the present invention includes the following steps: 
         [0017]    (2.1) for false triaxial test which adopts 50mm×50mm×100 mm hexahedron sample, first applying a hydrostatic pressure σ 11 =σ 22 =σ 33 =σ H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0 ε jj   H   , i,j ∈(1,3)   (7)
 
         [0018]    where C iijj   0  is an initial stiffness matrix. 
         [0019]    (2.2) applying a deviator stress q, q=σ 11 +σ 11   H −σ 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (8)
 
         [0000]      σ ii   H   =C   iijj   b ε jj   , i ∈(2,3),  j ∈(1,3)   (9)
 
         [0020]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space. 
         [0021]    For the material parameters in stress state exceeding yield limit stress space, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 1111   b , C 1122   b , C 1133   b  is calculated using the equation (8). C 2222   b , C 2233   b  is calculated using the equation (9). Using the feature of the false triaxial test σ 22   H =σ 33   H , i.e. C 2211   b ε 11 , C 2222   b ε 22 +C2233 b ε 33 =C 3311   b ε 11 +C 3322   b ε 22 +C 3333   b ε 33 , C 3333   b  is calculated. As known from the symmetry of the stiffness matrix, C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b , C 2211   b , C 3311   b , C 3322   b  are checked at the same time, that is, calculating the six material parameters and checking the three material parameters. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                     
                 
                  
                 
                   C 
                   
                     11 
                      
                     
                         
                     
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                       
                   
                    
                   
                     C 
                     
                       22 
                        
                       
                           
                       
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                       
                   
                    
                   
                     C 
                     
                       33 
                        
                       
                           
                       
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0022]    (2.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (8) and the equation (9) are expressed as 
         [0000]      σ 11 +σ 11   H −α 11   P=C   11jj   b ε jj    (10)
 
         [0000]      σ ii   H −α ii   P=C   iijj   b ε jj   , i ∈(2,3), j ∈(1,3)   (11)
 
         [0023]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ ii , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b  (Refer to ab in the FIGURE). The strain ε ii  is then sprung back from ε ii   a , to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =−Δε ii . The equation of increments for the equation (10) and the equation (11) are: 
         [0000]      α ii   ΔP=C   iijj   b Δε jj (α ii   ,i  ∈(1,3))   (12)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0024]    Three Biot coefficients α 11 , α 22 , α 33  can be calculated. 
         [0025]    The method of measurement of stress and strain whole process material parameter by using method for unloading hydrostatic pressure of the present invention includes the following steps: 
         [0026]    (3.1) for traditional false triaxial test using ø50 mm×100 mm cylinder sample, first applying a hydrostatic pressure σ 11 =σ 22 =σ 33 =σ H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0 ε jj   H   , i,j ∈(1,3)   (13)
 
         [0027]    where C iijj   0  is an initial stiffness matrix. 
         [0028]    (3.2) applying a deviator stress q, q=σ 11 +σ 11   H −ε 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (14)
 
         [0000]      σ 22   H   =C   22jj   b ε jj , σ 22   H =σ 33   H    (15)
 
         [0029]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space; 
         [0030]    For the material parameters in stress state exceeding yield limit stress, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 111   b , C 1122   b  is calculated using the equation (14). C 2222   b  is calculated using the equation (15). Using the symmetry of the stiffness matrix C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b  and the feature of deformation measurement of the traditional false triaxial test, then σ 22   H =σ 33   H  and ε 22 =ε 33 , C 2211   b =C 3311   b , i.e. C 2211   b ε 11 +C 2222   b ε 22 +C 2233   b ε 33 +C 3311   b ε 11 +C 3322   b ε 22 +C 3333   b ε 33 , then C 2222   b =C 3333   b , that is, the current radial deformation measurement by ring-shaped transducer in the traditional false triaxial test substantially assumes that the sample failure occurs symmetrically. The three material parameters can then be calculated. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   
                     11 
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       22 
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       33 
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0031]    (3.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (14) and the equation (15) are expressed as 
         [0000]      σ 11 +σ 11   H α 11   P=C   11jj   b ε jj    (16)
 
         [0000]      σ ii   H −α ii   P=C   iijj   b ε jj   , i ∈(2,3),  j ∈(1,3)   (17)
 
         [0032]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ 11 , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b . The strain ε ii  is then sprung back from ε ii   a , to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =Δε ii . The equation of increments for the equation (16) and the equation (17) are: 
         [0000]      α ii   P=C   iijj   b Δε jj (α ii   ,i ∈(1,2))   (18)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0033]    Two Biot coefficients α 11 , α 22  can be calculated. 
         [0034]    The peak and proportional limit yield stress spaces decrease with an increase in damage until connecting with the residual strength yield surface. A specific expression is f yield (σ yield )f D (D)=Const, where f yield (σ yield ) is the yield stress space, f D (D) is a function of damage variable (D), Const is a constant, that is a product of the yield stress space and the function of damage variable is constant. 
         [0035]    The advantages of the method of the present invention are: 
         [0036]    The present invention is directed to the cyclic test of loading and unloading. With the assumption that only the deviator stress generates damage to the sample, a test method of unloading the hydrostatic pressure is invented in order to calculate mechanical parameters in different stages of stress and strain. Nine mechanical parameters can be calculated in connection with hexahedral (50 mm×50 mm×100 mm) pores connecting samples in the true triaxial test. Six mechanical parameters can be calculated for non-pores connecting samples. Nine mechanical parameters can be calculated in connection with hexahedral (50 mm×50 mm×100 mm) pores connecting samples in the traditional triaxial test. Six mechanical parameters can be calculated for non-pores connecting samples. Six mechanical parameters can be calculated in connection with cylinder (ø50 mm×100 mm) pores connecting samples in the traditional triaxial test. Four mechanical parameters can be calculated for non-pores connecting samples. The specific expressions and test methods are provided. 
         [0037]    The present invention also points out that apart from a complete isotropic body and equal amounts loaded and unloaded in three directions at the same time, the slope of the average stress and volume strain relation curve is meaningless for geological materials with linear elasticity. In addition, the yield limit stress space changes with the variations of damage and the expressions are provided. It is inappropriate for the residual traction stress not to be considered for the simulation of cyclic mechanical behavior. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0038]    The sole FIGURE is a loading and unloading diagram of parameters of the present invention. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     The First Embodiment 
       [0039]    The method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading of the present invention includes the following steps: 
         [0040]    (1.1) applying a hydrostatic pressure first in a true triaxial test, assuming σ 11 =σ 11   H , σ 22 =σ 22   H , σ 33 =σ 33   H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0   , i, j ∈(1,3)   (1)
 
         [0041]    where C iijj   0  is an initial stiffness matrix. 
         [0042]    (1.2) applying a deviator stress q, q=σ 11 +σ 11   H −σ 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (2)
 
         [0000]      σ ii   H   =C   iijj   b ε jj   , i ∈(2,3),  j ∈(1,3)   (3)
 
         [0043]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space. 
         [0044]    For the material parameters in stress state exceeding yield limit stress space, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 1111   b , C 1122   b , C 1133   b  is calculated using the equation (2). C 2222   b , C 2233   b , C 3333   b  is calculated using the equation (3). As known from the symmetry of the stiffness matrix, C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b , C 2211   b , C 3311   b , C 3322   b  are checked at the same time, that is, calculating the six material parameters and checking the three material parameters. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   
                     11 
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       22 
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       33 
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0045]    (1.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (2) and the equation (3) are expressed as 
         [0000]      σ 11 +σ 11   H −α 11   P=C   11jj   b ε jj    (4)
 
         [0000]      σ ii   H −α ii   P=C   iijj   b ε jj   , i ∈(2,3), j ∈(1,3)   (5)
 
         [0046]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ ii , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b  (Refer to ab in the FIGURE). The strain ε ii  is then sprung back from ε ii   a  to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =−Δε ii . The equation of increments for the equation (4) and the equation (5) are: 
         [0000]      α ii   ΔP=C   iijj   b Δε jj (α ii   , i ∈(1,3))   (6)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0047]    Three Biot coefficients α 11 , α 22 , α 33  can be calculated from the equation (6). 
       The Second Embodiment 
       [0048]    The method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading of the present invention includes the following steps: 
         [0049]    (2.1) for false triaxial test which adopts 50 mm×50 mm×100 mm hexahedron sample, first applying a hydrostatic pressure σ 11 =σ 22 =σ 33 =σ H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0 ε jj   H   , i,j ∈(1,3)   ( 7 )
 
         [0050]    where C iijj   0  is an initial stiffness matrix. 
         [0051]    (2.2) applying a deviator stress q, q=σ 11 +σ 11   H −σ 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (8)
 
         [0000]      σ ii   H   =C   iijj   b ε jj   , i ∈(2,3),  j ∈(1,3)   (9)
 
         [0052]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space. 
         [0053]    For the material parameters in stress state exceeding yield limit stress space, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 1111   b , C 1122   b , C 1133   b  is calculated using the equation (8). C 2222   b , C 2233   b  is calculated using the equation (9). Using the feature of the false triaxial test σ 22   H =σ 33   H , i.e. C 2211   b ε 11 +C 2222   b ε 22 +C 2233   b ε 33 =C 3311   b ε 11 +C 3322   b ε 22 +C 3333   b ε 33 , C 3333   b  is calculated. As known from the symmetry of the stiffness matrix, C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b , C 2211   b , C 3311   b , C 3322   b  are checked at the same time, that is, calculating the six material parameters and checking the three material parameters. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   
                     11 
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       22 
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       33 
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0054]    (2.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (8) and the equation (9) are expressed as 
         [0000]      σ 11 +σ 11   H −α 11   P=C    11jj   b ε jj    (10)
 
         [0000]      σ ii   H −α ii   P=C   iijj   b ε jj   , i ∈(2,3), j ∈(1,3)   (11)
 
         [0055]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ ii , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b  (Refer to ab in the FIGURE). The strain ε ii  is then sprung back from ε ii   a  to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =−Δε ii . The equation of increments for the equation (10) and the equation (11) are: 
         [0000]      α ii   ΔP=C   iijj   b Δε jj (α ii   ,i ∈(1,3))   (12)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0056]    Three Biot coefficients α 11 , α 22 , α 33  can be calculated. 
       The Third Embodiment 
       [0057]    The method of measurement of stress and strain whole process material parameter by using method for hydrostatic pressure unloading of the present invention includes the following steps: 
         [0058]    (3.1) for traditional false triaxial test using ø50 mm×100 mm cylinder sample, first applying a hydrostatic pressure σ 11 =σ 22 =σ 33 =σ H , a relation between the hydrostatic pressure and an initial strain ε ii   H  is: 
         [0000]      σ ii   H   =C   iijj   0 ε jj   H   , i,j ∈(1,3)   (13)
 
         [0059]    where C iijj   0  is an initial stiffness matrix. 
         [0060]    (3.2) applying a deviator stress q, q=σ 11 +σ 11   H −σ 11   H ; when the applied deviator stress is greater than the proportional limit stress q Yield  (i.e. yield limit stress space), linear stress-strain relations are expressed as: 
         [0000]      σ 11 +σ 11   H   =C   11jj   b ε jj    (14)
 
         [0000]      σ 22   H =C 22jj   b ε jj , σ 22   H =σ 33   H    (15)
 
         [0061]    where C iijj   b  is a stiffness matrix after exceeding yield limit stress space; 
         [0062]    For the material parameters in stress state exceeding yield limit stress space, magnitudes thereof are calculated according to the linear segment of unloading curve and unloading is carried out until the hydrostatic pressure tends to zero (Refer to def in the FIGURE). C 1111   b , C 1122   b  is calculated using the equation (14). C 2222   b  is calculated using the equation (15). Using the symmetry of the stiffness matrix C 2211   b =C 1122   b , C 2233   b =C 3322   b , C 3311   b =C 1133   b  and the feature of deformation measurement of the traditional false triaxial test, then σ 22   H =σ 33   H  and ε 22 =ε 33 , C 2211   b =C 3311   b , i.e. C 2211   b ε 11 +C 2222   b ε 22 +C 2233   b ε 33 =C 3311   b ε 11 +C 3322   b ε 22 +C 3333   b ε 33 , then C 2222   b =C 3333   b , that is, the current radial deformation measurement by ring-shaped transducer in the traditional false triaxial test substantially assumes that the sample failure occurs symmetrically. The three material parameters can then be calculated. When the material is completely isotropic, 
         [0000]    
       
         
           
             
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   
                     11 
                      
                     jj 
                   
                   b 
                 
               
               = 
               
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       22 
                        
                       jj 
                     
                     b 
                   
                 
                 = 
                 
                   
                     ∑ 
                     
                       j 
                       = 
                       1 
                     
                     3 
                   
                    
                   
                     C 
                     
                       33 
                        
                       jj 
                     
                     b 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the volume modulus C V , 
         [0000]    
       
         
           
             
               C 
               V 
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   1 
                 
                 3 
               
                
               
                 
                   ∑ 
                   
                     j 
                     = 
                     1 
                   
                   3 
                 
                  
                 
                   C 
                   iijj 
                   b 
                 
               
             
           
         
       
     
         [0000]    can be calculated; or when equal amounts are unloaded in three directions at the same time, 1/C V  can be calculated. 
         [0063]    (3.3) closing a valve and carrying out non-drainage test for pores connecting material after the application of the hydrostatic pressure is completed under a condition of saturation. Assuming that Bishop effective stress exists, then the equation (14) and the equation (15) are expressed as 
         [0000]      σ 11 +σ 11   H −α 11   P=C   11jj   b ε jj    (16)
 
         [0000]      σ ii   H −α ii   P=C   iijj   b ε jj   , i ∈(2,3), j ∈(1,3)   (17)
 
         [0064]    Under the condition of saturation in which material stiffness parameters C iijj   b  are obtained and the external applied stress σ ii , i ∈(1,3) is kept constant, the water pressure of non-draining test is unloaded from P a  to P b . The strain ε ii  is then sprung back from ε ii   a  to ε ii   b . The corresponding amount of deformation spring back is ε ii   a −ε ii   b =−Δε ii . The equation of increments for the equation (16) and the equation (17) are: 
         [0000]      α ii   ΔP=C   iijj   b Δε jj (α ii   ,i ∈(1,2))   (18)
 
         [0000]    
       
      
       P 
       a 
       −P 
       b 
       =ΔP  
      
     
         [0065]    Two Biot coefficients α 11 , α 22  can be calculated. 
         [0066]    The reasons for providing the method of the present invention are: 
         [0067]    (1) due to the elastic behavior which damages and weakens the material. The material generates new damage after the loading stress is greater than the proportional limit stress space. That is, the proportional limit stress space decreases with an increase in damage until connecting with the residual strength stress space. 
         [0068]    (2) because the sample loading is greater than the proportional limit stress space. The residual irreversible deformation and elastic deformation exist in the sample. This residual deformation must generate tensile stress in the sample. During loading-unloading-reloading cycle, a part of deformation generated by the compressive stress must fill up with the deformation generated by the residual tensile stress, which results in the reloading curve difficult to reflect the characteristics of the material. Therefore, it is inappropriate for the residual traction stress not to be considered for the simulation of cyclic mechanical behavior.