Next Article in Journal
A Novel Decision-Making Approach under Complex Pythagorean Fuzzy Environment
Previous Article in Journal
Hedging Crop Yields Against Weather Uncertainties—A Weather Derivative Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

The Archimedes’ Constant, π Seen by Mechanical Engineers

1
Mechanical Engineering Dept, South Ural State University, Prospekt Lenina, 76, Chelyabinsk 454080, Chelyabinskaya Oblast, Russia
2
Center for Computational Mechanics South Ural State University, Ulitsa Ordzhonikidze, 50, Chelyabinsk 454091, Chelyabinskaya Oblast, Russia
3
Technology Center for Strength of Structures, The Agency for the Assessment and Application of Technology, Building #220, Region Puspiptek Serpong, Banten, Indonesia
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(3), 72; https://doi.org/10.3390/mca24030072
Submission received: 26 June 2019 / Revised: 20 July 2019 / Accepted: 22 July 2019 / Published: 23 July 2019
(This article belongs to the Section Engineering)

Abstract

:
Probabilistic engineering mechanics is used to relate the value of π with one of the main parameters in fracture mechanics. It proposes an engineering method to obtain the value of it from engineering data without involving any geometrical Euclidian’s circle’s data measurement or analysis. It is the first trial in studying the use of fracture mechanics to determine the value of ratio of circumference and diameter of Euclidean’s circles indirectly, and subsequently evaluate the number of the digits actually needed in fracture mechanics and engineering purposes.

1. Introduction

Probabilistic engineering mechanics is used in this exploratory article partly to commemorate the death of a Nobel Laurette on 1 February 2012, may she rest in peace. Born 2 July 1923 in Prowent, Poland, the Nobel Laurette, Wislawa Szymborska, had always fascinated with π and, therefore, she wrote a poem on it [1]. The beginning of the poem is:
π T h e   a d m i r a b l e   n u m b e r   π :   t h r e e   p o i n t   o n e   f o u r   o n e .   A l l   t h e   f o l l o w i n g d i g i t s   a r e   a l s o   j u s t   a   s t a r t ,   f i v e   n i n e   t w o   b e c a u s e   i t   n e v e r   e n d s .   I t   c a n t b e   g r a s p e d ,   s i x   f i v e   t h r e e   f i v e ,   a t   a   g l a n c e ,   e i g h t   n i n e ,   b y   c a l c u l a t i o n , s e v e n   n i n e ,   t h r o u g h   i m a g i n a t i o n ,   o r   e v e n   t h r e e   t w o   t h r e e   e i g h t   i n   j e s t , o r   b y   c o m p a r i s o n   f o u r   s i x   t o   a n y t h i n g   t w o   s i x   f o u r   t h r e e   i n   t h e   w o r l d . ( c o n t i n u e d )
Although it is commonly known as 3.14, the value of ‘ π ’is i r r a t i o n a l and can never be expressed exactly as a fraction of integers in both the numerator and denominator. It also can never be expressed as a finite sequence of algebraic operations on integers. It is a non-constructible, as it is impossible to square the circle as was proved by lots of great mathematicians, such as Lindemann back in 1882 [2]. Throughout history, there has been much effort to determine its value more accurately. The fascination with it is carried over into cultures other than mathematics. Possibly because of the simplicity of its definition, it became more entrenched in popular culture than almost any other mathematical concepts [3] since 1650 B.C, as it is recorded in the Egyptian rhind papyrus that shows π 4 × 8 9 2 [4]. It is also a common ground among researchers in other fields as well. It is truly the world’s most mysterious number [5]. Reports on the latest calculations of it are common news items, and some report the accuracy up to over several trillion decimal digits, and as you read this article it is probably even more because the computer is still calculating it [6].
There is no doubt that a lot of formulae used in mathematics, economics, science, and engineering use π . This emphasizes the importance and significance of the mathematical constant, π [7,8]. Researchers from almost any kind of field are fascinated to know the exact value of it. This short technical note discusses its basic concept briefly, historical chronology and scientific quest of research on it, followed by our original concept of the non-geometrical approach on the approximation of the π . Using the reverse engineering concept, we formulate the number of the digits of π needed for engineering computation.

2. Basic Concept

The letter π is the first letter in a Greek word of ‘ π ε ρ ι φ ε ´ ρ ε ι α ’, periphery, or ‘ π ε ρ ι μ ε τ ρ o ς ’, perimeter, or circumference. Figure 1 shows the basic concept of the traditional computation of constant π . It is defined as:
π = C 2 r
or
π = A r 2 ,
where C is the circumference, A is the area, and r is the radius of a Euclidean’s circle. Although the definition is so simple, the value of π can never be calculated analytically. The truncation of decimal of π to 11 is known to be good enough to estimate the circumference of any circle that fits inside the earth with an error of less than one millimeter, and the truncation to 39 decimals is known to be sufficient to estimate the circumference of any circle that fits in the observable universe with precision comparable to the radius of a hydrogen atom [9]. Therefore, from the engineering point of view, finding the exact value of π seems to be out of necessity. However, engineers are human too. Their satisfactions and fascinations are not necessarily expressible in words or in equations, although only a few engineers are capable of expressing their feeling on their fascinations like the Nobel Laurette in literacy, whose poem is quoted in the introduction part.

3. Compilation of the Values of π

Table 1 shows major findings of the empirical π , although in reality many more than those are available; here, the compilation from representative sources is tabulated in the chronological order [4,10,11,12,13,14,15]. It is worth noting that those empirical values are not necessarily derived based on the geometry anymore, instead they can be derived from other disciplines. Several studies deriving the value of π without a geometrical approach are also available, e.g., economics [10]. For our purpose, the value of π used for error estimation is based on:
π = 3.141592653589 ,
which is 13 significant numbers. This value was a computed value based on the formula developed by Viete [3].

4. Obtaining π from Engineering Data

The concepts of fracture mechanics are concerned with the basic ideas for developing methods to predict the load-carrying capabilities of structures containing cracks, although in here the concept is used solely for the purpose of deriving the value of ‘empirical- π ’ from fracture mechanics point of view. The approach is based on a mathematical description of the characteristic stress field surrounding a crack in a loaded body, see Figure 2. To explore the characteristics of the stress field surrounding a crack in a loaded body, one can start from the Westergaard function [16]. The Westergaard function is a complex solution to the Airy stress functions. To do this, consider a coordinate system x , y , z in a stressed solid. At each point ( x , y , z ) , one can define the stresses σ x , σ y , σ z , τ x y , τ x z , τ y z . Neglecting the body forces, for two-dimensional problems, the equilibrium equations are
σ x x + τ x y y and σ y y + τ x y x .
In addition, recall their relations with the elastic strain:
ϵ x = u x ϵ y = v y γ x y = u y + v x .
The equilibrium equation in (4) is automatically satisfied if
σ x = 2 ψ y 2 σ y = 2 ψ x 2 τ x y = 2 ψ x y
with the stress–strain relation
E ϵ x = σ x ν σ y E ϵ y = σ x ν σ x E γ x 2 ( 1 + ν ) = τ x y .
The function ψ is called the Airy stress function. Substitution of Equations (5) and (6) into Equation (7) and differentiating twice leads to the following compatibility:
2 ( 2 ψ ) = 0 .
One can then define a complex function
ψ = Re Z ¯ ¯ + y Im Z ¯ ,
where d Z ¯ ¯ d z = Z ¯ , d Z ¯ d z = Z , and d Z d z = Z . The position at which the equation valid is in front of the crack. With the Cauchy–Riemann equation, it follows that
2 Re Z ¯ ¯ = 2 Im Z = 0 .
Therefore, the stress becomes:
σ x = Re Z y Im Z σ y = Re Z + y Im Z . τ x y = y Re Z
The function is analytic except for ( a x a , y = 0 ) . The boundary stresses follow z , the results are σ x = σ y = σ and τ x y = 0 , and on the crack surface σ = τ x y = 0 , which means that the boundary conditions are satisfied. It is more convenient to convert to a coordinate system with the origin at the crack tip, hence z should be replaced by ( z + a ) . Turning to the general problem, which is shown in Figure 3 that has the form:
Z = f ( z ) z .
The required real and constant value of f ( z ) at the crack tip is given by the notation K I , hence
Z = K I 2 π z .
With some mathematical manipulation [16], the stress singularity in front of a crack tip can then be expressed by
σ i j = K I 2 π r f i j ( θ )
or
π = K I 2 2 r σ i j 2 f i j 2 ( θ ) ,
where r and θ are the polar coordinates of a point with respect to the crack tip, and K I is the stress intensity factor, see Figure 3. Taking the polar coordinate from the origin, Equation (15) can then be simplified for our purpose to only an element in the vertical direction:
π = K I 2 2 r σ y y c o s 2 θ 2 1 + s i n 2 θ 2 s i n 2 3 θ 2 2 .
Here, Equation (16) is our tool to evaluate the value of the π . The convergence of the value can then be used to evaluate to determine how many decimals needed from the engineering mechanics point of view.

5. Computational and Experimental Approaches

The computational approach is legitimate to be used here because K and σ y y are independently calculated without the need to use geometrical information pertaining to Eucledian’s circle. To make engineering sounds of the analysis, a computational model that is close to a real engineering experiment is constructed in addition to an experiment. The model here is homogeneous material having the shape of a CT (compact tension) specimen according to ASTM E647 and E08 [17,18] standards. The following steps were taken:
  • Model creation, load and boundary condition applications. For this purpose, two different meshing sizes were employed, namely rough (approximately 2000 nodes) and fine (approximately 6000 nodes).
  • Obtain the nodal/element information of θ and r to compute each individual value of K I 2 2 r c o s 2 θ 2 1 + s i n 2 θ 2 s i n 2 3 θ 2 2 .
  • Obtain the nodal/element values of σ y y .
  • Compute the value of π according to Equation (16) for each individual node/element.
  • Generate the graph to evaluate the convergence value of π using Equation (16).
  • Go back to step 1, mesh the model with different mesh size, and follow the same procedure.
The results are basically two different Gaussian like curves (Figure 5). Based on these curves, the FWHM ( full width at half maximum) can then be estimated. In our case here, the rough meshing, the FWHM is approximated by the following the lines shown in the figure:
  • Generalized estimation:
    f ( x ) = 1 σ 2 π · Exp ( x x 0 ) 2 2 σ 2 ,
    f ( x ) = C 1 · Exp ( x x 0 ) 2 C 2 · σ 2 ,
  • Curve fitting for data obtained from rough meshing (based on 2000 data points):
    f ( x ) = 52 · Exp ( x π ) 2 3.1 × 0.03 2 ,
    FWHM 0.086 . Therefore , our ( engineering π ) ( true π ) ± 0.043
    : 3.098 π 3.184 ,
  • Curve fitting for data obtained from fine meshing (based on 6000 data points):
    f ( x ) = 83 · Exp ( x π ) 2 1.5 × 0.03 2 ,
    FWHM 0.050 . Therefore , our ( engineering π ) ( true π ) ± 0.025
    : 3.116 π 3.166 .
Even though the FWHM here is based on the interpolation of the computational results, the experimental data within the framework of fracture mechanics also showed similarly due to the modeling that mimics the experiment. Therefore, the data obtained and used in the computational approach here simulated an ideal engineering experiment. Experimental data were not taken solely for this purpose; instead, it was for another purpose [16,19]. However, the data were then processed to relate π and K I 2 as it is described in Equation (16). The samples used were compact tension specimens. While the detailed data collection method is beyond the topic of this paper, and the readers are to consult to the reference directly [19], for the sake of the reader’s convenience, the method is briefly highlighted here. Figure 6 is the representative specimen used in the experiment. The stress intensity factor K for standard CT specimen geometry was calculated as follows: [17,18]:
K = P B W 2 + α ( 1 α ) 3 / 2 0.886 + 4.64 α 13.32 α 2 + 14.72 α 3 5.6 α 4 ,
where B and W are the specimen’s thickness and width, respectively. α is the relative crack length ( a / W ), and P is the applied load. The value of stress σ y y was obtained by the surface strain measurements outside the plastic zone area. The plastic zone sizes and definition are the ones defined and calculated in the references published previously [16], with the assumption of the plane stress governing rule. Figure 7 shows the distribution of the values of the π calculated from experimental results. In the figure, the curve fittings of the same distribution by computational approaches were also plotted.

6. Discussion

The research shows that π can be obtained empirically using a fracture mechanics concept too. This fact is not surprising as it can also be obtained by other data, including economics, and even societal human behaviors [7,8]. Our method is in a way similar to that of Buffon’s Needle [20]. It is one of the classic problems in the field of geometrical probability for finding π by the experiment. It involves dropping a needle on a lined sheet of paper and determining the probability of the needle crossing one of the lines on the page. The problem also can be solved using integral geometry. The remarkable result is that the probability is directly related to the value of π and can be used to design a Monte Carlo method for approximating the number of π .
Our observations of the answers clarify that, based on the rough meshing, 3.098 π 3.184 , while, based on the fine meshing, 3.116 π 3.166 . The difference of the result of rough meshing and fine meshing is only on the standard deviation. The fine meshing gives the smaller standard deviation, while the mean of both distributions remains the same. Thus, as the sample size increases, the standard deviation of the means decreases. This implies that actually engineering data only needs no more than three significant numbers of π . Experimental data show that an even lower digit number is actually needed. The fine meshing data shows that it requires three significant numbers, while the rough one requires only two significant numbers. This once again reiterated the fact that the truncation of decimal of π to 11 is good enough to estimate the circumference of any circle that fits inside the earth with an error of less than one millimeter, while modern engineering relies on extensive empirical equations based on experimental data that inherit large errors. The fact that many researchers and enthusiasts produce billions or even trillions of significant numbers is clearly beyond the engineering needs.

7. Conclusions

An engineering approach that does not directly involve any circle/geometrical analysis related to a circle was used to derive π here. Using the old school engineering rule of thumb, it has been agreed that conventional engineering fracture mechanics require no more than three digits of π , although nowadays computer usage automatically generates the number of digits depending on the type of the computer used as well as on the software used.

Author Contributions

This article was originally prepared by the Y.P. The A.S. also contributed substantially in the manuscript preparation.

Acknowledgments

Part of the work was performed during the first author’s stay at SUSU as visiting professor. The funds for covering the costs to publish in open access was solely the individual contributions from both authors out of their love in science and engineering.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nobelprize.org. Wislawa Szymborska—A Biography. Available online: http://www.nobelprize.org/nobelprizes/literature/laureates/1996/szymborska-bio.html (accessed on 22 July 2019).
  2. Lindenmann, F. Uber die zahl π. Math. Ann. 1982, 20, 213–215. [Google Scholar]
  3. Berggren, L.; Borwein, J.; Borwein, P. Pi: A Source Book; Springer: New York, NY, USA, 2004. [Google Scholar]
  4. Gillings, R. Mathematics in the Time of the Pharaohs; Dover: Mineola, NY, USA, 1982. [Google Scholar]
  5. Posamentier, A.S.; Lehmann, I. A Biography of the World’s Most Mysterious Number; Prometheus Books: Amherst, NY, USA, 2004. [Google Scholar]
  6. Takahashi, D. Parallel implementation of multiple-precision arithmetic and 2,576,980,370,000 decimal digits of calculation. Parallel Comput. 2010, 36, 439–448. [Google Scholar] [Green Version]
  7. Stanley, H. Introduction to Phase Transitions and Critical Phenomena; Oxford University Press: Oxford, UK, 1987. [Google Scholar]
  8. Cassimatis, P. A Concise Introduction to Engineering Economics; Taylor & Francis: Oxfordshire, UK, 1988. [Google Scholar]
  9. Young, R.M. Excursions in Calculus; Mathematical Association of America: Washington, DC, USA, 1992. [Google Scholar]
  10. Chong, T.T.L. The empirical quest for π. Comput. Math. Appl. 2008, 56, 2772–2778. [Google Scholar]
  11. Ponnusamy, S. Foundations of Mathematical Analysis; Birkhauser: New York, NY, USA, 2012. [Google Scholar]
  12. Miller, F.; Vandome, A.; John, M. π: Proof That π Is Irrational, Numerical Approximations of π, Piphilology, List of Formulae Involving P, Leibniz Formula for π, Chronology of Computation of P, Wallis Product; VDM Verlag Dr. Müller: Riga, Latvia, 2010. [Google Scholar]
  13. Newton, I.; Whiteside, D. The Mathematical Papers of Isaac Newton: Year 1684–1691; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  14. Rao, K. Srinivasa Ramanujan: A Mathematical Genius; East West Books: Madras, India, 1998. [Google Scholar]
  15. Eymard, P.; Lafon, J. The Number Pi; American Mathematical Society: Providence, RI, USA, 2004. [Google Scholar]
  16. Idris, R.; Prawoto, Y. Designing steel microstructure based on fracture mechanics approach. Mater. Sci. Eng. A 2009, 507, 74–86. [Google Scholar]
  17. ASTM. Standard test method for measurement of fatigue crack growth rates. In ASTM:E647; ASTM International Publisher: West Conshohocken, PA, USA, 2008. [Google Scholar]
  18. ASTM. Standard test methods for tension testing of metallic materials. In ASTM:E8; ASTM International Publisher: West Conshohocken, PA, USA, 2008. [Google Scholar]
  19. Prawoto, Y. Influence of ferrite fraction within martensite matrix on fatigue crack propagation: An experimental verification with dual phase steel. Mater. Sci. Eng. A 2012, 552, 547–554. [Google Scholar]
  20. Ramaley, J.F. Buffon’s Noodle Problem. Am. Math. Mon. 1969, 76, 916–918. [Google Scholar]
Figure 1. A Euclidean’s circle with a radius r, a circumference C, and an area A.
Figure 1. A Euclidean’s circle with a radius r, a circumference C, and an area A.
Mca 24 00072 g001
Figure 2. Mode I stress under bi-axial loading [16].
Figure 2. Mode I stress under bi-axial loading [16].
Mca 24 00072 g002
Figure 3. Basic usage of linear elastic fracture mechanics method [16].
Figure 3. Basic usage of linear elastic fracture mechanics method [16].
Mca 24 00072 g003
Figure 4. Model creation to compute the values of K of the entire system and σ y y , at corresponding values of θ .
Figure 4. Model creation to compute the values of K of the entire system and σ y y , at corresponding values of θ .
Mca 24 00072 g004
Figure 5. Convergence analysis of the value of π based on the computation of K I 2 2 r σ y y c o s 2 θ 2 1 + s i n 2 θ 2 s i n 2 3 θ 2 2 for a roughly meshed model (left) and a finely meshed model (right).
Figure 5. Convergence analysis of the value of π based on the computation of K I 2 2 r σ y y c o s 2 θ 2 1 + s i n 2 θ 2 s i n 2 3 θ 2 2 for a roughly meshed model (left) and a finely meshed model (right).
Mca 24 00072 g005
Figure 6. Test specimen for fatigue crack propagation [19].
Figure 6. Test specimen for fatigue crack propagation [19].
Mca 24 00072 g006
Figure 7. Experimental results using fracture mechanics data. The results are also compared with computational results.
Figure 7. Experimental results using fracture mechanics data. The results are also compared with computational results.
Mca 24 00072 g007
Table 1. Main equations pertaining to the value estimation of π .
Table 1. Main equations pertaining to the value estimation of π .
ValueInventor and RemarksSource
π 4 × ( 8 9 ) 2 Value used by ancient Egyptians (1650 B.C.)[4]
223 71 < π < 22 7 Archimedes (287–212 B.C.)[10]
3.141024 < π < 3.1415927 Zu (429–500 A.D.)[10]
π 2 × 2 2 × 2 2 + 2 . Viette (1593)[3]
π 4 k = 0 ( ( 1 ) k 2 k + 1 ) Madhava–Leibniz series. Slow to converge[11]
π 12 k = 0 ( ( 3 ) k 2 k + 1 ) Modified Madhava–Leibniz series. Better convergence, capable of producing accurate 11 digits[16]
π 2 k = 1 ( 2 k ) 2 ( 2 k ) 2 1 Walli’s product (1650)[12]
π 3 n = 0 2 n n 16 n ( 2 n + 1 ) Isaac Newton[13]
1 π 2 2 9801 k = 0 ( ( 4 k ) ! ( 1103 + 26390 k ( k ! ) 4 396 4 k ) Srinivasa Ramanujan[14]
π k = 0 1 16 k ( 4 8 k + 1 2 8 k + 4 1 8 k + 5 1 8 k + 6 ) Bailey et al.[15]

Share and Cite

MDPI and ACS Style

Prawoto, Y.; Suhartono, A. The Archimedes’ Constant, π Seen by Mechanical Engineers. Math. Comput. Appl. 2019, 24, 72. https://doi.org/10.3390/mca24030072

AMA Style

Prawoto Y, Suhartono A. The Archimedes’ Constant, π Seen by Mechanical Engineers. Mathematical and Computational Applications. 2019; 24(3):72. https://doi.org/10.3390/mca24030072

Chicago/Turabian Style

Prawoto, Yunan, and Agus Suhartono. 2019. "The Archimedes’ Constant, π Seen by Mechanical Engineers" Mathematical and Computational Applications 24, no. 3: 72. https://doi.org/10.3390/mca24030072

APA Style

Prawoto, Y., & Suhartono, A. (2019). The Archimedes’ Constant, π Seen by Mechanical Engineers. Mathematical and Computational Applications, 24(3), 72. https://doi.org/10.3390/mca24030072

Article Metrics

Back to TopTop