**1. Introduction**

In this chapter, we introduce a new methodology to measure the stress-train relation for ductile materials subjected to a continuum axial load; this loading was carried out by means of a universal testing machine. To obtain the stress-strain relation, materials were submitted to standardized compression tests, and to develop the technique, we used a laser beam impacting the cross section of the materials; the laser reflection was studied through Gaussian beam propagation equations and digital image processing. The scope of the present research is to determine a relation between the free propagation of a Gaussian beam and the stress-strain evolution of a ductile material.

Stress-strain diagrams are very important in understanding the behavior of materials under different loads [1]; these diagrams show the elastic, plastic, and rupture characteristics of materials. There are two methods to obtain these diagrams.

First, the invasive methods from which mechanics draw by doing physical tests, such as the test tube in which a standardized probe is placed on a testing machine; a continuum load is applied to it, and the resulting deformation is measured [2]. Other invasive methods are optical methods, which can determine residual stress, in-field displacements, and strain. Here, hole drilling is the most used technique.

Originally developed in 1930 by Mathar [3], it is nowadays a standardized technique by the ASTM [4].

function represented from Gaussian beams can be deduced by applying boundary conditions in the optical resonator, where the laser radiation is produced; this

½ � *ω*ð Þ*z*

where *E*<sup>0</sup> is the electric field amplitude, *ω*<sup>0</sup> is the beam waist, *ω*(*z*) is how the beam propagates, k is the wave number (*k* ¼ 2*π=λ*), *λ* is the wavelength, *R*(*z*) is the curvature radius of the spherical waves, and *ƞ*(*z*) is the beam phase angle [19]. Laser beams are able to pass through different media; the light reflection occurs when it arrives to the boundary separating two media of different optical densities, and some of the energy is reflected back into the first medium [20]; taking this outset, if a laser beam strikes a mirror, the reflection can be studied as a Gaussian propagation. In our case, the metallic surface will be modeled as a convex mirror. As it is well-known, there is a relation between the focal length and the curvature radius of

<sup>2</sup> � *kx*<sup>2</sup>

*kx*<sup>2</sup> *<sup>z</sup>* <sup>1</sup> <sup>þ</sup> *<sup>z</sup>*0<sup>2</sup> *z*2 � � þ

where *z*<sup>0</sup> is the initial Rayleigh distance and *f* is the focal length of the mirror [21]. To study the propagation through a lens phenomenon, the Gaussian propagation equations are established in which *ω*<sup>0</sup> and *z*<sup>0</sup> turns into *ω*<sup>1</sup> and *z*1, respectively,

*<sup>ω</sup>*<sup>1</sup> <sup>¼</sup> *<sup>ω</sup>*<sup>0</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

*<sup>z</sup>*<sup>1</sup> <sup>¼</sup> *<sup>f</sup>*

to calculate the new beam waist in the propagation axis, we have

*ω zp* � � <sup>¼</sup> *<sup>ω</sup>*<sup>1</sup>

compression test. Substituting Eqs. (3), (4) in (5) it can be deduced:

*<sup>z</sup>*0<sup>4</sup> � *<sup>x</sup>*<sup>2</sup> <sup>þ</sup> *<sup>x</sup>* <sup>1</sup> � *<sup>ω</sup> zp*

*zp* 2

of the focal length is determined applying

**31**

<sup>1</sup> <sup>þ</sup> *<sup>z</sup>*0<sup>2</sup> *f* 2

<sup>1</sup> <sup>þ</sup> *<sup>f</sup>* 2 *z*0<sup>2</sup>

This pair of equations involves how Gaussian beam propagates [21], so in order

where *ω*(*zp*) is the new beam waist at *zp* which is the propagation distance. The goal of this research is to reach a general relation between the strains caused in any material by measuring with one laser beam, its reflection and relating it to a focal length transformation. As we implicate the focal length, we need these Gaussian propagation equations to calculate the changes in the focal length during the

> � �<sup>2</sup> *ω*0<sup>2</sup>

The parameter *zp* is calculated doing DIP, taking each area increment per second during the compression test. Once the result from Eq. (6) is obtained, the variation

� *z*<sup>0</sup>

!

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þ *zp* 2 *z*1 2

" #

<sup>2</sup>*R z*ð Þ � *kz* <sup>þ</sup> *<sup>η</sup>*ð Þ*<sup>z</sup>*

*kx*<sup>2</sup>

r� � *,* (3)

� � *:* (4)

s� �*,* (5)

<sup>2</sup> <sup>¼</sup> <sup>0</sup>*:* (6)

*<sup>f</sup> ,* (2)

*,* (1)

amplitude is described by

*E x*ð Þ¼ *; z E*<sup>0</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

*Digital Image Processing Applied to Optical Measurements*

after a distance *f* and they are calculated by

*ω*0

*<sup>ω</sup>*ð Þ*<sup>z</sup>* � *exp* �*x*<sup>2</sup>

a mirror. Using this relation the phase of the transmitted wave is altered to

*φ*ð Þ¼ *x; z* 2*η*ð Þ� *z*

Second, nondestructive or noninvasive optical methods such as multiple laser displacement sensors which are applied to piping systems [5], stress measurements using crystal curvature technique during film growths [6], and another technique called deflectometry [7, 8], in which light passes through arrays that measure the curvature of objects; they are proposed as mirrors. The reflection caused by these objects is seen in a CCD camera, and the pattern is studied with standard phase shift techniques; most of the used objects are aspherical.

An image was defined as a two-dimensional function: *f x*ð Þ *; y* where *x* and *y* are plane coordinates and the amplitude *f* at any pair of coordinates ð Þ *x; y* is called the intensity referred as a gray level of the image at that point [9]. One of the first applications of digital images was in the newspaper industry through the introduction of the Bartlane cable picture transmission system in the early 1920s, where specialized printing equipment coded pictures for cable transmission and then reconstructed them at the receiving end [9]. However, the image processing method itself required storage and computational power; nevertheless, by this time, it was not possible to talk about processing since the development of the computer hadn't arrived. One of the first works in this field was done by Harry Andrews [10], whom introduced techniques in order to work with images as matrixes; another book about the treatment of pictures was introduced by 1976 [11].

Digital image processing has been applied widely, in such cases as in control level for industries [12], which is based on the treatment of a special light shape generated by a laser, also measuring displacements in infrastructures by a matrix of points [13]. Some researchers around the application in engineering are evaluating the curvature of a pipeline using a parabolic function to determine its deformation [14], analyzing the location of cracks in reinforced concrete by monitoring the surface deformation with two cameras [15], strain measurements during film stretching process following the digital changes of the marker [16], and measuring the strain in time for standardized fatigue tests [17], among others.

In this chapter we introduce a novel methodology to measure the stress-strain behavior of ductile materials; this technique uses a laser beam focused on the cross section of a material probe which is under a compression test. We propose that the materials (1018 steel and aluminum alloy) will act as optical spherical mirrors. This is due to their surface transformation; the laser focuses in a surface which is completely flat, and it will deform as a reaction of the compression test. Therefore, the laser strikes this surface, and acting as a "mirror," it will reflect and scatter the beam. Thus, the reflected laser area will increase as the deformation increases. This laser reflection is analyzed with Gaussian beam propagation equations by using digital image processing to measure each increasing area. This is the procedure we follow to obtain a relation between the laser-beam-free propagation and the strain, which we propose to be similar to the cross-section deformation of both materials. We show the accuracy, error, and sensitivity of the methodology, and we theoretically demonstrate how the laser beam spreads; the present research can be a suitable technique as it only uses a single laser beam for stress-strain characterization.

### **2. Theory**

#### **2.1 Gaussian beam analysis**

The Gaussian beam is the simplest, the most known, and the most worked with, because its characteristics and evolution are well-known [18]. The amplitude *Digital Image Processing Applied to Optical Measurements DOI: http://dx.doi.org/10.5772/intechopen.88704*

Originally developed in 1930 by Mathar [3], it is nowadays a standardized tech-

Second, nondestructive or noninvasive optical methods such as multiple laser displacement sensors which are applied to piping systems [5], stress measurements using crystal curvature technique during film growths [6], and another technique called deflectometry [7, 8], in which light passes through arrays that measure the curvature of objects; they are proposed as mirrors. The reflection caused by these objects is seen in a CCD camera, and the pattern is studied with standard phase shift

An image was defined as a two-dimensional function: *f x*ð Þ *; y* where *x* and *y* are plane coordinates and the amplitude *f* at any pair of coordinates ð Þ *x; y* is called the intensity referred as a gray level of the image at that point [9]. One of the first applications of digital images was in the newspaper industry through the introduction of the Bartlane cable picture transmission system in the early 1920s, where specialized printing equipment coded pictures for cable transmission and then reconstructed them at the receiving end [9]. However, the image processing method itself required storage and computational power; nevertheless, by this time, it was not possible to talk about processing since the development of the computer hadn't arrived. One of the first works in this field was done by Harry Andrews [10], whom introduced techniques in order to work with images as matrixes; another

Digital image processing has been applied widely, in such cases as in control level for industries [12], which is based on the treatment of a special light shape generated by a laser, also measuring displacements in infrastructures by a matrix of points [13]. Some researchers around the application in engineering are evaluating the curvature of a pipeline using a parabolic function to determine its deformation [14], analyzing the location of cracks in reinforced concrete by monitoring the surface deformation with two cameras [15], strain measurements during film stretching process following the digital changes of the marker [16], and measuring

In this chapter we introduce a novel methodology to measure the stress-strain behavior of ductile materials; this technique uses a laser beam focused on the cross section of a material probe which is under a compression test. We propose that the materials (1018 steel and aluminum alloy) will act as optical spherical mirrors. This is due to their surface transformation; the laser focuses in a surface which is completely flat, and it will deform as a reaction of the compression test. Therefore, the laser strikes this surface, and acting as a "mirror," it will reflect and scatter the beam. Thus, the reflected laser area will increase as the deformation increases. This laser reflection is analyzed with Gaussian beam propagation equations by using digital image processing to measure each increasing area. This is the procedure we follow to obtain a relation between the laser-beam-free propagation and the strain, which we propose to be similar to the cross-section deformation of both materials. We show the accuracy, error, and sensitivity of the methodology, and we theoretically demonstrate how the laser beam spreads; the present research can be a suitable technique as it only uses a single laser beam for stress-strain characterization.

The Gaussian beam is the simplest, the most known, and the most worked with, because its characteristics and evolution are well-known [18]. The amplitude

nique by the ASTM [4].

*Digital Imaging*

**2. Theory**

**30**

**2.1 Gaussian beam analysis**

techniques; most of the used objects are aspherical.

book about the treatment of pictures was introduced by 1976 [11].

the strain in time for standardized fatigue tests [17], among others.

function represented from Gaussian beams can be deduced by applying boundary conditions in the optical resonator, where the laser radiation is produced; this amplitude is described by

$$E(\mathbf{x}, z) = E\_0 \frac{a \nu\_0}{o(z)} \cdot \exp\left[\frac{-\mathbf{x}^2}{[o(z)]^2} - \frac{k \mathbf{x}^2}{2R(z)} - kz + \eta(z)\right],\tag{1}$$

where *E*<sup>0</sup> is the electric field amplitude, *ω*<sup>0</sup> is the beam waist, *ω*(*z*) is how the beam propagates, k is the wave number (*k* ¼ 2*π=λ*), *λ* is the wavelength, *R*(*z*) is the curvature radius of the spherical waves, and *ƞ*(*z*) is the beam phase angle [19]. Laser beams are able to pass through different media; the light reflection occurs when it arrives to the boundary separating two media of different optical densities, and some of the energy is reflected back into the first medium [20]; taking this outset, if a laser beam strikes a mirror, the reflection can be studied as a Gaussian propagation. In our case, the metallic surface will be modeled as a convex mirror. As it is well-known, there is a relation between the focal length and the curvature radius of a mirror. Using this relation the phase of the transmitted wave is altered to

$$\rho(\mathbf{x}, z) = 2\eta(z) - \frac{k\mathbf{x}^2}{z\left(1 + \frac{x\_0}{x^2}\right)} + \frac{k\mathbf{x}^2}{f},\tag{2}$$

where *z*<sup>0</sup> is the initial Rayleigh distance and *f* is the focal length of the mirror [21]. To study the propagation through a lens phenomenon, the Gaussian propagation equations are established in which *ω*<sup>0</sup> and *z*<sup>0</sup> turns into *ω*<sup>1</sup> and *z*1, respectively, after a distance *f* and they are calculated by

$$\alpha\_1 = \frac{\alpha\_0}{\sqrt{\left(\mathbf{1} + \frac{x\_0}{f^2}\right)}},\tag{3}$$

$$z\_1 = \frac{f}{\left(1 + \frac{f^2}{x\_0^2}\right)}.\tag{4}$$

This pair of equations involves how Gaussian beam propagates [21], so in order to calculate the new beam waist in the propagation axis, we have

$$
\alpha(\mathbf{z}\_p) = \alpha\_1 \sqrt{\left(\mathbf{1} + \frac{\mathbf{z}\_p^{-2}}{\mathbf{z}\_1^2}\right)},\tag{5}
$$

where *ω*(*zp*) is the new beam waist at *zp* which is the propagation distance.

The goal of this research is to reach a general relation between the strains caused in any material by measuring with one laser beam, its reflection and relating it to a focal length transformation. As we implicate the focal length, we need these Gaussian propagation equations to calculate the changes in the focal length during the compression test. Substituting Eqs. (3), (4) in (5) it can be deduced:

$$\frac{{}^{2}\mathbf{z}\_{p}^{2}}{{}^{2}\mathbf{z}\_{0}{}^{4}}\cdot\mathbf{x}^{2} + \mathbf{x}\left(\mathbf{1} - \frac{{}^{g}\mathbf{z}\_{p}{}^{2}\right) {}^{2}\mathbf{z}\_{0}{}^{2} = \mathbf{0}.\tag{6}$$

The parameter *zp* is calculated doing DIP, taking each area increment per second during the compression test. Once the result from Eq. (6) is obtained, the variation of the focal length is determined applying

$$
\xi = \frac{f\_f - f\_i}{f\_i},
\tag{7}
$$

where *ff* is the final focal length and *fi* is the initial focal length; following this procedure, we obtained a dimensionless variable.

### **2.2 Compression tests**

Mechanical tests are used to know the mechanical properties of materials; a compression test enables the user to understand the behavior of a material under a continuum axial load; from this test it is possible to obtain the stress-strain diagram of any material [2]. In this chapter we worked with samples of steel with dimensions of 2.5 � 2.5 � 2.2 cm and aluminum samples of 1 � 1 � 0.9 cm according to the small-type standard dimensions [22]; all samples undergo into compression test according to ASTM E-9 [22]. The tests were performed with a speed ratio of 0.2 cm/s up to 90 MPa for steel and up to 20 MPa for aluminum. Since we are working on the elastic part of the diagram, therefore we can apply Hooke's law:

$$
\sigma = E \cdot \varepsilon,\tag{8}
$$

and

**Figure 1.**

where

<sup>A</sup> � 0.009 m<sup>2</sup>

**33**

**3. Methodology**

.

coordinates, we obtain

*Gaussian propagation: intensity profile.*

*Digital Image Processing Applied to Optical Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

*U x*ð Þ¼ *; <sup>y</sup> exp* � *<sup>x</sup>*2þ*y*<sup>2</sup> ½ � ð Þ *:* (12)

� �*,* (14)

*<sup>f</sup> ,* (15)

*,* (13)

Replacing Eqs. (11) and (12) into Eq. (10), expanding formally and changing

*ik* 2*f* þ*b*2 4*a* � ��*ρ*2 h i

*iπ <sup>a</sup><sup>λ</sup> <sup>f</sup> exp*

*<sup>a</sup>* <sup>¼</sup> <sup>1</sup> � *ik*

*<sup>b</sup>* <sup>¼</sup> *ik*

According to the present research, the initial area of the reflected laser was

As we established before, we use two ductile materials: steel and aluminum; in total, seven samples were hand-sanded and polished in the same manner to mirror grade in the laser-impact section, and they were placed in the machine perpendicular to the piston to begin the compression test. A diagram of this method is shown schematically in **Figure 2a**. The output of a He-Ne laser (L), with a wavelength

2*f*

*U*ð Þ¼ *ρ*

where *σ* is the engineering stress, *E* is the Young's module of the material (200 and 70 GPa) for steel and aluminum, respectively, and *ε* is the engineering strain. We propose that the experimental stress is going to behave similar to the Hooke's law; thus, we relate Eq. (8) to the focal length by

$$
\sigma = E \cdot K \cdot \xi \,\tag{9}
$$

where *K* is a material coefficient proposed in this chapter and *ξ* is the dimensionless value obtained in Eq. (7). The coefficient *K* is obtained from the relation between the slopes of both graphs of interest: stress-strain diagram and DIP plot.

#### **2.3 Theoretical analysis**

In this research, we used a He-Ne laser with specific parameters such as initial intensity, beam waist, and Rayleigh distance that are shown in its handbook [23].

In **Figure 1** we plot Eqs. (1) and (2) substituting the results from Eq. (7), the initial parameters of our laser and propagating formally this laser beam. In this figure, we can observe that when the focal length increases, the beam waist has also an increase, in this theoretical demonstration; we propose that these increments are proportional as we assume that *ε* ≈ *f*.

By using Fresnel diffraction and the treatment of lens transmission [24], we have

$$U\_f(\boldsymbol{u}, \boldsymbol{v}) = \frac{\exp\left[\frac{\bar{\boldsymbol{u}}\left(\boldsymbol{u}^2 + \boldsymbol{v}^2\right)}{\mathcal{T}}\right]}{i\lambda\!\!f} \, \* \, \left[\int\_{-\infty}^{\infty} U(\boldsymbol{x}, \boldsymbol{y}) \cdot \boldsymbol{t}\_A \cdot \exp\left[\frac{\bar{\boldsymbol{u}}\_{(\rm X^4 \boldsymbol{y} + \boldsymbol{y}\boldsymbol{v})}}{\mathcal{T}}\right] d\mathbf{x} d\mathbf{y},\tag{10}$$

where

$$t\_A = \exp\left[\overleftarrow{\gamma}^{\mathbb{A}}(\mathsf{x}^2 + \mathsf{y}^2)\right] \tag{11}$$

*Digital Image Processing Applied to Optical Measurements DOI: http://dx.doi.org/10.5772/intechopen.88704*

**Figure 1.** *Gaussian propagation: intensity profile.*

and

*<sup>ξ</sup>* <sup>¼</sup> *ff* � *fi fi*

where *ff* is the final focal length and *fi* is the initial focal length; following

Mechanical tests are used to know the mechanical properties of materials; a compression test enables the user to understand the behavior of a material under a continuum axial load; from this test it is possible to obtain the stress-strain diagram of any material [2]. In this chapter we worked with samples of steel with dimensions of 2.5 � 2.5 � 2.2 cm and aluminum samples of 1 � 1 � 0.9 cm according to the small-type standard dimensions [22]; all samples undergo into compression test according to ASTM E-9 [22]. The tests were performed with a speed ratio of 0.2 cm/s up to 90 MPa for steel and up to 20 MPa for aluminum. Since we are working on the elastic part of the diagram, therefore we can apply Hooke's law:

where *σ* is the engineering stress, *E* is the Young's module of the material (200 and 70 GPa) for steel and aluminum, respectively, and *ε* is the engineering strain. We propose that the experimental stress is going to behave similar to the Hooke's

where *K* is a material coefficient proposed in this chapter and *ξ* is the dimensionless value obtained in Eq. (7). The coefficient *K* is obtained from the relation between the slopes of both graphs of interest: stress-strain diagram and DIP plot.

In this research, we used a He-Ne laser with specific parameters such as initial intensity, beam waist, and Rayleigh distance that are shown in its handbook [23]. In **Figure 1** we plot Eqs. (1) and (2) substituting the results from Eq. (7), the initial parameters of our laser and propagating formally this laser beam. In this figure, we can observe that when the focal length increases, the beam waist has also an increase, in this theoretical demonstration; we propose that these increments are

By using Fresnel diffraction and the treatment of lens transmission [24],

*U x*ð Þ� *; y tA* � *exp*

*ik* <sup>2</sup>*<sup>f</sup> <sup>x</sup>*2þ*y*<sup>2</sup> ð Þ h i

*ik* 2*f* ð Þ *xu*þ*yv* h i

dxdy*,* (10)

(11)

ð ð∞

�∞

*tA* ¼ *exp*

this procedure, we obtained a dimensionless variable.

law; thus, we relate Eq. (8) to the focal length by

**2.2 Compression tests**

*Digital Imaging*

**2.3 Theoretical analysis**

we have

where

**32**

proportional as we assume that *ε* ≈ *f*.

*ik* <sup>2</sup>*<sup>f</sup> <sup>u</sup>*2þ*v*<sup>2</sup> ð Þ h i

*<sup>i</sup>λ<sup>f</sup>* <sup>∗</sup>

*Uf*ð Þ¼ *<sup>u</sup>; <sup>v</sup> exp*

*,* (7)

*σ* ¼ *E* � *ε,* (8)

*σ* ¼ *E* � *K* � *ξ,* (9)

$$U(\mathbf{x}, \mathbf{y}) = \exp^{-\left[\left(\mathbf{x}^2 + \mathbf{y}^2\right)\right]}.\tag{12}$$

Replacing Eqs. (11) and (12) into Eq. (10), expanding formally and changing coordinates, we obtain

$$U(\rho) = \frac{i\pi}{a\lambda f} \exp\left[\left(\frac{\bar{u} + \frac{\rho^2}{4a}}{2}\right) \cdot \rho^2\right],\tag{13}$$

where

$$a = \left(\mathbf{1} - \frac{ik}{\mathcal{Y}}\right),\tag{14}$$

$$b = \frac{ik}{f},\tag{15}$$

According to the present research, the initial area of the reflected laser was <sup>A</sup> � 0.009 m<sup>2</sup> .

### **3. Methodology**

As we established before, we use two ductile materials: steel and aluminum; in total, seven samples were hand-sanded and polished in the same manner to mirror grade in the laser-impact section, and they were placed in the machine perpendicular to the piston to begin the compression test. A diagram of this method is shown schematically in **Figure 2a**. The output of a He-Ne laser (L), with a wavelength

**Figure 2.** *Experimental setup: (a) material with no strain; (b) the strained material will scatter light.*

λ = 632 nm and a power of 1 mW, is placed 15 cm from the sample (Sa) in order to irradiate the polished face, and then, the reflected beam impacts a screen (Sc) which is placed 10 cm beside the laser beam and parallel to the sample. The material is first completely flat, with no stress, and as it is compressed with the universal machine (M), the size of the cross section will decrease as a reaction (**Figure 2b**). During the tests, the deformation rate was 0.2 cm/s with duration of 5 min approximately, and while they were taking place, the reflected light was recorded with a video camera (Vc) which is placed in front and parallel to the screen aside the sample. All tests were done at room temperature, and samples were compressed around 2 mm.

*im2bw* to do a binary code of the image. At this step we can apply the function *bwarea* to know the area of this "circle." We can observe in **Figure 3d** and **e** how the beam diameter increases, for a *ti* = 0 s and *tf* = Δt s. By following the methodology described and using Eq. (6), we can obtain the change in the focal length and relate

*Image processing: (a) how the image is loaded; (b) grayscale image; (c) binary image; (d) initial beam*

The code was written to repeat this process iteratively for every image. We propose that the calculated area is the parameter *ω*(*zp*). Therefore Eq. (6) is used to calculate the focal length. We know that as a result of doing a mechanical test, we can obtain a stress-strain diagram; it is calculated through strain gages and special software incorporated to the machine. At the end of the iterative process, we plot

Along this chapter, we have mentioned that we prepared four steel samples and three aluminum probes, to test them under standardized compression test; for each one we obtained two graphs of interest, the stress-strain diagram and the DIP plot; **Figure 4** shows these two graphs of interest. **Figures 4**–**8** correspond to steel

In **Figure 4a** we show the stress-strain diagram obtained from the software of the machine; in **Figure 4b** we can observe the related stress-focal length diagram obtained by following the experimental procedure and using Eq. (9). In the graphs we can observe a linear behavior. In section 2.2, we propose an equation related to Hooke's law, where the *K* value is calculated using the slopes of each graph and *Δε* is taken from 0.1 to 0.3%*ε*. Therefore the calculated slope of the real stress-strain

it to the deformation in the material.

*diameter; (e) beam diameter at t = ti + Δt.*

**4. Results**

**Figure 3.**

samples.

**35**

the focal length differences following Eq. (7).

*Digital Image Processing Applied to Optical Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

When we finish the mechanical test, a post-processing is carried out:


To do the image processing, we used some Matlab® toolbox functions, mention, and describe them:


Along **Figure 3**, we can observe how the code works: in **Figure 3a** it is seen how *imread* works; when we call it, the image is loaded; for **Figure 3b** we can observe how *rgb2grey* converts the colored image into a grayscale; in **Figure 3c** we called

*Digital Image Processing Applied to Optical Measurements DOI: http://dx.doi.org/10.5772/intechopen.88704*

**Figure 3.**

λ = 632 nm and a power of 1 mW, is placed 15 cm from the sample (Sa) in order to irradiate the polished face, and then, the reflected beam impacts a screen (Sc) which is placed 10 cm beside the laser beam and parallel to the sample. The material is first completely flat, with no stress, and as it is compressed with the universal machine (M), the size of the cross section will decrease as a reaction (**Figure 2b**). During the tests, the deformation rate was 0.2 cm/s with duration of 5 min approximately, and while they were taking place, the reflected light was recorded with a video camera (Vc) which is placed in front and parallel to the screen aside the sample. All tests were done at room temperature, and samples were compressed

*Experimental setup: (a) material with no strain; (b) the strained material will scatter light.*

When we finish the mechanical test, a post-processing is carried out:

3.A program is written in Matlab® for digital image processing.

• *rgb2grey*: this code turns the colored image into a gray scale level.

• *im2bw*: this is necessary to convert the image matrix into a binary code.

• *bwarea*: we call this function at the end of the code to calculate by binary code

Along **Figure 3**, we can observe how the code works: in **Figure 3a** it is seen how *imread* works; when we call it, the image is loaded; for **Figure 3b** we can observe how *rgb2grey* converts the colored image into a grayscale; in **Figure 3c** we called

videos with an average of 54 frames per second.

• *imread*: this function loads the image as a matrix.

the total area of the reflected laser beam.

1.First, we need to transfer the video from the video card to the computer.

2.Then, we use the free software called *Video to JPG converter*; it helps us to divide the entire video into each frame. The software enabled us to turn our

To do the image processing, we used some Matlab® toolbox functions, mention,

around 2 mm.

**Figure 2.**

*Digital Imaging*

and describe them:

**34**

*Image processing: (a) how the image is loaded; (b) grayscale image; (c) binary image; (d) initial beam diameter; (e) beam diameter at t = ti + Δt.*

*im2bw* to do a binary code of the image. At this step we can apply the function *bwarea* to know the area of this "circle." We can observe in **Figure 3d** and **e** how the beam diameter increases, for a *ti* = 0 s and *tf* = Δt s. By following the methodology described and using Eq. (6), we can obtain the change in the focal length and relate it to the deformation in the material.

The code was written to repeat this process iteratively for every image. We propose that the calculated area is the parameter *ω*(*zp*). Therefore Eq. (6) is used to calculate the focal length. We know that as a result of doing a mechanical test, we can obtain a stress-strain diagram; it is calculated through strain gages and special software incorporated to the machine. At the end of the iterative process, we plot the focal length differences following Eq. (7).

### **4. Results**

Along this chapter, we have mentioned that we prepared four steel samples and three aluminum probes, to test them under standardized compression test; for each one we obtained two graphs of interest, the stress-strain diagram and the DIP plot; **Figure 4** shows these two graphs of interest. **Figures 4**–**8** correspond to steel samples.

In **Figure 4a** we show the stress-strain diagram obtained from the software of the machine; in **Figure 4b** we can observe the related stress-focal length diagram obtained by following the experimental procedure and using Eq. (9). In the graphs we can observe a linear behavior. In section 2.2, we propose an equation related to Hooke's law, where the *K* value is calculated using the slopes of each graph and *Δε* is taken from 0.1 to 0.3%*ε*. Therefore the calculated slope of the real stress-strain

**Figure 4.** *Stress-strain diagrams: (a) obtained from universal machine, (b) experimentally from optical treatment.*

**Figure 5.** *Comparison between stress-strain diagrams for steel's sample 1.*

diagram is 0.00351 and for the experimental data is 0.001401; *K* coefficient between them is 2.5.

In **Figure 5** we show the first steel sample comparison between the real strain (Real Def.) and the experimental related-strain (Exp. Def.); both are represented in % as strain has no dimensions, and also some statistical results are shown in **Table 1**. From here, the following graphs are going to be presented the same way.

is due to the material surface, as it was hand-sanded and polished; light is not

In **Table 5** we summarize the accuracy and mean error for this four strain measurements; we can observe that the mean accuracy of our methodology is 89%. With this novel technique, taking the focal length as a parameter for measuring the stress-strain relation of a material, we proved an 89% effectiveness following standardized compression tests using steel samples; also we obtained a *K* value for this material of 2.5, which establishes the relation between the slopes of the real and

In this section, we present the statistical results and stress-strain diagrams obtained from aluminum samples. We are following the same process by comparing the DIP graph with the real stress-strain obtained from the universal machine

reflecting equally in every sample of material.

*Comparison between stress-strain diagrams for steel's sample 3.*

*Comparison between stress-strain diagrams for steel's sample 2.*

*Digital Image Processing Applied to Optical Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

the experimental graphs [25].

**Figure 6.**

**Figure 7.**

**37**

We repeat this procedure for the next three samples, obtaining *K* factor and comparing the real stress-strain diagram with the experimental plot and showing the statistical results for each one.

For the second test, we can see in **Table 2** that *K* factor remains in a value of 2.5; in this test we obtained an accurate result of 88% as compared to the real graph. In the third and fourth tests, we see in **Tables 3** and **4** that the value of the proposed *K* factor decreases to 2.4, but we obtained the same approach in reproducing the stress-strain diagram.

In **Figures 5** and **7**, the behavior of the experimental graph is above the real stress-strain diagram; for **Figures 6** and **8**, the behavior has an intersection—it first starts below then crosses above it. We propose that the behavior of the area increase *Digital Image Processing Applied to Optical Measurements DOI: http://dx.doi.org/10.5772/intechopen.88704*

**Figure 6.** *Comparison between stress-strain diagrams for steel's sample 2.*

**Figure 7.** *Comparison between stress-strain diagrams for steel's sample 3.*

is due to the material surface, as it was hand-sanded and polished; light is not reflecting equally in every sample of material.

In **Table 5** we summarize the accuracy and mean error for this four strain measurements; we can observe that the mean accuracy of our methodology is 89%.

With this novel technique, taking the focal length as a parameter for measuring the stress-strain relation of a material, we proved an 89% effectiveness following standardized compression tests using steel samples; also we obtained a *K* value for this material of 2.5, which establishes the relation between the slopes of the real and the experimental graphs [25].

In this section, we present the statistical results and stress-strain diagrams obtained from aluminum samples. We are following the same process by comparing the DIP graph with the real stress-strain obtained from the universal machine

diagram is 0.00351 and for the experimental data is 0.001401; *K* coefficient

*Stress-strain diagrams: (a) obtained from universal machine, (b) experimentally from optical treatment.*

From here, the following graphs are going to be presented the same way.

In **Figure 5** we show the first steel sample comparison between the real strain (Real Def.) and the experimental related-strain (Exp. Def.); both are represented in % as strain has no dimensions, and also some statistical results are shown in **Table 1**.

We repeat this procedure for the next three samples, obtaining *K* factor and comparing the real stress-strain diagram with the experimental plot and showing

In **Figures 5** and **7**, the behavior of the experimental graph is above the real stress-strain diagram; for **Figures 6** and **8**, the behavior has an intersection—it first starts below then crosses above it. We propose that the behavior of the area increase

For the second test, we can see in **Table 2** that *K* factor remains in a value of 2.5; in this test we obtained an accurate result of 88% as compared to the real graph. In the third and fourth tests, we see in **Tables 3** and **4** that the value of the proposed *K* factor decreases to 2.4, but we obtained the same approach in reproducing the

between them is 2.5.

**Figure 5.**

**Figure 4.**

*Digital Imaging*

stress-strain diagram.

**36**

the statistical results for each one.

*Comparison between stress-strain diagrams for steel's sample 1.*

#### **Figure 8.**

*Comparison between stress-strain diagrams for steel's sample 4.*


software. **Figure 9** and **Table 6** show the results for the first aluminum sample; with this material we obtained a *K* value of 4.52 and an approach of 90%.

**Sample Accuracy (%)** *K* 1 89 2.5 2 89 2.5 3 89 2.4 4 89 2.4

*K* factor 2.4 Young's modulus *E* 187 13 GPa % difference from 200 GPa 11% Standard deviation 0.013

*Digital Image Processing Applied to Optical Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

*Statistical results for sample 4 during the compression tests and Gaussian beam analysis.*

Using aluminum samples, we obtained similar approximations to the steel ones. In **Table 7** we can observe a *K* value of 4.54 and a slightly less approach than sample 1, obtaining 89%. In **Table 8** we observe a *K* value of 4.5 and an approach of 88%. By comparing **Figures 9**–**11** to steel results, we can see that aluminum alloys behave differently than steel ones. In these figures we can appreciate that the experimental graph goes above the real diagram at the same point (0.08%). In **Figures 10** and **11**, the experimental result goes first below, then above the real stress-strain diagram. These differences, such as steel ones, are due to the sample surface; during the process of making mirror-grade samples, we can accidentally cause nonuniform surfaces, and thus light reflection is not always the same.

In **Table 9** we summarize the accuracy and mean error for these three measure-

ments; it is seen that the mean accuracy of our method is 89%.

*Comparison between stress-strain diagrams for aluminum's sample 1.*

**Fourth sample**

**Table 4.**

**Table 5.**

**Figure 9.**

**39**

*Accuracy and error of the four tests.*

#### **Table 1.**

*Statistical result for steel's sample 1 during the compression tests and Gaussian beam analysis.*


#### **Table 2.**

*Statistical results for steel's sample 2 during the compression tests and Gaussian beam analysis.*


#### **Table 3.**

*Statistical results for steel's sample 3 during the compression tests and Gaussian beam analysis.*


**Table 4.**

*Statistical results for sample 4 during the compression tests and Gaussian beam analysis.*


**Table 5.**

**Figure 8.**

*Digital Imaging*

**Table 1.**

**Table 2.**

**Table 3.**

**38**

**First sample**

**Second sample**

**Third sample**

*Comparison between stress-strain diagrams for steel's sample 4.*

*K* factor 2.5 Young's modulus E 188 12 GPa % difference from 200 GPa 11% Standard deviation 0.011

*Statistical result for steel's sample 1 during the compression tests and Gaussian beam analysis.*

*K* factor 2.5 Young's modulus E 188 12 GPa % difference from 200 GPa 11% Standard deviation 0.010

*Statistical results for steel's sample 2 during the compression tests and Gaussian beam analysis.*

*K* factor 2.4 Young's modulus *E* 188 12 GPa % difference from 200 GPa 11% Standard deviation 0.012

*Statistical results for steel's sample 3 during the compression tests and Gaussian beam analysis.*

*Accuracy and error of the four tests.*

software. **Figure 9** and **Table 6** show the results for the first aluminum sample; with this material we obtained a *K* value of 4.52 and an approach of 90%.

Using aluminum samples, we obtained similar approximations to the steel ones. In **Table 7** we can observe a *K* value of 4.54 and a slightly less approach than sample 1, obtaining 89%. In **Table 8** we observe a *K* value of 4.5 and an approach of 88%.

By comparing **Figures 9**–**11** to steel results, we can see that aluminum alloys behave differently than steel ones. In these figures we can appreciate that the experimental graph goes above the real diagram at the same point (0.08%). In **Figures 10** and **11**, the experimental result goes first below, then above the real stress-strain diagram. These differences, such as steel ones, are due to the sample surface; during the process of making mirror-grade samples, we can accidentally cause nonuniform surfaces, and thus light reflection is not always the same.

In **Table 9** we summarize the accuracy and mean error for these three measurements; it is seen that the mean accuracy of our method is 89%.

**Figure 9.** *Comparison between stress-strain diagrams for aluminum's sample 1.*


#### **Table 6.**

*Statistical results for aluminum's sample 1 during the compression tests and Gaussian beam analysis.*


**Table 7.**

*Statistical results for aluminum's sample 2 during the compression tests and Gaussian beam analysis.*


Summarizing these tests, steel and aluminum, we observed that the present methodology can measure the stress-strain relation of a ductile material with an average error of 11%. We established that by recording the video, we obtained 54 images in average; therefore, we calculate the mean differences between measurement in a 10-s time interval, thus to obtain the sensitivity of this methodology,

**Probe Accuracy (%)** *K* 1 90 4.54 2 89 4.52 3 88 4.5

The stress-strain relation of two ductile materials under standardized compres-

sion test was measured using one laser beam, its reflection, and digital image processing treatment; we obtained 89% of accuracy and a relation of 27

ison between the reflected area and the experimental results is congruent.

The authors declare no conflict of interest for this chapter.

microdeformation. Also it was demonstrated that the diameter of the laser reflection has an increment according to the increase in the focal length; we also showed following Fresnel diffraction and the transmission through a lens that the compar-

Alonso Saldaña Heredia wants to thank CONACYT for the grant no. 360140.

which resulted to 27 microdeformations.

*Accuracy and error of the three tests.*

*Comparison between stress-strain diagrams for aluminum's sample 3.*

*Digital Image Processing Applied to Optical Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.88704*

**5. Conclusion**

**Figure 11.**

**Table 9.**

**Acknowledgements**

**Conflict of interest**

**41**

#### **Table 8.**

*Statistical results for aluminum's sample 3 during the compression tests and Gaussian beam analysis.*

With this developed technique, taking the focal length as a parameter for measuring the stress-strain relation of a material, we proved an 89% effectiveness following standardized compression tests using aluminum samples; we also obtained a *K* value for this material of 4.5, which stablishes the relation between the slopes of the real and the experimental graphs [26].

*Digital Image Processing Applied to Optical Measurements DOI: http://dx.doi.org/10.5772/intechopen.88704*

#### **Figure 11.**

*Comparison between stress-strain diagrams for aluminum's sample 3.*


#### **Table 9.**

*Accuracy and error of the three tests.*

Summarizing these tests, steel and aluminum, we observed that the present methodology can measure the stress-strain relation of a ductile material with an average error of 11%. We established that by recording the video, we obtained 54 images in average; therefore, we calculate the mean differences between measurement in a 10-s time interval, thus to obtain the sensitivity of this methodology, which resulted to 27 microdeformations.

## **5. Conclusion**

The stress-strain relation of two ductile materials under standardized compression test was measured using one laser beam, its reflection, and digital image processing treatment; we obtained 89% of accuracy and a relation of 27 microdeformation. Also it was demonstrated that the diameter of the laser reflection has an increment according to the increase in the focal length; we also showed following Fresnel diffraction and the transmission through a lens that the comparison between the reflected area and the experimental results is congruent.

### **Acknowledgements**

Alonso Saldaña Heredia wants to thank CONACYT for the grant no. 360140.

### **Conflict of interest**

The authors declare no conflict of interest for this chapter.

With this developed technique, taking the focal length as a parameter for mea-

obtained a *K* value for this material of 4.5, which stablishes the relation between the

suring the stress-strain relation of a material, we proved an 89% effectiveness following standardized compression tests using aluminum samples; we also

slopes of the real and the experimental graphs [26].

*Comparison between stress-strain diagrams for aluminum's sample 2.*

**First sample**

*Digital Imaging*

**Second sample**

**Third sample**

**Table 6.**

**Table 7.**

**Table 8.**

**Figure 10.**

**40**

*K* factor 4.52 Young's modulus *E* 63 7 GPa % difference from 70 GPa 10% Standard deviation 0.015

*K* factor 4.54 Young's modulus *E* 62 8 GPa % difference from 70 GPa 11% Standard deviation 0.013

*K* factor 4.5 Young's modulus *E* 61 9 GPa % difference from 70 GPa 12% Standard deviation 0.018

*Statistical results for aluminum's sample 1 during the compression tests and Gaussian beam analysis.*

*Statistical results for aluminum's sample 2 during the compression tests and Gaussian beam analysis.*

*Statistical results for aluminum's sample 3 during the compression tests and Gaussian beam analysis.*
