Next Article in Journal
Hybrid Geopolymeric Foams for the Removal of Metallic Ions from Aqueous Waste Solutions
Next Article in Special Issue
On the Room-Temperature Creep Behavior and Its Correlation with Length Scale of a LiTaO3 Single Crystal by Spherical Nanoindentation
Previous Article in Journal
Electrophoretic Deposition of Quantum Dots and Characterisation of Composites
Previous Article in Special Issue
TiN-Nanoparticulate-Reinforced ZrO2 for Electrical Discharge Machining
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Functionally Step-Variable Graded Plate Under In-Plane Compression

by
Leszek Czechowski
* and
Zbigniew Kołakowski
Department of Strength of Materials, Lodz University of Technology, 90-924 Lodz, Poland
*
Author to whom correspondence should be addressed.
Materials 2019, 12(24), 4090; https://doi.org/10.3390/ma12244090
Submission received: 16 October 2019 / Revised: 30 November 2019 / Accepted: 5 December 2019 / Published: 7 December 2019
(This article belongs to the Special Issue Mechanical Properties and Applications of Advanced Ceramics)

Abstract

:
A study of the pre- and post-buckling state of square plates built from functionally graded materials (FGMs) and pure ceramics is presented. In contrast to the theoretical approach, the structure under consideration contains a finite number of layers with a step-variable change in mechanical properties across the thickness. An influence of ceramics content on a wall and a number of finite layers of the step-variable FGM on the buckling and post-critical state was scrutinized. The problem was solved using the finite element method and the asymptotic nonlinear Koiter’s theory. The investigations were conducted for several boundary conditions and material distributions to assess the behavior of the plate and to compare critical forces and post-critical equilibrium paths.

1. Introduction

Functionally graded materials (FGMs) are treated as modern materials which due to their varying properties for thickness or other dimensions are assigned to work in ultrahigh hard environments under high and other temperature fields. Gradual changes in the volume fraction of constituents and the non-homogenous structure provide theoretically continuous graded macroscopic properties (hardness, wear resistance, thermal conductivity, specific heat, mass density, etc.). At present, well-known techniques for the manufacturing of FGMs have been strongly developed, however taking into account different properties or different combinations of the components used, none of them can be found to be excellent and devoid of defects [1]. While taking a look at methods of FGM fabrication, one can distinguish: gas-based methods (chemical vapor deposition, thermal spray, and the surface reaction process), liquid phase processes (chemical solution deposition, laser deposition, and electro-chemical gradation) or solid-phase processes (spark plasma sintering, powder metallurgy) [1]. It should be mentioned that actual FGMs differ from one another in contrast to the theoretical descriptions of these materials. Thus, it is impossible to obtain a perfect structure with fluently varying properties. Indeed, moderately adequate structures of FGMs with appropriate constituents can be attained using some techniques, but many methods enable producing mere multilayer (several-layer) moderately thick structures of FGMs (using spark plasma sintering, powder metallurgy, etc.). The basic concept of FGMs was presented by Niino and the other studies referring to FGMs can be found in references [2,3,4,5]. Further research on FGMs showed that many works were devoted to an analysis of manufacturing processing and its influence on various properties [6,7,8,9]. On the other hand, in many papers based on original assumptions, perfect structures were often analyzed (with theoretical and assumed material distributions as well as assumed material properties). Those works concerned theoretical studies on the stability of thin-walled structures under mechanical or/and thermal load (plates, boxes, etc.) or investigations of the frequency of natural vibrations. Such results of buckling and post-buckling of the FGM plate under mechanical loads were included in references [10,11,12,13,14,15,16,17] or were developed for plate structures in reference [18,19,20]. Analyses of natural vibrations of FGM structures were reported in references [21,22,23]. Other works were focused on the investigation of the structure behavior under static or dynamic thermal loads [24,25,26,27]. Kumar et al. in reference [28] analyzed the buckling and post-buckling state of an FGM plate using two new higher-order transverse shear deformation theories (NHSDTs). Trabelsi et al. in reference [29] studied the response of FGM shell structures (plates and cylindrical shells separately) due to thermal load with the first-order shear deformation theory (FOSDT). In reference [30], a study of the thick simply supported FGM plate under bending was presented within the displacement potential function (DPF). Nguyen et al. in their papers [31,32] proposed analytical modeling of thin-walled open-section beams made of functionally graded materials on the basis of Vlasov’s assumptions. Moreover, the same authors also investigated functionally graded open-section beams with different types of material distributions by means of a two-node beam element with 14 degrees of freedom in reference [33]. Xu et al. [34] analyzed elastoplastic buckling behaviors of rectangular plates built of functionally graded materials with a homogenization method of the Tamura-Tomota-Ozawa model. Besides studies on strictly functionally graded materials, recently more and more analyses devoted to printed functional materials have been developed which can be found in references [35,36,37,38], and elsewhere.
As assessed in the literature, almost all of these papers were based on theoretical or mathematical mechanical properties of perfect FGM structures. Because it seems impossible to achieve an idealized material distribution (on the basis of the literature survey), the authors of the present paper conducted investigations and showed the results for a plate with step-variable gradation of mechanical properties and with a finite number of layers across the wall. Moreover, a combination of the FGM with additional components (in the presented case, additional thickness of ceramics was studied) was taken into consideration. The non-linear problem of stability was solved with the semi-analytical method (SAM) based on non-linear Koiter’s theory [27] and with the finite element method (FEM) code ANSYS® [39]. The full Green’s strain tensor, the second Piola-Kirchhoff’s stress tensor, and the transition matrix using Godunov’s orthogonalization were used in the description of the problem. In reference [17], on the basis of Koiter’s theory, FGM plates have non-symmetric stable post-buckling equilibrium paths. This feature explains differences in the plate response dependence on the imperfection sign. An FGM plate has a non-trivial coupling matrix B and the coupling between extensional and bending deformations exists as is in the case of non-symmetric laminated plates. On the basis of both of the above-mentioned methods, the present paper reveals the results of influence on the critical and post-critical state of plates made of an alumina-FGM with a finite number of layers (a total of five to 15 layers were assumed). Furthermore, several boundary conditions of the edge support and different contents of alumina in relation to the FGM were taken into account. Considered variants of material distributions were assumed to reflect real FGM.

2. Problem Description

The object of investigation was a square plate (a = b = 1 m) subjected to the mechanical compression load, as shown in Figure 1a.
The length and the total thickness of the plate was equal to 200 mm and 2 mm ( t t ), respectively.
The thickness of pure ceramics t c ranges from 0.2 mm to 1 mm but the thickness of FGM t F G M is equal from 1.0 mm to 1.8 mm (the total thickness of wall is a sum t t = t F G M + t c ). The description of considered variants were posted in Table 1. The total thickness of a plate comprises the FGM (Al-Al203) and the pure ceramics Al203 (See Figure 1b).
To solve the problem, two methods were employed: an analytical-numerical method based on the asymptotic Koiter’s approach [15,16,17] and a numerical one based on the finite element method [39]. To verify an influence of the material distribution (see Table 1) for Var_1 from 5 up to 15 (5, 7, 11, 15), layers of the FGM with a composition shown in Figure 2 were considered. In the case of Var_2-Var_5, 11 layers of the FGM were assumed. Basic material properties for the considered constituents have been listed in Table 2, but in the case of a composition of different materials in a dependence on their content, the mixture law was applied.

2.1. FE Model

Numerical simulations based on the finite element method were conducted with the ANSYS 18.2® software [39]. To generate an adequate numerical model, an 8-node 281 shell element was assumed. The plate was divided into 10 k finite elements (100 elements along the edge—Figure 3). A nonlinear analysis for large deflections was performed on the basis of Green-Lagrangian equations. During computations, nonlinear calculations were conducted in accordance with the Newton-Raphson algorithm. The critical- and post-critical state of the plate with different boundary conditions was analyzed (See Table 3). The initial deflection of the plate in all cases was assumed to be 0.01tt, which referred to the first buckling mode.

2.2. Koiter’s Asymptotic Approach

The equilibrium equations for FGM plate structures can be written as [15,16,17,18]:
( 1 σ σ r ) ζ r + a p q r ζ p ζ q + b r r r r ζ r 3 σ σ r ζ r * + = 0    for r = 1 , J
where: σ r is the critical stress of the r-th buckling mode, ζ r is the dimensionless amplitude of the r-th buckling mode, ζ r * is the dimensionless amplitude of the initial imperfections related to the r-th buckling mode, σ is the compressive stress, and a p q r and b r r r r are the coefficients, respectively. The range of indices are p , q , r is from 1 to J , where J is the number of interacting modes. The summation was determined based on the repeated indices. For the case of the uncoupled buckling mode (i.e., for the one-mode approach), J = 1 has to be satisfied. The first order coefficients (i.e., a p q r ) were found with the analytical-numerical method based on Koiter’s theory [1518,27] etc. The second order coefficients (i.e., b r r r r ) were calculated with the semi-analytical method (SAM) [16,17]. When using that method, it was necessary to determine approximate values of b r r r r on the basis of the linear buckling problem. The co-author of SAM is also the co-author of this article. In the present work, a SAM one-modal approach has been applied ( J = 1 in Equation (1)). This means one degree of freedom in nonlinear analysis was taken into account. The imperfection of the plate was assumed as in the numerical model: ζ 1 * = w 0 / t t = 1.0 . The finite element method allows us to simulate a complex numerical model but boundary conditions and element types also play a major role. In comparison to the SAM included in the present paper, a single calculation based on FEM requires between one hour and a couple of days to compute. Using the semi-analytical method, one can obtain the results within 5–15 min, which can provide a major benefit. Moreover, differences in results between two methods amount to a few percent (referring to critical forces) or to several dozen percent (compared to the post-buckling paths), at most.

3. Results and Discussion

3.1. Buckling Forces

In this subsection, critical forces F c r are determined within both of the methods. The critical forces for the first buckling mode related to a number of layers are given in Table 4. It was noticed that with an increasing number of layers and a step-variable gradation (according to distributions as in Figure 2), the critical forces grew slightly (0.2–0.3% at most). This trend was also retained in the case of the SAM. In addition, the SAM for five layers showed the greatest discrepancy in the critical forces, which almost disappeared for more layers (even below 0.1%). The critical forces for the clamped plate were only estimated with the FEM. Taking into consideration an increase in the ceramics thickness (changing from 0.2 mm to 1 mm), the critical forces grew as well (Table 5).
This is caused by a significantly higher Young’s modulus and a growing thickness of ceramics. Nevertheless, a higher modulus of ceramics than the modulus of aluminium results in a non-symmetric distribution of the force with regard to the neutral and geometric axis, which can faster yield critical loads. Besides, it should be mentioned that very close results were achieved for both of the methods applied. As expected, the clamped plate can withstand more than twice as high load as in the case of the simply supported plate.

3.2. Post-Buckling Behavior of the Plate

In this subsection, the results concern an assessment of the plate behavior in the post-buckling state with reference to different parameters. Figure 4a,b depict a deflection in the center of the plate vs. the static compression load for the SSSS plate (in the first stage, one half-wave was noticed). The curves are plotted for both the methods under consideration. Figure 4a presents an influence of the layer number on the plate stability (from five to 15 layers with a step-variable gradation of material properties according to Figure 2). In the case of a change in the number of layers, any significant difference cannot be seen in the obtained curves. It seems that the overall bending stiffness for the considered cases provides the same effect. In a comparison of the results based on the two methods, a major discrepancy can be noticed. First of all, the FEM results provide a considerable deflection of the plate before the first critical force is reached. In contrast, in the SAM, a deflection appears just in the vicinity only and grows very fast (compared to the FEM curves). Moreover, in the SAM results, a change in the buckling mode was not observed after applying a higher load, which takes place in the FEM analysis in contrast. Furthermore, deflections in the SAM are greater than deflections in the FEM. The same response of the plate was seen for different variants (Figure 4b) by increasing the thickness of alumina and decreasing the thickness of the FGM. Nevertheless, the trends of curves seem to be very logical because with an increase in ceramics, the deflections of the plate slightly drop—this is due to higher resistance. Apart from these aspects, a jump in the buckling mode occurs on the diagrams at higher loads. Minor differences between the results of the two methods were observed for the SCSC plate (clamped unloaded edges—Figure 5a,b) but the plate analyzed with the FEM deflected faster than in the SAM.
The curves obtained with the FEM are almost identical regardless of a number of layers (Figure 5a), or remain close to each other as far as the considered variants (Var_1-Var_5) are concerned. In the case of the FEM analysis, some shifts of the curves are seen due to a growth in the number of layers. This occurs for the force approximately three times greater than the critical force. Taking into account different variants with an increase of ceramics content, a slightly higher resistance of the plate was observed (as in the previous case). However, for all variants, both for the SSSS plate or for the SCSC plate, a change in buckling modes in the plate occurs in the FEM analysis (in contrast, it does not happen in the SAM). The behavior for the plate clamped on all its edges (CCCC plate) based on the FEM alone was investigated (Figure 6a,b). In the first stage of the plate deflection (Figure 6a), the curves coincide with each other but under the load equal to double critical loads, the plots diverge slightly. For this case, a change in the buckling mode from one half-wave to two half-waves was noticed. For all the cases under consideration (Figure 7a,b), the same similarity can be seen. First of all, the deflection begins at a considerably lower force before the critical load is reached. Secondly, as Figure 7b depicts (the static load was referred to the critical load F c r for the given boundary conditions), for the load corresponding to the three-fold critical force, curves showing the middle deflection are almost the same. Afterwards, the further behavior of the plate results perhaps from the type of support available. The maps of plate deformation (displacements in nodes in direction of Z-axis) for three different boundary conditions and for VAR_1 were presented in Table 6. The color “blue” represents maximum values of displacements in the opposite direction to the Z-axis, but the color “red” denotes the maximum value of displacements in the same direction as the Z-axis (See Figure 3a). For the SSSS plate, a natural deflection (in the middle of the plate) was observed. With an increase in the compression force, the deflection was growing along unloaded edges until two half-waves appeared (under the load corresponding to the approximately six-fold critical force). This buckling mode in the plate did not last for a long time because a fast jump of deflection in three half-waves followed rapidly. In the center, the plate changed the direction of the maximal deflection. Afterwards, the buckling mode of three half-waves was continued until the maximal load was reached. In the case of the SCSC plate, more longitudinal deflection from the beginning due to clamped unloaded edges can be seen. The plate behaves very similarly to the SSSS plate. A change in the buckling mode from one half-wave to two half-waves was noted as well (at this point, the static load was equal to triple critical forces). In the case of the CCCC plate, the deflection is concentrated in the middle of the structure under analysis. In general, the deflection changes from the round shape to an elliptical one along the loaded edges. Finally, under a higher load, one half-wave changes violently into three half-waves.

4. Summary

An analytical-numerical and numerical analysis of the ceramics-FGM plate subjected to static load was performed. The functionally step-variable plate has a finite number of layers to reflect a real FGM structure. In addition, an influence of thickness of the ceramic wall on the overall stability of the plate was examined. The problem was solved with two methods, namely: the finite element method and the nonlinear asymptotic Koiter’s method. Moreover, the analysis was extended to consider several boundary conditions of the plate. On the basis of the achieved results, it was concluded that:
  • The gradation of layers in the number of plates from five to 15 revealed a slight influence on the plate stability because the obtained curves in both the methods ran almost identically. It indicates that the overall bending stiffness remains on a comparable level, though a non-symmetrical distribution of normal forces with respect to the neutral axis of the plate exists;
  • in a comparison of critical forces based on the two applied methods, a sufficiently good agreement was achieved (a few percent difference, at most). The SAM gave slightly lower values;
  • in the cases under consideration, the growth in ceramics thickness (from 0.2 mm to 1 mm) played an insignificant role in post-buckling paths. Indeed, the differences in curves are visible but they differ only slightly from one another;
  • a higher discrepancy could be seen by comparing the behavior of the plate obtained using the two methods. Firstly, in the FEM analysis, the plate deflects earlier than in the SAM analysis, but after exceeding the critical load, the SAM indicates a larger deflection. In addition, there is a change in the buckling mode during the plate compression. In contrast to the SAM, the FEM reveals a transformation of the defection function from one half-wave to two or three half-waves;
  • when analyzing the curves obtained for three different boundary conditions, a similarity in the plate deflection up to three-fold overloads was noticed if static loads were referred to their critical buckling loads (see Figure 7b).

Author Contributions

Z.K., L.C. developed this study. L.C. performed the numerical simulations. Z.K. carried out the analytical-numerical calculations. L.C. prepared the first version of this manuscript. Both the authors participated in the discussion and the revision of this manuscript.

Funding

This research was funded by the National Science Centre of Poland (UMO-2017/25/B/ST8/01046).

Acknowledgments

The investigations were carried out thanks to the financial support of the National Science Centre of Poland (UMO-2017/25/B/ST8/01046).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Naebe, M.; Shirvanimoghaddam, K. Functionally graded materials: A review of fabrication and properties. Appl. Mater. Today 2016, 5, 223–245. [Google Scholar] [CrossRef]
  2. Koizumi, M. FGM activities in Japan. Compos. Part B 1997, 28, 1–4. [Google Scholar] [CrossRef]
  3. Koizumi, M.; Niino, M. Overview of FGM Research in Japan. MRS Bull. 1995, 20, 19–21. [Google Scholar] [CrossRef]
  4. Kayser, W.A.; Ilschner, B. FGM research activities in Europe. MRS Bull. 1995, 20, 22–26. [Google Scholar] [CrossRef]
  5. Sasaki, M.; Wang, Y.; Hirano, T.; Hirai, T. Design of SiC/C functionally gradient material and its preparation by chemical vapor deposition. J. Ceram. Soc. Jpn. 1989, 97, 539–543. [Google Scholar] [CrossRef] [Green Version]
  6. Rubio, W.M.; Paulino, G.H.; Silva, E.C.N. Analysis, manufacture and characterization of Ni/Cu functionally graded structures. Mater. Des. 2012, 41, 255–265. [Google Scholar] [CrossRef]
  7. Cai, Y.; Cheng, Y.; Yin, H.; Yin, X.; Tian, Y.; Chen, J.; Wang, N. Preparation and mechanical properties of Ti3SiC2/SiC functionally graded materials. Ceram. Int. 2017, 43, 6648–6658. [Google Scholar] [CrossRef]
  8. Hangai, Y.; Morita, T.; Utsunomiya, T. Functionally graded aluminium foam consisting of dissimilar aluminium alloys fabricated by sintering and dissolution process. Mater. Sci. Eng. A 2017, 696, 544–551. [Google Scholar] [CrossRef]
  9. Liu, H.; Wang, Z.; Wang, Y.; Li, H. Effect of technological parameters on the process performance of pure Al2O3 layer of Ni–Al2O3 FGMs by self-induced EDM. Int. J. Adv. Manuf. Technol. 2017, 90, 3633–3641. [Google Scholar] [CrossRef]
  10. Shi, G. A new simple third-order shear deformation theory of plates. Int. J. Solids Struct. 2007, 44, 4399–4417. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, J.; Shen, H.S. Nonlinear bending analysis of shear deformable functionally graded plates subjected to thermo-mechanical loads under various boundary conditions. Compos. Part B Eng. 2003, 34, 103–115. [Google Scholar] [CrossRef]
  12. Samsam Shariat, B.A.; Eslami, M.R. Thermal buckling of imperfect functionally graded plates. Int. J. Solids Struct. 2006, 43, 4082–4096. [Google Scholar] [CrossRef] [Green Version]
  13. Yin, S.H.; Hale, J.S.; Yu, T.T.; Bui, Q.T.; Bordas, S.P.A. Isogeometric locking-free plate element: A simple first order shear deformation theory for functionally graded plates. Compos. Struct. 2015, 118, 121–138. [Google Scholar] [CrossRef] [Green Version]
  14. Bhardwaj, G.; Singh, I.V.; Mishra, B.K.; Bui, Q.T. Numerical simulation of functionally graded cracked plates using NURBS based XIGA under different loads and boundary conditions. Compos. Struct. 2015, 126, 347–359. [Google Scholar] [CrossRef]
  15. Kolakowski, Z.; Mania, R.J. Dynamic response of thin FG plates with a static unsymmetrical stable postbuckling path. Thin Walled Struct. 2015, 86, 10–17. [Google Scholar] [CrossRef]
  16. Kolakowski, Z.; Mania, J.R. Semi-analytical method versus the FEM for analysis of the local post-buckling. Compos. Struct. 2013, 97, 99–106. [Google Scholar] [CrossRef]
  17. Kolakowski, Z.; Mania, R.J.; Grudziecki, J. Local nonsymmetric post-buckling equilibrium path in thin FGM plate. Eksploat. Niezawodn. Maint. Reliab. 2015, 17, 135–142. [Google Scholar] [CrossRef]
  18. Czechowski, L.; Kolakowski, Z. The study of buckling and post-buckling of a step-variable FGM box. Materials 2019, 16, 918. [Google Scholar] [CrossRef] [Green Version]
  19. Kołakowski, Z.; Teter, A. Static interactive buckling of functionally graded columns with closed cross-sections subjected to axial compression. Compos. Struct. 2015, 123, 257–262. [Google Scholar] [CrossRef]
  20. Kołakowski, Z.; Teter, A. Load carrying capacity of functionally graded columns with open cross-sections under static compression. Compos. Struct. 2015, 129, 1–7. [Google Scholar] [CrossRef]
  21. Huang, X.L.; Shen, H.S. Nonlinear vibration and dynamic response of functionally graded plates in thermal environments. Int. J. Solids Struct. 2004, 41, 2403–2427. [Google Scholar] [CrossRef]
  22. Wattanasakulpong, N.; Prusty, G.B.; Kelly, D.W. Thermal buckling and elastic vibration of third-order shear deformable functionally graded beams. Int. J. Mech. Sci. 2011, 53, 734–743. [Google Scholar] [CrossRef]
  23. Czechowski, L.; Kowal-Michalska, K. Static and dynamic buckling of rectangular functionally graded plates subjected to thermal loading. Strength Mater. 2013, 45, 666–673. [Google Scholar] [CrossRef]
  24. Czechowski, L. Analysis of dynamic response of functionally graded plate due to temperature pulse load. Compos. Struct. 2017, 160, 625–634. [Google Scholar] [CrossRef]
  25. Bui, T.Q.; Do, T.V.; Ton, L.H.T.; Doan, D.H.; Tanaka, S.; Pham, D.T.; Nguyen-Van, T.-A.; Yu, T.; Hirose, S. On the high temperature mechanical behaviors analysis of heated functionally graded plates using FEM and a new third-order shear deformation plate theory. Compos. Part B 2016, 92, 218–224. [Google Scholar] [CrossRef]
  26. Praveen, G.N.; Reddy, J.N. Nonlinear transient thermal elastic analysis of functionally graded ceramic-metal plates. Int. J. Solids Struct. 1998, 35, 4457–4476. [Google Scholar] [CrossRef]
  27. Kolakowski, Z.; Królak, M. Modal coupled instabilities of thin-walled composite plate and shell structures. Compos. Struct. 2006, 76, 303–313. [Google Scholar] [CrossRef]
  28. Kumar, R.; Lal, A.; Singh, B.N.; Singh, J. New transverse shear deformation theory for bending analysis of FGM plate under patch load. Compos. Struct. 2019, 208, 91–100. [Google Scholar] [CrossRef]
  29. Trabelsi, S.; Frikha, A.; Zghal, S.; Dammak, F. Thermal post-buckling analysis of functionally graded material structures using a modified FSDT. Int. J. Mech. Sci. 2018, 144, 74–89. [Google Scholar] [CrossRef]
  30. Vafakhah, Z.; Navayi Neya, B. An exact three dimensional solution for bending of thick rectangular FGM plate. Compos. Part B 2019, 156, 72–87. [Google Scholar] [CrossRef]
  31. Nguyen, T.T.; Kim, N.I.; Lee, J. Analysis of thin-walled open-section beams with functionally graded materials. Compos. Struct. 2016, 138, 75–83. [Google Scholar] [CrossRef]
  32. Nguyen, T.-T.; Thang, P.T.; Lee, J. Lateral buckling analysis of thin-walled functionally graded open-section beams. Compos. Struct. 2017, 160, 952–963. [Google Scholar] [CrossRef]
  33. Nguyen, T.-T.; Thang, P.T.; Lee, J. Flexural-torsional stability of thin-walled functionally graded open-section beams. Thin Walled Struct. 2017, 110, 88–96. [Google Scholar] [CrossRef]
  34. Xu, G.; Huang, H.; Chen, B.; Chen, F. Buckling and postbuckling of elastoplastic FGM plates under inplane loads. Compos. Struct. 2017, 176, 225–233. [Google Scholar] [CrossRef]
  35. Qin, Z.; Compton, B.G.; Lewis, J.A.; Buehler, M.J. Structural Optimization of 3D-Printed Synthetic Spider Webs for High Strength. Nat. Commun. 2015, 6, 7038. [Google Scholar] [CrossRef] [Green Version]
  36. Dimas, L.S.; Bratzel, G.H.; Eylon, I.; Buehler, M.J. Tough Composites Inspired by Mineralized Natural Materials: Computation, 3D printing, and Testing. Adv. Funct. Mater. 2013, 23, 4629–4638. [Google Scholar] [CrossRef]
  37. Zhang, X.-Y.; Fang, G.; Leeflanf, S.; Zadpoor, A.A.; Zhou, J. Topological design, permeability and mechanical behavior of additively manufactured functionally graded porous metallic biomaterials. Acta Biomater. 2019, 85, 437–452. [Google Scholar] [CrossRef]
  38. Studart, A.R. Additive manufacturing of biologically-inspired materials. Chem. Soc. Rev. 2016, 45, 359–376. [Google Scholar] [CrossRef]
  39. User’s Guide ANSYS® 18.2. Available online: http://www.pmt.usp.br/academic/martoran/notasmodelosgrad/ANSYS%20Fluent%20Users%20Guide.pdf (accessed on 16 November 2019).
Figure 1. A functionally graded material (FGM) plate with its dimensions and a coordinate system (a) and a schematic view of the material distribution (b).
Figure 1. A functionally graded material (FGM) plate with its dimensions and a coordinate system (a) and a schematic view of the material distribution (b).
Materials 12 04090 g001
Figure 2. Material distribution versus a number of layers.
Figure 2. Material distribution versus a number of layers.
Materials 12 04090 g002
Figure 3. Discrete model with boundary conditions (a) and a view of the plate corner (b).
Figure 3. Discrete model with boundary conditions (a) and a view of the plate corner (b).
Materials 12 04090 g003
Figure 4. Static load vs. the normalized deflection in the middle of the plate for the plate SSSS: (a) a different number of layers (b) different variants.
Figure 4. Static load vs. the normalized deflection in the middle of the plate for the plate SSSS: (a) a different number of layers (b) different variants.
Materials 12 04090 g004
Figure 5. Static load vs. the normalized deflection in the middle of the plate for the plate SCSC: (a) a different number of layers (b) different variants.
Figure 5. Static load vs. the normalized deflection in the middle of the plate for the plate SCSC: (a) a different number of layers (b) different variants.
Materials 12 04090 g005
Figure 6. Static load vs. the normalized deflection (FEM) in the middle of the plate for the plate CCCC: (a) a different number of layers (b) different variants.
Figure 6. Static load vs. the normalized deflection (FEM) in the middle of the plate for the plate CCCC: (a) a different number of layers (b) different variants.
Materials 12 04090 g006
Figure 7. Comparison of the plots (FEM) for Var_1: (a) static load vs. the normalized deflection and (b) the critical load of static load vs. the normalized deflection.
Figure 7. Comparison of the plots (FEM) for Var_1: (a) static load vs. the normalized deflection and (b) the critical load of static load vs. the normalized deflection.
Materials 12 04090 g007
Table 1. Thicknesses for individual variants.
Table 1. Thicknesses for individual variants.
Description of Varianttt
[mm]
tc
[mm]
tFGM
[mm]
Var_120.21.8
Var_220.41.6
Var_320.61.4
Var_420.81.2
Var_5211
Table 2. Material properties of basic constituents.
Table 2. Material properties of basic constituents.
ComponentsYoung’s Modulus [GPa]Poisson’s Ratio [-]
Al700.33
Al2033930.25
Table 3. Boundary conditions assumed in the numerical model.
Table 3. Boundary conditions assumed in the numerical model.
Type of BCEdge 1Edge 2Edge 3Edge 4
SSSSuz = 0
ux = moveable
applied load
Couple degree of freedom for nodes in x-directions
uy,uz = 0ux,uz = 0uz = 0
uy = moveable
Couple degree of freedom for nodes in y-directions
SCSCuz = 0
ux = moveable
applied load
Couple degree of freedom for nodes in x-directions
uy,uz = 0
rotx = 0
ux,uz = 0uz = 0
uy = moveable
rotx = 0
Couple degree of freedom for nodes in y-directions
CCCCuz = 0
ux = moveable
roty = 0
applied load
Couple degree of freedom for nodes in x-directions
uy,uz = 0
rotx = 0
ux,uz = 0
roty = 0
uz = 0
uy = moveable
rotx = 0
Couple degree of freedom for nodes in y-directions
Table 4. Critical forces versus a number of layers or boundary conditions.
Table 4. Critical forces versus a number of layers or boundary conditions.
Number of Layers in the FGM (Var_1)Type of BCFEM [N]SAM [N]
5SSSS28,99427,876
7SSSS29,17928,776
11SSSS29,66329,392
15SSSS29,83029,672
5SCSC53,49150,888
7SCSC53,89552,820
11SCSC54,97054,160
15SCSC55,33954,772
5CCCC68,584------
7CCCC69,152------
11CCCC70,667------
15CCCC71,186------
Table 5. Critical forces versus variants of distributions or boundary conditions.
Table 5. Critical forces versus variants of distributions or boundary conditions.
VariantType of BCFEM [N]SAM [N]
Var_1SSSS29,66329,392
Var_2SSSS30,77730,524
Var_3SSSS31,91031,716
Var_4SSSS33,24433,140
Var_5SSSS34,95434,956
Var_1SCSC54,97054,160
Var_2SCSC57,16156,384
Var_3SCSC59,51958,836
Var_4SCSC62,37461,840
Var_5SCSC66,03065,688
Var_1CCCC70,667------
Var_2CCCC73,571------
Var_3CCCC76,776------
Var_4CCCC80,704------
Var_5CCCC85,734------
Table 6. Maps of displacements in direction of Z-axis (Var_1) for the boundary conditions under consideration.
Table 6. Maps of displacements in direction of Z-axis (Var_1) for the boundary conditions under consideration.
SSSSSCSCCCC
Fst = 16.5 kN Materials 12 04090 i001Fst = 16.5 kN Materials 12 04090 i002Fst = 16.5 kN Materials 12 04090 i003
Fst = 88.5 kN Materials 12 04090 i004Fst = 88.5 kN Materials 12 04090 i005Fst = 88.5 kN Materials 12 04090 i006
Fst = 154.5 kN Materials 12 04090 i007Fst = 154.5 kN Materials 12 04090 i008Fst = 154.5 kN Materials 12 04090 i009
Fst = 190.5 kN Materials 12 04090 i010Fst = 183 kN Materials 12 04090 i011Fst = 190.5 kN Materials 12 04090 i012
Fst = 300 kN Materials 12 04090 i013Fst = 300 kN Materials 12 04090 i014Fst = 300 kN Materials 12 04090 i015

Share and Cite

MDPI and ACS Style

Czechowski, L.; Kołakowski, Z. Analysis of the Functionally Step-Variable Graded Plate Under In-Plane Compression. Materials 2019, 12, 4090. https://doi.org/10.3390/ma12244090

AMA Style

Czechowski L, Kołakowski Z. Analysis of the Functionally Step-Variable Graded Plate Under In-Plane Compression. Materials. 2019; 12(24):4090. https://doi.org/10.3390/ma12244090

Chicago/Turabian Style

Czechowski, Leszek, and Zbigniew Kołakowski. 2019. "Analysis of the Functionally Step-Variable Graded Plate Under In-Plane Compression" Materials 12, no. 24: 4090. https://doi.org/10.3390/ma12244090

APA Style

Czechowski, L., & Kołakowski, Z. (2019). Analysis of the Functionally Step-Variable Graded Plate Under In-Plane Compression. Materials, 12(24), 4090. https://doi.org/10.3390/ma12244090

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