Acessibilidade / Reportar erro

Modeling of Cracked Beams by the Experimental Design Method

Abstract

The understanding of phenomena, no matter their nature is based on the experimental results found. In the most cases, this requires an important number of tests in order to put a reliable and useful observation served into solving the technical problems subsequently. This paper is based on independent and variables combination resulting from experimentation in a mathematical formulation. Indeed, mathematical modeling gives us the advantage to optimize and predict the right choices without passing each case by the experiment. In this work we plan to apply the experimental design method on the experimental results found by (Deokar, A, 2011Deokar, A. A.V, Wakchaure, B. V.D. (2011). Experimental Investigation of Crack Detection in Cantilever Beam Using Natural Frequency as Basic Criterion, INSTITUTE OF TECHNOLOGY, NIRMA UNIVERSITY, AHMEDABAD - 382 481), concerning the effect of the size and position of a crack on the measured frequency of a beam console, and validating the mathematical model to predict other frequencies

Keywords:
Parameters; experimental design method; modeling; frequency; crack

INTRODUCTION

Several scientific works has been conducted on the experimental design method each in its field of application and its objective, (Bounazef M, Goupy J and Castro. L, 2009Bounazef, M., Djeffal, A., Serier, M.and Adda, BEA. (2009). Optimization by behavior modeling of a protective porous material. Journal of the Computational Materials Science. 44(3):921-928.,1925Fisher,R. (1925). Statistical Methods for Research Workers. Oliver and Boyd;. PMid: 17246289,2004Castro.L, Viéville.P, and Lipinski.P, Phys. J(2006).IV France 134, 1273.). Other authors (Aît Yala, 2009A. Ait Yala , A. Megueni, (2009) Optimisation of composite patches repairs with the design of experiments method, Materials and Design pp 200-205) have implemented experimental design using numerical results. Experimental design were used for the first time by (Fisher R, 1925Fisher,R. (1925). Statistical Methods for Research Workers. Oliver and Boyd;. PMid: 17246289), in the field of agriculture where the experimental parameters are numerous and significant which leads to mathematical modeling and therefore to optimize the model sought. In mechanical terms, for example, (Serier et al, 2013aRefassi,K. Serier, M. and Lousdad ,A. (2013). DOE Analysis of Feed and Cutting Speed Effects on Frequency During Machining Operations, Advanced Science Letters Vol. 19, 942-945.) used experimental design method for optimizing the machining parameters in order to increase the life of the tool. These same authors (2013bSerier, M. Lousdad , A. Refassi, K.and Megueni,A. (2013). Analysis of Parameters Effects on Crac Breathing and Propagation in Shaft of Rotor Dynamic Systems, Materials Research. 2013; 16(4): 867-873 (Brazil).) performed a work in the field of vibration where they showed the importance of the geometrical properties, the speed of rotation of a shaft (tree) on the opening and closing mechanism of crack in a rotating machine. It is clear that the phenomenon of cracking beams has negative consequences on the large industrial constructions, (Thomas, M, 1995Thomas, M., Nguyen, H., Hamidi,L. , Massoud ,M. and Piaud J.B. (1995). Detection of rotor cracks by experimental modal analysis, Transactions of the Canadian Society of Mechanical Engineering, 19(2): pp 155-174.). The work of (Deokar, A, 2011Deokar, A. A.V, Wakchaure, B. V.D. (2011). Experimental Investigation of Crack Detection in Cantilever Beam Using Natural Frequency as Basic Criterion, INSTITUTE OF TECHNOLOGY, NIRMA UNIVERSITY, AHMEDABAD - 382 481) has focused on the experimental investigation for the detection of a crack on a cantilever used a criteria based on natural frequencies. Our work consisted on the use of experimental design method for modeling and prediction of frequencies.

2 EXPERIMENTAL RESULTS FOR EACH FREQUENCY MODE

In the tables below, the experimental results of the variation of frequency ratio as a function of position and depth of modes 1, 2 and 3 are presented. These results came from the work of (Deokar, 2011Deokar, A. A.V, Wakchaure, B. V.D. (2011). Experimental Investigation of Crack Detection in Cantilever Beam Using Natural Frequency as Basic Criterion, INSTITUTE OF TECHNOLOGY, NIRMA UNIVERSITY, AHMEDABAD - 382 481).

Table 1
Frequency ratio of the various tests of the first mode and experience matrix.

Table 2
Frequency ratio of the various tests of the second mode and experience matrix.

Table 3
Frequency ratio of the various tests of the third mode and experience matrix.

3 CALCULATIONS OF THE EFFECTS OF FACTORS

Each factor xi is affected (acted on) the behavior of the beam and it's defined by the effect of ai . It is possible that the factors interact with each other, this is the case in our work, and therefore we are left with three factors, instead of two, including the average of these latter. In other words any response yi depends on the action of all factors together xi. Analytically, and the relationship between the response factor can exist only when a certain proportionality exists between them. This leads us to write:

The Solve of system of equations is based on the least squares method, and the solution is noted a.

This solution is given by the following formula derived from the theory of matrix calculation.

Therefore

Table 4
Coefficients that make up the mathematical model

The models are thus written in the following forms

  1. Mode 1

    :y 1 = 0.92 - 0.090(X 1) + 0.076(X 2)+0.081(X 1 X 2)

  2. Mode 2

    :y 2 = 0.93 - 0.70(X 1) - 0.026(X 2) - 0.028(X 1 X 2)

  3. Mode 3

    :y 3 =0.93 - 0.70(X 1) - 0.036(X 2) - 0.037(X 1 X 2)

4 ANALYSES WITH ONE VARIABLE FACTOR

4.1 Effect of Each Factor for the Three Modes

The figures below represent the effect of the two most principals factors (the crack depth and position of the crack) for each mode of vibration. We note that the ratio of the frequency is proportional to the crack depth; the slope is negative regardless of the active mode. This confirms that the propagation of a crack causes a decrease in the ratio of the frequency. In the case of the effect of crack position on the frequency, we observed a change of the upward frequency in the first mode. In other words, the more crack is remote from the recess the more important is the frequency, however in the case of other modes (2 and 3) an opposite change is noticed with slopes more or less important. This is due to the values of the amplitudes of modes 2 and 3.

Figure 1
Effect of the depth of the crack on the frequency for the three modes.

Figure 2
Effect of the position of the crack on the frequency for the three modes.

5 INTERACTION ANALYSE

The analyze of the interaction between the depth and crack location was done by the help of a representation of iso-courbes

Table 5
The results of combined effects (mode1).

Thefigure 3arepresents the effects of the two factors versus the ratio of frequency in mode 1. We observed that there is just one near the abscissa of embedding and at low depths frequency ratio reached the value of 0.98. This ratio decreases to great depths. For abscissa near the edge of the beam, the depth of the crack has not effect. In the second and third mode analysis of iso-curves (fig.3b,3c) shows that the ratio reaches its maximum when the two factors are the lowest values (tableaux 6and7)

Figure 3
Effect of the interaction of two factors on the frequencies for the three modes.

Table 6
The results of combined effects (mode 2).

Table 7
The results of combined effects (mode 3).

6 ANALYZE WITH THREE FACTORS

The mathematic models already established, they allowed us to calculate the predicted frequencies and the residues for each mode (table 8,9and10).

Table 8
Residues for the mode 1.

Table 9
Residues for the mode 2.

Table 10
Residues for the mode 3.

With the help of the statistics calculations, we were able to define the effects the most significant of the factors and their gaps of confidence (trust), all with calculating the residues ei. The residues are the difference between the experimental value and predicted value by the mathematic model (figure 4) and they are linked by the linear regression.

Figure 4
Distribution of the experimental points from the mathematical model for each mode.

7 REALIZATION OF THE TEST OF EFFECTS SIGNIFICANCE

The test used is the Student test << t>>. An effect is said significant (that is to say that the interaction or the variable associated with it has an influence on the response), if, for a given risk significantly different from 0. So we test the hypothesis

Hypothesis:

Hi = << ai ≠ 0>>

For this, we calculate

Student table is then used to v= n - p degrees of freedom (n is the number of experiments and p is the number of effects including the constant). The risk of a first species is selected (usually 1% or 5%) and is read in the table of the Student t value, using the part of the table related to a bilateral test.

Mode 1:

The variance calculated for each effect is then:

Therefore: ti x si = 0.037

  • The effect of crack depth : |-0.090| > 0.037 → significatif

  • The effect of crack location : |+0.076| > 0.037 → significatif

  • The effect of general average : |+0.920| > 0.037 → significatif

  • The effect of the interaction : |-0.810| < 0.037 → significatif

Confidence interval of model effects of the first mode

Table 11
Confidence interval of model effects of the first mode.

Figure 5
Confidence interval of model effects of the first mode.

Mode 2:

s2 =Σ

s2 =Σ(-0.000008)2 + (-0.000004)2 + (0.00000009)2 + (0.0003)2 + (-0.00006)2 + (-0.0005)2 + (-0.0013)2 + (-0.0005)2

s2 =(0.0028)

s2 = (0.00057)

The variance calculated for each effect is then:

s2 =(0.00057)

s2 = 0.000060

si = 0.0077

Therefore: ti x si = 0.030

Then:

  • The effect of crack depth: |+0.070| > 0.032 → significatif

  • The effect of crack location: |-0.026| < 0.032 → non significatif

  • The effect of general average: |+0.930| > 0.032 → significatif

  • The effect of the interaction: |-0.028| < 0.032 → non significatif

Confidence interval of model effects of the first mode

[ai- t (α, ν) si ; ai + t (α, ν) si]

Table 12
Confidence interval of model effects of the second mode.

Figure 6
confidence interval of model effects of the second mode.

Mode 3:

s2 =Σ

s2 =Σ(-0.0004)2 + (-0.00002)2 + (0.000068)2 + (0.0003)2 + (-0.00026)2 + (-0.00007)2 + (-0.00082)2 + (-0.00059)2 + (-0.0017)2 + (-0.00015)2

s2 =(0.0038)

s2 = (0.00076)

The variance calculated for each effect is then:

s2 =(0.00076)

s2 = 0.000084

si = 0.0092

Therefore: ti x si = 0.0370

Then:

  • The effect of crack depth: |-0.070| > 0.037 → significatif

  • The effect of crack location: |-0.036| < 0.037 → non significatif

  • The effect of general average: |+0.930| > 0,037 → significatif

  • The effect of the interaction: |-0.037| > 0.032 → significatif

Confidence interval of model effects of the first mode

[ai- t (α,ν) si ; ai + t (α,ν) si]

Table 13
Confidence interval of model effects of the third mode.

Figure 7
Confidence interval of model effects of the third mode.

8 CONCLUSIONS

Using the method of experimental design allows to extract the maximum information at a reasonable cost and minimum time. In our case, we have:

  1. 1

    ) Modeled the vibration behavior of a cracked beam console

  2. 2

    ) Plotted the iso-curves for different modes of vibration

These curves allow the choice according to the needs of the operating state (condition) of the beam.

It further notes that:

The depth of the crack remains significantly regardless of the current mode vibration, while this is not the case for the location of the crack which is significant only in the first mode.

The interaction between the two factors of an unexplained way varies from one mode to another; in the first mode of interaction is significant, it is only slightly in the third mode when it no longer is in the second mode

References

  • Deokar, A. A.V, Wakchaure, B. V.D. (2011). Experimental Investigation of Crack Detection in Cantilever Beam Using Natural Frequency as Basic Criterion, INSTITUTE OF TECHNOLOGY, NIRMA UNIVERSITY, AHMEDABAD - 382 481
  • Bounazef, M., Djeffal, A., Serier, M.and Adda, BEA. (2009). Optimization by behavior modeling of a protective porous material. Journal of the Computational Materials Science. 44(3):921-928.
  • Fisher,R. (1925). Statistical Methods for Research Workers. Oliver and Boyd;. PMid: 17246289
  • Castro.L, Viéville.P, and Lipinski.P, Phys. J(2006).IV France 134, 1273.
  • A. Ait Yala , A. Megueni, (2009) Optimisation of composite patches repairs with the design of experiments method, Materials and Design pp 200-205
  • Thomas, M., Nguyen, H., Hamidi,L. , Massoud ,M. and Piaud J.B. (1995). Detection of rotor cracks by experimental modal analysis, Transactions of the Canadian Society of Mechanical Engineering, 19(2): pp 155-174.
  • Refassi,K. Serier, M. and Lousdad ,A. (2013). DOE Analysis of Feed and Cutting Speed Effects on Frequency During Machining Operations, Advanced Science Letters Vol. 19, 942-945.
  • Serier, M. Lousdad , A. Refassi, K.and Megueni,A. (2013). Analysis of Parameters Effects on Crac Breathing and Propagation in Shaft of Rotor Dynamic Systems, Materials Research. 2013; 16(4): 867-873 (Brazil).

Publication Dates

  • Publication in this collection
    Sept 2015

History

  • Received
    11 June 2014
  • Accepted
    17 Feb 2015
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br