Next Article in Journal
Effect of Heat Treatment on the Mechanical and Corrosion Properties of Mg–Zn–Ga Biodegradable Mg Alloys
Previous Article in Journal
A Recurrent Neural Network-Based Method for Dynamic Load Identification of Beam Structures
Previous Article in Special Issue
Fabrication of Functionally Graded Diamond/Al Composites by Liquid–Solid Separation Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation of Stress Distributions in 3D Printed Graded Plates with a Circular Hole

Jiangsu Province Key Laboratory of Advanced Manufacturing Technology, Huaiyin Institute of Technology, Huai’an 223003, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(24), 7845; https://doi.org/10.3390/ma14247845
Submission received: 21 August 2021 / Revised: 12 December 2021 / Accepted: 15 December 2021 / Published: 18 December 2021

Abstract

:
An experimental investigation is presented for the stress distributions in functionally graded plates containing a circular hole. On the basis of the authors’ previously constructed theoretical model, two kinds of graded plates made of discrete rings with increasing or decreasing Young’s modulus were designed and fabricated in virtue of multi-material 3D printing. The printed graded plates had accurate size, smooth surface, and good interface. The strains of two graded plates under uniaxial tension were measured experimentally using strain gages. The stresses were calculated within the range of linear elastic from the measured strains and compared with analytical theory. It is found that the experimental results are consistent with the theoretical results, and both of them indicate that the stress concentration around the hole reduces obviously in graded plates with radially increasing Young’s modulus, in comparison with that of perforated homogenous plates. The successful experiment in the paper provides a good basis and support for the establishment of theoretical models and promotes the in-depth development of the research field of stress concentration in functionally graded plates.

1. Introduction

Functionally graded materials (FGMs) are a class of advanced composite materials in which the material properties change uninterruptedly along one or more directions. FGMs are characterized by a compositional gradient of one material into another, which is totally different from the conventional composite materials. The continuous change in microstructure of FGMs avoids the mismatch of material properties across the interface, and thus stress concentration can be effectively reduced compared with that existing at material interfaces [1,2]. Especially when the structures contain various holes, the stress concentrations around the holes can be decreased by choosing the properly radial changes of the elastic properties [3,4].
A lot of research has been carried out on the problem of stress concentration around holes in FGM plates with radially varying elastic properties in the past 20 years. Two analytical approaches to this problem are mainly employed. One is the numerical investigations on the basis of the finite element method. The other is rigorous theoretical analyses according to the linear elasticity theory.
In the aspect of numerical investigations, Venkataraman et al. [5] proposed a two-dimensional numerical model of the plate with a hole having a radially varying elastic modulus, inspired by the special material distributions near the blood vessel holes in bone. They explained how bones design holes by optimizing the structure of the perforated plates with biological variables. To reduce stress concentration and increase load-carrying capacity of the plate with holes, Huang et al. [6,7] also optimized the material distribution near the hole by mimicking bones through axisymmetric and nonaxisymmetric FGMs. On the basis of the finite element method, Kubair and Bhanu-Chandar [8] and Wang et al. [9] studied the stress concentration coefficient (SCF) in FGM plates with circular and elliptical holes, respectively, under uniform tensile traction. The elastic properties varying along three directions (X, Y, and radial direction) are discussed. Goyat et al. [10,11] analyzed the reduction of SCF around a rectangular hole and a pair of circular holes in homogenous plates by using a FGM layer based on the extended finite element method. Berezhnoi et al. [12,13] presented the exact solutions of stress distribution in the graded layers of super flywheel and discussed the effect of the relative radius of the holes in the flywheel on the specific energy.
In the aspect of theoretical analyses, Zhang et al. [14] proposed an analytical model of radial FGM plates with a circular hole and derived exact thermal stress solutions for the plates under axisymmetric thermal loading. Yang et al. [15,16] developed the theoretical models for the cases of non-axisymmetric loads, including arbitrary uniform tension and elastic wave. They decomposed the continuous FGM plates into a homogeneous plate containing multi-rings with different elastic constants and then solved the elastic fields on the basis of complex function theory. Mohammadi et al. [17] presented the general solution of SCF around a circular hole in a radially inhomogeneous plate under uniaxial tension, biaxial tension, and pure shear. Sburlati et al. [18] proposed a method to relieve the SCF around holes in homogeneous plates by using radial FGM ring and presented an analytical solution of SCF, as Young’s modulus varies with a monotonic power law. Kubair [19,20] considered the SCF and stress-gradients due to a circular hole in a radial FGM plate under anti-plane shear and presented the exact expressions for the elastic fields. Nie et al. [21] studied the problems of material tailoring for reducing SCF at circular holes in a FGM plate. Nie and Batra [22] further studied the reduction of SCF in homogeneous plates with holes by inserting a functionally graded incompressible material layer.
However, few experiments have been performed on the topic compared with present theoretical and numerical investigations. It is well known that the experiments play a fundamental and supporting role in scientific research. The first-hand experimental data can provide basis and support for the establishment of theoretical models and promote the in-depth development of the corresponding research field. The activities in the experimental research on stress concentration in FGM plates are limited mainly by the complexity and difficulty of the manufacturing technology of specimens. Buskirk et al. [23] experimentally studied the strength of a FGM biomimetic plates with a hole, which is designed based on the functional gradation in elastic modulus properties, as observed in the bone. The FGM plate was fabricated by different foam types, such as cellular, solid, and high density. Combined with finite element analyses, they verified that the FGM biomimetic plate has a higher load carrying capacity than a homogenous plate. In this paper, we present the experimental investigation for the stress distributions in graded plates fabricated by multi-material 3D printing. The surface morphology of the printed plates was observed by optical microscope. The objective of this paper is to present the design, fabrication, analysis, and results from mechanical tests of FGM plates with a circular hole and compare the results to the theoretical predictions given by the authors’ previously theoretical model.

2. Theoretical Model

In the authors’ previous work [15], the theoretical model was developed for the stress concentration in an infinite FGM plate with a circular hole subjected to remote uniform loads σ x , σ y , and τ x y , as shown in Figure 1. σ x , σ y , and τ x y are the normal stress and shear stress at infinity, respectively. The radius of the hole is denoted by r 0 . The elastic properties are assumed to change continuously and arbitrarily in the radial direction. With the method of piece-wise homogeneous layers, the domains of FGM plate can be decomposed into N homogeneous rings Ω ( 1 ) , Ω ( 2 ) , Ω ( j ) , Ω ( N ) with equal width and an outer homogeneous plate. When the number of rings N is large enough, the radially continuous Young’s modulus E ( r ) and Poisson’s ratio ν ( r ) is approximately regarded as constant E j and ν j ( j = 1 , 2 , N ) in each ring. Therefore, the boundary condition of force and displacement at the surface of the hole and the interface between the rings Ω ( j ) and Ω ( j + 1 ) can be expressed as
X n = Y n = 0 ,
X n ( j ) = X n ( j + 1 ) ,   Y n ( j ) = Y n ( j + 1 ) ,
u ( j ) = u ( j + 1 ) ,   v ( j ) = v ( j + 1 ) ,
where X n , Y n , u , and v symbolize the components of force and displacement.
Based on the theory of the complex variable functions, the boundary conditions for a plane problem of a homogeneous solid in a fixed rectangular coordinate system ( x , y ) can be expressed as [24]
φ ( z ) + z φ ( z ) ¯ + ψ ( z ) ¯ = ± i 0 s ( X n + i Y n ) d s ,   z = x + i y ,
κ φ ( z ) z φ ( z ) ¯ ψ ( z ) ¯ = 2 G ( u + i v ) ,
where φ ( z ) and ψ ( z ) stand for the potential functions. κ and G are elastic constants with different values for the cases of plane stress and plane strain.
According to Equations (4) and (5), Equations (1)–(3) can be expressed as
φ 1 ( z ) + z φ 1 ( z ) ¯ + ψ 1 ( z ) ¯ = 0 ,
φ j ( z ) + z φ j ( z ) ¯ + ψ j ( z ) ¯ = φ j + 1 ( z ) + z φ j + 1 ( z ) ¯ + ψ j + 1 ( z ) ¯ ,
1 G j [ κ j φ j ( z ) z φ j ( z ) ¯ ψ j ( z ) ¯ ] = 1 G j + 1 [ κ j + 1 φ j + 1 ( z ) z φ j + 1 ( z ) ¯ ψ j + 1 ( z ) ¯ ] ,
The complex potential functions in each circular ring Ω ( j ) ( j = 1 , 2 , N ) and the plate can be expressed as [24]
φ j ( z ) = a k ( j ) ( R 0 z ) k ,   ψ j ( z ) = b k ( j ) ( R 0 z ) k ,
φ N + 1 ( z ) = τ z + k = 0 a k ( N + 1 ) ( R 0 z ) k ,   ψ N + 1 ( z ) = τ z + k = 0 b k ( N + 1 ) ( R 0 z ) k ,
where a k ( j ) , b k ( j )   a k ( N + 1 ) , b k ( N + 1 ) are the unknown coefficient and R 0 is a reference radius. τ and τ are the constants dependent on the remote stresses as τ   =   ( σ x + σ y ) / 4 , τ   =   ( σ y σ x ) / 2 + i τ x y .
Substituting Equations (9) and (10) into Equations (6)–(8), we can get a set of linear equations containing all unknown coefficients a k ( j ) , b k ( j ) , a k ( j ) , b k ( j ) and a k ( N + 1 ) , b k ( N + 1 ) . These equations can be applied to determine these unknown coefficients. In this case, the general solution of the stress fields in each ring and the homogenous plate can be derived based on the following field equations for the two-dimensional problem of the solid [24]
σ y + σ x = 4 Re [ φ ( z ) ] ,
σ y σ x + 2 i τ x y = 2 [ z ¯ φ ( z ) + ψ ( z ) ] ,
where σ x , σ y , and τ x y are the components of stresses.
In the work [15], authors discussed the effect of different functions of the Young’s modulus on the stress distribution in the plate. On the base of analysis and comparison, the desired optimal distribution of Young’s modulus in the selected functions is the following exponential form:
E ( r ) = E 0 ( 1 0.8 e e r / r 0 ) ,
where E 0 is a constant. As the Young’s modulus changes in the above function, there is no stress concentration around the hole, and the stress distribution is almost uniform in the plate.

3. Experimental Procedures

In this paper, the experiment was designed and conducted based on the theoretical model. Two kinds of FGM plates made of discrete rings were fabricated for the cases of Young’s modulus, increasing and decreasing along the radial direction in virtue of multi-material 3D printing. The PolyJet 3D Printer named Stratasys Objet260 Connex3 (Stratasys Ltd., Rehovot, Israel) was used for printing the specimens of the FGM plate. This 3D Printer had a system size of 870 × 735 × 1200 mm and maximum build size of (XYZ) 255 × 252 × 200 mm. The printing accuracy reached up to 0.2 mm for full model size, and minimum build layer thickness could be as fine as 16 μm.

3.1. Tension Test of the Printing Materials

Two types of base materials (VeroWhitePlus and TangoBlackPlus, Stratasys Ltd., Rehovot, Israel) were used by the 3D Printer. VeroWhitePlus was rigid at room temperature and made from isobornyl acrylate, acrylic monomer, urethane acrylate, epoxy acrylate, acrylic monomer, acrylic oligomer, and a photoinitiator. TangoBlackPlus was rubbery at room temperature and formed from urethaneacrylate oligomer, Exo-1,7,7-trimethylbicyclo hept-2-yl acrylate, methacrylate oligomer, polyurethane resin, and photoinitiator [25,26]. The 3D printer mixed two base materials to form a broad range of intermediate materials with different strengths, which were predefined by the manufacturer. Six materials were selected to design and fabricate the required FGM plate for the experiment. Material 1 was the first base material VeroWhitePlus. Materials 2–6 were intermediate materials with graded strengths. The manufacturer only reported Young’s modulus of VeroWhitePlus as E = 2000–3000 MPa. In order to obtain the specific Young’s moduli of six materials, their tensile properties were tested according to the Standard Test Method of Plastics D638-14 made by the American Society for Testing Material (ASTM). The dog-bone-shaped specimen for tension test in the Standard is shown in Figure 2 [27].
The 3D models of the tension specimens were firstly built with commercial software Solidworks (sw2018, 2018, Concord, Massachusetts, America); then, with the models imported into the printing software Objet Studio, 3D objects could be automatically printed by 3D Printer. The six obtained specimens can be found in Figure 3. The tensile properties of specimens were tested by electronic universal testing machine. Stress-strain curves are shown in Figure 3. It can be found that the strength of material 1 was the highest and material 6 was the lowest. The strengths of material 2–5 were between those of material 1 and 6. The Young’s moduli of six printing materials were also measured and listed in Table 1. It is found that six Young’s moduli showed a well-graded change, which was in accordance with our expectation. In the following experiment, the specimens of FGM plates were designed and fabricated with using six printing materials.

3.2. Design and Fabrication of FGM Plate with a Circular Hole

On the basis of the theoretical model in Section 2, two FGM plates with the same geometrical sizes but different material distribution were designed, as shown in Figure 4. In terms of geometric design in Figure 4a, the width and thickness of the FGM plate were taken as 120 mm and 3 mm, respectively. The diameter of the center hole was 10 mm. There were five rings, and the width of each ring was 4 mm. The width of the plate was more than 10 times the diameter of the hole. According to the results shown in Figures 10 and 11 in the paper by Yang et al. [28], who discussed the effect of FGM plate size on the stress concentration, the plate can be approximately considered infinite in the case and it causes few errors, which can be neglected. Therefore, the following experimental results for the designed plates can be compared with the theoretical solution, which was derived based on the assumption of the infinite plate in Section 2.
In terms of material design, two different cases, in which Young’s modulus decreases and increases in the radial direction, were designed, as shown in Figure 4b. The materials in the five rings and outer plate were set in the order of materials 1-6 and 6-1, respectively. They are named the decreasing FGM plate and increasing FGM plate, respectively, throughout the paper. During the design of increasing FGM plate, we matched the desired optimal distribution of Young’s modulus given by Equation (13) as best we could, given the available materials.
In order to fabricate two FGM plates, 3D models were firstly built by Solidworks, according to the geometric sizes of two FGM plates. It is worth emphasizing that in order to set different material types in the software of 3D Printer, each part with different material in the plate should be separately built during the modeling process. The complete 3D model of FGM plate was obtained by assembling each separate part. The specimens of FGM plate were then fabricated, similar to the printing process of previous tension specimen. The obtained specimens are shown in Figure 5. According to different material colors in Figure 5a,b, it can be found that the material distributions in two FGM plates were in accordance with material design. Figure 5c,d shows the surface macromorphology of two FGM plates observed by microscope DSX 100 (Olympus Corporation, Tokyo, Japan). It is clear that the surface was smooth and there were no macroscopic cracks, inclusions, or other defects in two plates. Figure 5e,f shows the interface morphology observed by Carl Zeiss Axio Imager (Carl Zeiss AG, Oberkochen, Germany). It is found that two FGM plates had high density, good bonding interface, and forming quality.

3.3. Strain Measurement of FGM Plates

In this work, strain gages were utilized to determine the strain distribution along the radial direction in the FGM plate. The physical dimensions of strain gauges were 3.5 × 2.2 mm. Seven measuring points were chosen, and their radial positions on the x-axis are listed in Table 2. The points with numbers #1–5 were located in five rings, and points #6 and #7 were located in the outer plate. For the two-dimensional plane stress problems in the paper, stress calculation required knowledge of normal strain in both X and Y directions. Therefore, two strain gages were pasted parallel to the X and Y axis, respectively, for each measuring position. Moreover, in order to monitor and minimize the out-of-plane warping deformation of the plate after fabrication or any other possible out-of-plane bending during loading, a further seven groups of strain gages were mounted on the back faces of the plate corresponding to each measured position. Two strain gauges with the same radial position on two faces were connected to the half bridge circuit in the experiment.
The incremental loading method was employed in the tension test. The load was increased with the same scale as Δ F = 100 N each time and from 100 N up to 700 N in seven loads. The strains ε x and ε y at seven measured points were read out for each load level, as listed in Table 3 and Table 4. Since the measurement was made on the half bridge connection, the original measured strains in Table 3 and Table 4 are two times the real strain. The real strains in the paper were expressed as ε i j , where the superscript i represents the loading times and the j represents the number of measured points. It is worth noting that all materials in the whole FGM plate were within the range of linear elastic deformation under the maximum tensile load 700 N. In order to show the measured strains more visibly, the variations of the strains ε i j along the x direction were drawn, as in Figure 6 for two FGM plates under seven different loads. It can be seen that as the loads increase, the strains ε x i j and ε y i j at each measured point decrease or increase with equal magnitude. In particular, two strains ε y i j for the loads F = 600 N and 700 N fluctuated obviously at point #6 in the decreasing plate. The reason for fluctuations may be the instability of resistance strain gauge or the problem of bridge balance during measurement.

4. Results and Discussions

4.1. Theoretical Results

The theoretical solution of the stress field was firstly presented via the analytical method [15] before discussing the experimental results. According to the dimensions of printed FGM plates, we took the radius of circular hole as r 0   =   5   mm in the theoretical model and decomposed the FGM plate into five homogenous rings in the domain r   =   5 ~ 25   mm and an outer homogenous plate. The Young’s moduli E j ( j   =   1 , 2 , 5 ) in the rings and E 6 in outer plate were chosen as those listed in Table 1. Moreover, the remote loading stresses were taken as σ x = 0 , τ x y = 0 , and
σ y = F i / A ,
where the cross-sectional area of the plate is A   =   360   mm 2 . Since the variation of Poisson’s ratio has little influence on the stress distribution [15], Poisson’s ratio in the whole plate is assumed to be a constant ν ( r )   =   0.3 , which is the common value for materials.
The theoretical results of stress σ y in two FGM plates under seven load levels are shown in Figure 7. The solid lines represent the full field stress distributions. In order to make a more intuitive comparison with the experimental results, the theoretical results of stress at the measured points are marked with asterisks in Figure 7. It can be seen that the full-field stress in two FGM plates increased in a regular manner, with the load increasing from 100 N to 700 N. Due to the difference of Young’s modulus between each ring, the hoop stress had an obvious abrupt change at the interface of each ring. By comparing the results in Figure 7a,b, it can be found that the stresses σ y of two FGM plates differ small away from the hole but differ greatly near the hole. For the decreasing FGM plates, the stresses near the hole were very high and about 4–5 times of those in the far field. On the contrary, for the increasing FGM plates, the stresses near the hole were very low and even lower than those in the far field.

4.2. Experimental Results

Based on the experimental strains ε x i j and ε y i j in Table 3 and Table 4, the stresses σ y i j at seven measured points can be calculated according to the generalized Hooke’s law due to the tensile deformation within the range of linear elastic
σ y i j = E j 1 - ν 2 ( ε y i j + ν ε x i j ) ,
where E j is Young’s modulus corresponding to each measured point.
The variations of the experimental stress σ y along the x direction were drawn as in Figure 8 for two FGM plates under seven different loads. It can be seen that, as the load F i increased from 100 N to 700 N, the experimental stress at most positions in two plates increased by the same amount. In particular, there were certain fluctuations at the point of strain gage #6 in the decreasing plate in Figure 8a. It can be also found from Figure 8a that the stress value decreased rapidly away from the hole. The stresses near the hole were significantly higher than those in the far field. For the increasing FGM plate in Figure 8b, the stresses near the hole were lower than those in the far field, and the maximum stress appeared at the position of strain gage #5. The stress concentration did not occur at the edge of hole. The experimental results here were consistent with the theoretical analysis in Figure 7 and further prove an important conclusion proposed by Yang et al. [15]. That is, the decreasing Young’s modulus along the radial direction in the plate results in more severe stress concentration at the hole, while the increasing Young’s modulus can alleviate or even eliminate the stress concentration near the hole.

4.3. Comparison of Experimental and Theoretical Results

In order to make a comparison of experimental and theoretical results, their average stresses under seven loads were calculated. For the experimental results, we firstly calculated the strain increments Δ ε x j and Δ ε y j , corresponding to the loading increment Δ F = 100 N at each point, and then calculated their average strain increments Δ ε x - ave j and Δ ε y - ave j . Finally, the average stresses σ y - ave j can be derived at each point in the linear range
σ y - ave j = E j 1 - ν 2 ( Δ ε y - ave j + ν Δ ε x - ave j ) .
The comparison of experimental and theoretical results of dimensionless hoop stresses along x direction is shown in Figure 9. The influences of the radius of a circular hole and loads were eliminated in the x-coordinate and y-coordinate, respectively. The referring stress σ 0 in y-coordinate was taken as σ 0 = Δ F / A . The black solid curve in Figure 9 denotes the stress distribution of a homogeneous plate containing a circular hole under uniaxial tension, which is one of the classical solutions in the elastic theory [24,25,26,27,28,29]. It is well known that the SCF is 3 for the case. The blue and red curves represent the stress distributions in decreasing and increasing FGM plates, and the theoretical and experimental results of two FGM plates are denoted by solid and dash lines, respectively. It can be found that, among three types of plates, the peak stress or SCF was highest in the decreasing FGM plate but lowest in the increasing FGM plate. It can also be found that, although the overall varying trend of experimental and theoretical results was consistent, three experimental points only of 14 points lay on the theoretically calculated curves. The other experimental points were more or less out of the theoretical predictions.
The errors between the experimental and theoretical results of dimensionless stresses σ y ave / σ 0 were calculated and listed in Table 5. It can be observed that the error was more than 15% at measured point #1 and more than 20% at points #6 and #7 in both plates. For measured point #1, the error was mainly caused by the limitation of the strain gauge method. As we all know, one strain gage can only measure one “point”, and the measured strain was just the average value measured over the active area of the strain gage. The active area of the strain gage covered a steep gradient of stress concentrations at the edge of hole (point #1). Therefore, its accuracy was much less. The errors at points #6 and #7 were mainly caused by the size of chucks in the testing machine. The width of the chucks was smaller than that of the plate, as shown in Figure 8. Therefore, the chucks only held the middle part of the plate during loading. The marginal area of the plate was not loaded enough. Therefore, the experimental stresses at points #6 and #7 were much lower than the theoretical predictions. The errors at other points mainly came from some assumptions in theoretical calculation and the operation of the experimental process; for example, (1) neglect of the variation of Poisson’s ratio in the theoretical analysis, (2) errors in the measurement of Young’s modulus of printing materials, (3) human operations, including mounting of strain gauge, clamping of specimen, etc. With the above factors taken into account, it can be considered that the experimental results verify the theoretical analysis within the reliable accuracy of the test.

5. Conclusions

Based on the authors’ previous theoretical model, two kinds of FGM plates containing holes with radially increasing and decreasing Young’s modulus were designed and fabricated by means of modern multi-material 3D printing. The strains at different radial points in two FGM plates were experimentally measured using strain gages under uniaxial tension. The stresses were calculated within the range of linear elastic from the measured strains. The results were compared with the theoretical solutions derived on the basis of the authors’ previously theoretical model. It is found that the experimental results agree with the theoretical solutions. They both demonstrate that the stresses near the hole in the decreasing FGM plate were much higher than that of homogeneous plate, while the stress concentration reduced obviously in the increasing FGM plate. The successful experiment and fabrication of the specimens in the paper are of great significance to the application of FGM in reducing stress concentration in the future.

Author Contributions

Conceptualization, Q.Y. and Y.L.; methodology, H.C. and Y.T.; formal analysis, Q.Y. and X.C.; investigation, X.C. and Y.L.; data curation, Y.T. and H.C.; writing, Q.Y. and H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NATURAL SCIENCE RESEARCH PROJECT OF JIANGSU HIGHER EDUCATION INSTITUTIONS, grant numbers 19KJA150007 and 17KJA460003, and PROJECT OF “SIX TALENT PEAKS” IN JIANGSU PROVINCE, grant number GDZB-014.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available on request to the corresponding author.

Acknowledgments

Q.Y. is grateful to the support of Funding of Jiangsu Overseas Visiting Scholar Program for University Prominent Young and Middle-Aged Teachers and Presidents. Laboratory Technician Feng He skillfully mounted the strain gages.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tsukamoto, H. Design of functionally graded thermal barrier coatings based on a nonlinear micromechanical approach. Comput. Mater. Sci. 2010, 50, 429–436. [Google Scholar] [CrossRef]
  2. Noda, N.; Guo, L.C. Thermal shock analysis for a functionally graded plate with a surface crack. Acta Mech. 2008, 195, 157–166. [Google Scholar] [CrossRef]
  3. Yang, Q.; Gao, C.F. Reduction of the stress concentration around an elliptic hole by using a functionally graded layer. Acta Mech. 2016, 227, 2427–2437. [Google Scholar] [CrossRef]
  4. Goyat, V.; Verma, S.; Garg, R.K. On the reduction of stress concentration factor in an infinite panel using different radial functionally graded materials. Int. J. Mater. Prod. Technol. 2018, 57, 109–131. [Google Scholar] [CrossRef]
  5. Venkataraman, S.; Haftka, R.T.; Rapoff, A.J. Structural optimization using biological variables to help understand how bones design holes. Struct. Multidiscip. Optim. 2003, 25, 19–34. [Google Scholar] [CrossRef]
  6. Huang, J.; Venkataraman, S.; Rapoff, A.J.; Haftka, R.T. Optimization of axisymmetric elastic modulus distributions around a hole for increased strength. Struct. Multidiscip. Optim. 2003, 25, 225–236. [Google Scholar] [CrossRef]
  7. Huang, J.; Rapoff, A.J. Optimization design of plates with holes by mimicking bones through nonaxisymmetric functionally graded material. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2003, 217, 23–27. [Google Scholar] [CrossRef]
  8. Kubair, D.V.; Bhanu-Chandar, B. Stress concentration factor due to a circular hole in functionally graded panels under uniaxial tension. Int. J. Mech. Sci. 2008, 50, 732–742. [Google Scholar] [CrossRef]
  9. Wang, W.; Yuan, H.; Li, X.; Shi, P. Stress concentration and damage factor due to central elliptical hole in functionally graded panels subjected to uniform tensile traction. Materials 2019, 12, 422. [Google Scholar] [CrossRef] [Green Version]
  10. Goyat, V.; Verma, S.; Garg, R.K. Reduction of stress concentration for a rounded rectangular hole by using a functionally graded material layer. Acta Mech. 2017, 228, 3695–3707. [Google Scholar] [CrossRef]
  11. Goyat, V.; Verma, S.; Garg, R.K. Reduction in stress concentration around a pair of circular holes with functionally graded material layer. Acta Mech. 2018, 229, 1045–1060. [Google Scholar] [CrossRef]
  12. Berezhnoi, D.V.; Chickrin, D.E.; Gajnulina, L.R.; Kokunin, P.A.; Kurchatov, E.Y.; Sozutov, I.G. Calculation of stresses in superflywheel. Contemp. Eng. Sci. 2016, 9, 1591–1598. [Google Scholar] [CrossRef]
  13. Berezhnoi, D.V.; Sachenkov, A.A.; Xu, Y. Investigation of the effect of elastic-plastic loading modes in the manufacture of metal flywheels on their performance properties. IOP Conf. Ser. Mater. Sci. Eng. 2020, 709, 044018. [Google Scholar] [CrossRef]
  14. Zhang, X.Z.; Kitipornchai, S.; Liew, K.M.; Lim, C.W.; Peng, L.X. Thermal stresses around a circular hole in a functionally graded plate. J. Therm. Stresses 2003, 26, 379–390. [Google Scholar] [CrossRef]
  15. Yang, Q.; Gao, C.F.; Chen, W.T. Stress analysis of a functional graded material plate with a circular hole. Arch. Appl. Mech. 2010, 80, 895–907. [Google Scholar] [CrossRef]
  16. Yang, Q.; Gao, C.F. Dynamic stress analysis of a functionally graded material plate with a circular hole. Meccanica 2013, 48, 91–101. [Google Scholar] [CrossRef]
  17. Mohammadi, M.; Dryden, J.R.; Jiang, L.Y. Stress concentration around a hole in a radially inhomogeneous plate. Int. J. Solids Struct. 2011, 48, 483–491. [Google Scholar] [CrossRef]
  18. Sburlati, R.; Atashipour, S.R.; Atashipour, S.A. Reduction of the stress concentration factor in a homogeneous panel with hole by using a functionally graded layer. Compos. Part B 2014, 61, 99–109. [Google Scholar] [CrossRef]
  19. Kubair, D.V. Stress concentration factors and stress-gradients due to circular holes in radially functionally graded panels subjected to anti-plane shear loading. Acta Mech. 2013, 224, 2845–2862. [Google Scholar] [CrossRef]
  20. Kubair, D.V. Stress concentration factor in functionally graded plates with circular holes subjected to anti-plane shear loading. J. Elast. 2014, 114, 179–196. [Google Scholar] [CrossRef]
  21. Nie, G.J.; Zhong, Z.; Batra, R.C. Material tailoring for reducing stress concentration factor at a circular hole in a functionally graded material (FGM) panel. Compos. Struct. 2018, 205, 49–57. [Google Scholar] [CrossRef]
  22. Nie, G.J.; Batra, R.C. Reducing stress concentration factor by strengthening circular hole with functionally graded incompressible material layer. Thin Wall Struct. 2019, 144, 106223. [Google Scholar] [CrossRef]
  23. Buskirk, S.R.; Venkataraman, S.; Ifju, P.G.; Rapoff, A.J. Functionally graded biomimetic plate with hole. In Proceedings of the AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Denver, CO, USA, 22–25 April 2002; p. 1330. [Google Scholar]
  24. Muskhelishvili, N.I. Some Basic Problem of Mathematical Theory of Elasticity; Noordhoff: Groningen, The Netherlands, 1975. [Google Scholar]
  25. Ge, Q.; Dunn, C.K.; Qi, H.J.; Dunn, M.L. Active origami by 4D printing. Smart Mater. Struct. 2014, 23, 094007. [Google Scholar] [CrossRef]
  26. Teoh, J.E.M.; An, J.; Chua, C.K.; Lv, M.; Krishnasamy, V.; Liu, Y. Hierarchically selfmorphing structure through 4D printing. Virtual Phys. Prototyp. 2017, 12, 61–68. [Google Scholar] [CrossRef]
  27. American Society of Testing and Materials. D638-14. In Standard Test Method for Tensile Properties of Plastics; American Society of Testing and Materials: West Conshohocken, PA, USA, 2014. [Google Scholar]
  28. Yang, Q.; Gao, C.F.; Chen, W.T. Stress concentration in a finite functionally graded material plate. Sci. China Phys. Mech. 2012, 55, 1263–1271. [Google Scholar] [CrossRef]
  29. Savin, G.N. Stress Concentration around Holes; Pergamon Press: London, UK, 1961. [Google Scholar]
Figure 1. The theoretical model of a FGM plate with a circular hole developed in the authors’ previous work [15].
Figure 1. The theoretical model of a FGM plate with a circular hole developed in the authors’ previous work [15].
Materials 14 07845 g001
Figure 2. Dog-bone-shaped specimen for tension test (ASTM D638-14).
Figure 2. Dog-bone-shaped specimen for tension test (ASTM D638-14).
Materials 14 07845 g002
Figure 3. Stress-strain curves for six printing materials in tension.
Figure 3. Stress-strain curves for six printing materials in tension.
Materials 14 07845 g003
Figure 4. Design for FGM plates with a circular hole: (a) geometric design (mm), (b) material design.
Figure 4. Design for FGM plates with a circular hole: (a) geometric design (mm), (b) material design.
Materials 14 07845 g004
Figure 5. Specimens of two FGM plates fabricated by 3D printing: (a) photo of decreasing FGM plate, (b) photo of increasing FGM plate, (c) macromorphology of decreasing FGM plate, (d) macromorphology of increasing FGM plate, (e) interface morphology of decreasing FGM plate, (f) interface morphology of increasing FGM plate.
Figure 5. Specimens of two FGM plates fabricated by 3D printing: (a) photo of decreasing FGM plate, (b) photo of increasing FGM plate, (c) macromorphology of decreasing FGM plate, (d) macromorphology of increasing FGM plate, (e) interface morphology of decreasing FGM plate, (f) interface morphology of increasing FGM plate.
Materials 14 07845 g005
Figure 6. The measured strains ε i j along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Figure 6. The measured strains ε i j along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Materials 14 07845 g006
Figure 7. Theoretical results of stresses σ y along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Figure 7. Theoretical results of stresses σ y along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Materials 14 07845 g007
Figure 8. Experimental result of stresses σ y along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Figure 8. Experimental result of stresses σ y along the x-axis in two FGM plates: (a) decreasing FGM plate, (b) increasing FGM plate.
Materials 14 07845 g008
Figure 9. Comparison of experimental and theoretical results of dimensionless stresses σ y ave / σ 0 along the x direction (the blue curve shows the decreasing FGM plate and red curve shows the increasing FGM plate).
Figure 9. Comparison of experimental and theoretical results of dimensionless stresses σ y ave / σ 0 along the x direction (the blue curve shows the decreasing FGM plate and red curve shows the increasing FGM plate).
Materials 14 07845 g009
Table 1. Elastic modulus E j of six printing materials (MPa).
Table 1. Elastic modulus E j of six printing materials (MPa).
Material 1Material 2Material 3Material 4Material 5Material 6
2154188115761174880445
Table 2. Radial positions of seven measured points on the x-axis (mm).
Table 2. Radial positions of seven measured points on the x-axis (mm).
Point #1Point #2Point #3Point #4Point #5Point #6Point #7
7111519233545
Table 3. Measured strains 2 ε i j ( 10 - 6 ) in the decreasing FGM plate under different loads F i .
Table 3. Measured strains 2 ε i j ( 10 - 6 ) in the decreasing FGM plate under different loads F i .
Load  F i Point #1Point #2Point #3Point #4Point #5Point #6Point #7
2 ε y i 1 2 ε x i 1 2 ε y i 2 2 ε x i 2 2 ε y i 3 2 ε x i 3 2 ε y i 4 2 ε x i 4 2 ε y i 5 2 ε x i 5 2 ε y i 6 2 ε x i 6 2 ε y i 7 2 ε x i 7
100 N484.2−59.2370.2−38.4422.7−75.1379.2−81.4502.5−153.7309.7−151.7263.4−142.9
200 N997.4−174772.0−96.2714.5−145.3778.8−163.9951.3−309.5605.4−287.9602.5−289.4
300 N1530.1−285.31177.3−155.21173.5−217.81192.3−239.41389.2−454.3910.3−414.1955.5−431.1
400 N2151.3−389.41662.5−211.41698.7−285.81669.3−308.81854.2−588.71235.7−529.51387.3−562.4
500 N2746.2−490.42129.4−256.92193.6−345.72122.2−379.72281.4−730.61610.7−6571812.1−694.1
600 N3290.3−602.52605.7−311.62601.2−415.62589.7−4472712.9−870.43200.5−813.22239.6−843.8
700 N3828.0−717.53043.2−361.92891.7−476.52989.0−510.23072.8−999.33619.9−947.42604.9−986.1
Table 4. Measured strains 2 ε i j ( 10 6 ) in the increasing FGM plate under different loads F i .
Table 4. Measured strains 2 ε i j ( 10 6 ) in the increasing FGM plate under different loads F i .
Load   F i Point #1Point #2Point #3Point #4Point #5Point #6Point #7
2 ε y i 1 2 ε x i 1 2 ε y i 2 2 ε x i 2 2 ε y i 3 2 ε x i 3 2 ε y i 4 2 ε x i 4 2 ε y i 5 2 ε x i 5 2 ε y i 6 2 ε x i 6 2 ε y i 7 2 ε x i 7
100 N659.4−78.1496.5−58.8435.2−49.2388−52.7341.6−74.2198−33.7127.3−39.9
200 N1238.6−121.9971.4−116.1872.5−92.7746.3−96.9688.3−120.9410.7−59.3286.4−69.2
300 N1795−170.11455.7−1741346.5−141.21157−148.71055−165.9640−88.2464.9−98.3
400 N2376.8−212.31948.9−234.91824.5−189.21568.9−194.51446.1−214.6860−113.7651.2−129.9
500 N2965.8−258.12407.9−294.32294.7−235.41978.4−242.71825.6−260.91074.6−139.2833.5−158.5
600 N3567.6−300.22876−354.22764.9−283.32380.8−292.72210.6−307.61285−169.51014−187.9
700 N4238.4−340.33293.6−412.33187.2−327.92740.5−3432557.1−349.31476.9−194.61182.6−219.4
Table 5. The errors for experimental and theoretical results of dimensionless stresses σ y ave / σ 0 .
Table 5. The errors for experimental and theoretical results of dimensionless stresses σ y ave / σ 0 .
Point NumberTheoretical ResultExperimental ResultError
Decreasing PlateIncreasing PlateDecreasing PlateIncreasing PlateDecreasing PlateIncreasing Plate
Point #12.6310.6822.2340.514−15.1%−24.6%
Point #21.6350.9211.5970.781−2.3%−15.2%
Point #31.3420.9951.2201.033−9.1%3.8%
Point #41.1011.1460.9601.177−12.8%2.7%
Point #50.9541.2070.6721.323−29.6%9.6%
Point #60.8351.1140.4510.874−46.0%−21.5%
Point #70.9211.0580.3060.711−66.8%−32.8%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Q.; Cao, H.; Tang, Y.; Li, Y.; Chen, X. Experimental Investigation of Stress Distributions in 3D Printed Graded Plates with a Circular Hole. Materials 2021, 14, 7845. https://doi.org/10.3390/ma14247845

AMA Style

Yang Q, Cao H, Tang Y, Li Y, Chen X. Experimental Investigation of Stress Distributions in 3D Printed Graded Plates with a Circular Hole. Materials. 2021; 14(24):7845. https://doi.org/10.3390/ma14247845

Chicago/Turabian Style

Yang, Quanquan, He Cao, Youcheng Tang, Yun Li, and Xiaogang Chen. 2021. "Experimental Investigation of Stress Distributions in 3D Printed Graded Plates with a Circular Hole" Materials 14, no. 24: 7845. https://doi.org/10.3390/ma14247845

APA Style

Yang, Q., Cao, H., Tang, Y., Li, Y., & Chen, X. (2021). Experimental Investigation of Stress Distributions in 3D Printed Graded Plates with a Circular Hole. Materials, 14(24), 7845. https://doi.org/10.3390/ma14247845

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop