# Patent

## Title
Method for determining mechanical properties of a material

## Abstract

A tool having at least one vibration sensor and at least one standoff is disposed in a contact with a material and the at least one standoff of the tool is pushed into the material. Vibration is excited by the at least one vibration source and at least one coupling frequency of the tool is measured by the at least one vibration sensor. Based on the determined coupling frequency determining a contact stiffness of the at least one standoff and the mechanical properties of the material are determined taking into account mechanical properties of the at least one standoff.

## Background
FIELD OF THE INVENTION
This invention relates to methods for determination of mechanical properties of materials, namely Young modulus and Poisson ratio, and can be used, for example, for studying a formation surrounding a borehole, a ground floor or any surface for which material properties are to be measured.
BACKGROUND OF THE INVENTION
Many none destructive methods for determining mechanical properties of materials are known, such as wave propagation methods and micro/nano indentation methods (see, for example, General relationship between contact stiffness, contact depth, and mechanical properties for indentation in linear visco elastic solids using axisymetric indenters of arbitrary profile, Yang-Tse Chenga, Che-Min Cheng. 2005, APPLIED PHYSICS LETTERS 87, 111914, or An analysis of nanoindentation in linearly elastic solids. B. Poon a, D. Rittel b, G. Ravichandran. s.l.: International Journal of Solids and Structures 45, 2008, pp. 6018-6033). The known method provides for measuring hardness based on loading and unloading of specimens and can be accomplished only in laboratories.
The proposed method provides for easy and fast determination of mechanical properties of materials which requires only information on dynamic properties of a tool being used and can be done in situ.
SUMMARY OF THE INVENTION
The method for determining mechanical properties of a material comprises disposing a tool having at least one vibration sensor and at least one standoff in a contact with a material. Then the at least one standoff of the tool is pushed into the material and vibration is excited by at least one vibration source. At least one coupling frequency of the tool is measured by the at least one vibration sensor and a contact stiffness of the at least one standoff is determined based on the determined coupling frequency. The mechanical properties of the material are determined based on the determined contact stiffness of the at least one standoff and on dynamic properties of the at least one standoff.
The vibration source can be disposed inside or outside the tool.
The material can be a formation surrounding a borehole.
The tool can be clamped to the material.
The mechanical properties of the material and of the at least one standoff are Poisson ratio and Young modulus of the material and of the at least one standoff.



BRIEF DESCRIPTION OF THE DRAWINGS
The invention is explained by drawings where FIG. 1(a) shows an example of a tool with three standoffs, FIG. 2 shows a dynamic diagram for the tool showed in FIG. 1, FIG. 3 shows two variants of possible shapes of the standoffs.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
On FIG. 1 it is shown an example of a tool according to the proposed invention. The tool has standoffs (contact points) 1 and holes 2 for vibration sensors. The vibration sensors can be geophones or accelerometers that can measure direction of vibration. At least one vibration source also can be disposed inside the tool. The vibration source can be a shaker type or any vibration device capable of generate controlled vibration in a defined wave band. All those vibration sensors and vibration sources can be located in different positions in the tool and with different orientation. This is with the purpose of having the possibility of excitation of different vibration modes. It is especially important for addressing anisotropy materials and for quality control of measurements.
Then the at least one standoff of the tool is pushed into the material and vibration is excited by the vibration source. Vibration can be excited by any external or internal vibration source. At least one coupling frequency (a frequency at which the tool starts to vibrate due to the contact (coupling)) of the tool is measured by the at least one vibration sensor disposed inside the tool. The coupling frequencies are measured by examining a spectrum of the acquired readings of the vibration sensors.
A clamping force and a mass and moment of inertia of the tool are known, a shape and mechanical properties of the standoffs are also known. A combination of dynamic equations and equations for contact can be solved to find the unknowns—Young modulus and Poisson ratio of the material where the tool is located (and if necessary clamped). For the case of an isotropy material, the configuration can be freely chosen and can be used two or more standoffs for quality control. Knowing two first coupling frequencies of a particular configuration of the tool allows to calculate the Young modulus and Poisson coefficient of the material where the tool is located. Even with only the first coupling frequency registered for any configuration, it is possible to have a good approximation of the Young modulus regarding a good initial estimation of Poisson ratio. For the case of anisotropic material the larger the number of the frequencies, the more is the number of material properties that can be determined.
Below it is described a procedure for obtaining the equations for the case of using a tool with three standoffs shown on FIG. 1.
Equations of motion give the relationship to calculate coupling frequencies (Lagrange equation of the system) for the tool. We obtain:




  
    
      
        
          
            Ω
            
              coup
              ⁢
              
                
              
              ⁢
              1
            
          
          =
          
            
              
                
                  k
                  1
                
                +
                
                  k
                  2
                
                +
                
                  k
                  3
                
              
              m
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            Ω
            
              coup
              ⁢
              
                
              
              ⁢
              2
            
          
          =
          
            
              
                
                  
                    k
                    2
                  
                  ⁢
                  
                    d
                    2
                  
                
                +
                
                  
                    k
                    3
                  
                  ⁢
                  
                    d
                    
                      
                        
                      
                      ⁢
                      2
                    
                  
                
              
              
                I
                2
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            Ω
            
              coup
              ⁢
              
                
              
              ⁢
              3
            
          
          =
          
            
              
                
                  
                    k
                    2
                  
                  ⁢
                  
                    c
                    2
                  
                
                +
                
                  
                    k
                    3
                  
                  ⁢
                  
                    c
                    2
                  
                
                +
                
                  
                    k
                    1
                  
                  ⁢
                  
                    a
                    
                      
                        
                      
                      ⁢
                      2
                    
                  
                
              
              
                I
                3
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




where Ωcoup1, Ωcoup2, Ωcoup3 are the coupling frequencies of the tool, k1 k2, k3-contact stiffness of three standoffs. Based on specific design it is possible to achieve very similar values of the contact stiffness so we can assume them equal: k1=k2=k3. Parameters d, c and a are distances from a center of mass to the location of the standoffs (see FIG. 2); m is a mass of the tool, and I2, I3 are moments of inertia of the tool body/package respect to an axis of rotation. Making several measurements of the frequencies, a more robust value of the frequencies can be obtained.
Below, based on contact theory several relationships for the contact stiffness are derived, the main relationship is




  k
  =
  
    2
    ⁢
    
      
    
    ⁢
    
      E
      r
    
    ⁢
    
      
        A
      
      
        π
      
    
  




where k is a contact stiffness and A is a projected contact area.
Er is a reduced stiffness modulus that is defined as




  
    1
    
      E
      r
    
  
  =
  
    
      
        (
        
          1
          -
          υ
        
        )
      
      E
    
    +
    
      
        (
        
          1
          -
          
            υ
            ′
          
        
        )
      
      
        E
        ′
      
    
  




Here, υ, E are Poisson ratio and Young modulus of the material to be tested respectively. υ′, E′ are Poisson ratio and Young modulus of the standoffs.
The definition of the contact area A depends on a shape of the standoff. It is calculated using Hertz contact theory.
Bellow there are examples of calculation for two shapes of the standoffs (FIGS. 3a and 3b).
Semi Spherical Stand Off
Using Hertz contact theory, is possible to relate a clamping force of the tool to the contact area. For a shape as indicated in the FIG. 3a, an area of contact will be

A2=π·R·h 

where R is a radius of a round standoff, h is a height of penetration of the standoff into the material. A force is related to the reduced modulus and area of contact by




  F
  =
  
    
      A
      ·
      
        E
        r
      
      ·
      
        R
        
          1
          2
        
      
      ·
      
        h
        
          3
          2
        
      
    
    3
  




From this the following relationship is obtained




  
    A
  
  =
  
    
      (
      
        
          3
          ·
          F
          ·
          R
        
        
          4
          ·
          
            E
            r
          
        
      
      )
    
    
      1
      6
    
  




Using this the relationship for the stiffness contact for the round standoff (with υ, E) pushed with a force F into a sample (with υ′, E′) is obtained




  K
  =
  
    
      
        2
        ·
        
          E
          r
        
      
      
        π
      
    
    ⁢
    
      
        (
        
          
            3
            ·
            F
            ·
            R
          
          
            4
            ·
            
              E
              r
            
          
        
        )
      
      
        1
        6
      
    
  




The obtained value of the contact stiffness can be introduced in the formulas of the frequencies (1):




  {
  
    
      
        
          
            
              
                Ω
                
                  coup
                  1
                
              
              =
              
                
                  
                    α
                    m
                  
                
                =
                
                  
                    
                      
                        2
                        ⁢
                        
                          k
                          1
                        
                      
                      +
                      
                        2
                        ⁢
                        
                          k
                          2
                        
                      
                      +
                      
                        2
                        ⁢
                        
                          k
                          3
                        
                      
                    
                    m
                  
                
              
            
          
        
        
          
            
              
                Ω
                
                  coup
                  2
                
              
              =
              
                
                  
                    ξ
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
                =
                
                  
                    
                      
                        
                          k
                          2
                        
                        ·
                        
                          d
                          2
                        
                      
                      +
                      
                        
                          k
                          3
                        
                        ·
                        
                          d
                          2
                        
                      
                    
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
              
            
          
        
        
          
            
              
                Ω
                
                  coup
                  3
                
              
              =
              
                
                  
                    η
                    
                      m
                      ·
                      
                        r
                        3
                        2
                      
                    
                  
                
                =
                
                  
                    
                      
                        
                          k
                          2
                        
                        ·
                        
                          c
                          2
                        
                      
                      +
                      
                        
                          k
                          3
                        
                        ·
                        
                          c
                          2
                        
                      
                      +
                      
                        
                          k
                          1
                        
                        ·
                        
                          a
                          2
                        
                      
                    
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
              
            
          
        
      
      ⁢
      
        

      
      ⁢
      with
      ⁢
      
        

      
      ⁢
      
        K
        
          1
          ,
          2
          ,
          3
        
      
    
    =
    
      
        
          2
          ·
          
            E
            r
          
        
        
          π
        
      
      ⁢
      
        
          (
          
            
              3
              ·
              F
              ·
              R
            
            
              4
              ·
              
                E
                r
              
            
          
          )
        
        
          1
          3
        
      
    
  




Er is found and then Young modulus and Poisson ratio of the tested material are determined.
Conical Stand Offs (Spikes)
As for the previous case Hertz contact theory is used to relate a clamping force to a contact area. For a shape as indicated in FIG. 3b, the area of contact will be




  
    A
    2
  
  =
  
    
      
        
          h
          2
        
        ·
        tan
      
      ⁢
      
        
      
      ⁢
      
        
          (
          
            
              π
              2
            
            -
            Θ
          
          )
        
        2
      
    
    π
  




Where θ is an angle defining a sharpness of the standoffs, h-a depth of penetration into the material. The force is related to the reduced modulus and area of contact by




  F
  =
  
    
      A
      ·
      
        E
        r
      
    
    
      
        2
        ·
        tan
      
      ⁢
      
        
      
      ⁢
      
        (
        
          
            π
            2
          
          -
          θ
        
        )
      
    
  




From this the relationship is obtained




  
    A
  
  =
  
    
      (
      
        
          
            F
            ·
            tan
          
          ⁢
          
            
          
          ⁢
          
            (
            
              
                π
                2
              
              -
              θ
            
            )
          
        
        
          2
          ·
          
            E
            r
          
        
      
      )
    
    
      1
      2
    
  




Using this the relationship for the stiffness contact for a conical standoff (with υ, E) pushed with a force F into a casing (with υ′, E′) is obtained




  K
  =
  
    
      
        2
        ·
        
          E
          r
        
      
      
        π
      
    
    ⁢
    
      
        (
        
          
            
              2
              ·
              F
              ·
              tan
            
            ⁢
            
              
            
            ⁢
            
              (
              
                
                  π
                  2
                
                -
                θ
              
              )
            
          
          
            E
            r
          
        
        )
      
      
        1
        2
      
    
  




The obtained value of the contact stiffness can be introduced in the formulas of the frequencies (1):




  {
  
    
      
        
          
            
              
                Ω
                
                  coup
                  1
                
              
              =
              
                
                  
                    α
                    m
                  
                
                =
                
                  
                    
                      
                        2
                        ⁢
                        
                          k
                          1
                        
                      
                      +
                      
                        2
                        ⁢
                        
                          k
                          2
                        
                      
                      +
                      
                        2
                        ⁢
                        
                          k
                          3
                        
                      
                    
                    m
                  
                
              
            
          
        
        
          
            
              
                Ω
                
                  coup
                  2
                
              
              =
              
                
                  
                    ξ
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
                =
                
                  
                    
                      
                        
                          k
                          2
                        
                        ·
                        
                          d
                          2
                        
                      
                      +
                      
                        
                          k
                          3
                        
                        ·
                        
                          d
                          2
                        
                      
                    
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
              
            
          
        
        
          
            
              
                Ω
                
                  coup
                  3
                
              
              =
              
                
                  
                    η
                    
                      m
                      ·
                      
                        r
                        3
                        2
                      
                    
                  
                
                =
                
                  
                    
                      
                        
                          k
                          2
                        
                        ·
                        
                          c
                          2
                        
                      
                      +
                      
                        
                          k
                          3
                        
                        ·
                        
                          c
                          2
                        
                      
                      +
                      
                        
                          k
                          1
                        
                        ·
                        
                          a
                          2
                        
                      
                    
                    
                      m
                      ·
                      
                        r
                        2
                        2
                      
                    
                  
                
              
            
          
        
      
      ⁢
      
        

      
      ⁢
      with
      ⁢
      
        

      
      ⁢
      
        K
        
          1
          ,
          2
          ,
          3
        
      
    
    =
    
      
        
          2
          ·
          
            E
            r
          
        
        
          π
        
      
      ⁢
      
        
          (
          
            
              2
              ·
              F
              ·
              
                tan
                ⁡
                
                  (
                  
                    
                      π
                      2
                    
                    -
                    θ
                  
                  )
                
              
            
            
              E
              r
            
          
          )
        
        
          1
          2
        
      
    
  




Then it is solved for Er and then Young modulus and Poisson ratio of the tested material are determined.
The described procedure is general, for any shape or mechanical design of the tool and any shape and material of the standoffs.

## Claims

1. A method for determining mechanical properties of a material comprising:
disposing a tool having at least one vibration sensor and at least one standoff in a contact with a material,
pushing the at least one standoff of the tool into the material,
exciting vibration by the at least one vibration source,
measuring at least one coupling frequency of the tool by the at least one vibration sensor,
determining a contact stiffness of the at least one standoff based on the determined coupling frequency and determining the mechanical properties of the material based on the determined contact stiffness of the at least one standoff and on mechanical properties of the at least one standoff.

2. The method of claim 1 wherein the at least one vibration source is disposed inside the tool.

3. The method of claim 1 wherein the at least one vibration source is disposed outside the tool.

4. The method of claim 1 wherein the mechanical properties of the material and of the at least one standoff are Poisson ratio and Young modulus of the material and of the at least one standoff.

5. The method of claim 1 wherein material is a formation surrounding a borehole.

6. The method of claim 1 wherein the tool is clamped to the material.