Next Article in Journal
Development and Vibration Control of Frequency Adjustable Tuned Mass Damper Based on Magnetorheological Elastomer
Next Article in Special Issue
Influence of Asymmetric Rolling Process on the Microstructure Properties of Bimetallic Sheet Metals
Previous Article in Journal
Shell-to-Beam Numerical Homogenization of 3D Thin-Walled Perforated Beams
Previous Article in Special Issue
Microstructural Evolution and Mechanical Properties of Pure Aluminum upon Multi-Pass Caliber Rolling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Element and Sensitivity Analysis of Anisotropic Thermoelastic Metal and Alloy Discs with Holes

by
Mohamed Abdelsabour Fahmy
1,2,* and
Mohammed Owaidh Alsulami
3
1
Department of Mathematics, Jamoum University College, Umm Al-Qura University, Makkah 25371, Saudi Arabia
2
Faculty of Computers and Informatics, Suez Canal University, New Campus, Ismailia 41522, Egypt
3
Department of Mathematical Sciences, Faculty of Applied Sciences, Umm Al-Qura University, Makkah 24381, Saudi Arabia
*
Author to whom correspondence should be addressed.
Materials 2022, 15(5), 1828; https://doi.org/10.3390/ma15051828
Submission received: 10 February 2022 / Revised: 24 February 2022 / Accepted: 24 February 2022 / Published: 28 February 2022
(This article belongs to the Special Issue Advanced Processing Methods for Metals and Their Alloys)

Abstract

:
The main aim of this paper was to develop an advanced processing method for analyzing of anisotropic thermoelastic metal and alloy discs with holes. In the boundary element method (BEM), the heat impact is expressed as an additional volume integral in the corresponding boundary integral equation. Any attempt to integrate it directly will necessitate domain discretization, which will eliminate the BEM’s most distinguishing feature of boundary discretization. This additional volume integral can be transformed into the boundary by using branch-cut redefinitions to avoid the use of additional line integrals. The numerical results obtained are presented graphically to show the effects of the transient and steady-state heat conduction on the quasi-static thermal stresses of isotropic, orthotropic, and anisotropic metal and alloy discs with holes. The validity of the proposed technique is examined for one-dimensional sensitivity, and excellent agreement with finite element method and experimental results is obtained.

1. Introduction

Thermoelastic analysis is a critical topic in engineering that has sparked a lot of attention in recent years. Thermoelastic research can be carried out using experimental, analytical, and numerical solutions. Only problems with simple geometry and specified boundaries can be solved analytically. To solve problems with complicated boundaries, numerical methods such as the finite element method (FEM) or the boundary element method (BEM) must be utilized.
When thermal effects are considered, several methods for solving the volume integral equation in the boundary element formulation have been presented over the years [1,2,3]. These methods include the dual reciprocity method [4] and multiple reciprocity method [5], particular integral boundary element method [6], and the exact boundary integral transformation method (EBITM) [7]. The EBITM is the most attractive of these boundary element methods since it maintains the BEM’s notion of border discretization without any additional internal treatments or numerical approximations. The EBITM has been successfully employed to transform the volume integral to the boundary in isotropic thermoelasticity [8] and anisotropic thermoelasticity [9].
The main purpose of the considered boundary element analysis is to convert the additional volume integral into the boundary by using branch-cut redefinitions to avoid utilizing additional line integrals.

2. Formulation of the Problem

The thermoelastic field is estimated by solving the elastostatic boundary integral equation with associated thermal data on a region Ω bounded by a surface S , generated first by the boundary element analysis of the following heat conduction equations:
k i j x 1 , x 2 Θ , i j = c ρ Θ ˙ ,   k i j = k j i   ,   k 12 2 k 11 k 22 < 0 ,   i , j = 1 , 2 ,
where k i j , c , ρ , and Θ are anisotropic thermal conductivity coefficients, specific heat capacity, density, and temperature, respectively.
The plane-stress constitutive equations, which describe the relationship between stress σ i j and strain ε i j for a homogeneous, anisotropic solid in the x 1 x 2 plane ( σ 13 = σ 23 = σ 33 = 0 ), can be expressed as follows:
σ 11 σ 22 σ 12 = c 11 c 12 c 16 c 12 c 22 c 26 c 16 c 26 c 66 ε 11 ε 22 2 ε 12 , ε 11 ε 22 2 ε 12 = a 11 a 12 a 16 a 12 a 22 a 26 a 16 a 26 a 66 σ 11 σ 22 σ 12 ,
where c m n and a m n are the stiffness and compliance constants, respectively.
The generalized variable that describes the field point Q x 1 , x 2 position can be expressed as
z j = x 1 x p 1 + μ j x 2 x p 2 ,
where μ j is the material complex constant, and x p 1 ,   x p 2 are the global coordinates of the source point P .

3. Boundary Element Implementation

In the presence of the thermal effect of Equation (1) in Equation (2), we can write the following boundary integral equation [10]:
Θ P + S Θ * q q * Θ d S = q = 1 N Θ q P + S Θ * q q q *   Θ q d S α ¯ q ,
where α ¯ q are unknown coefficients, n is the outward unit normal vector, q is the heat flux vector, Θ q   and   q q are particular solutions, and Θ *   and   q * are anisotropic fundamental solutions.
Equation (4) can be expressed as the following linear algebraic system [11]:
H Θ = G Θ n ,
where H and G are nonsymmetric and symmetric matrices, respectively.
The thermoelastic field is estimated in a progressively coupled technique by solving the elastostatic BIE with accompanying thermal data, generated first by the boundary element solution of heat conduction.
The boundary integral equation based on the thermal effect can now be written as
c i j P u i P = S U i j * P Q t i Q d S S T i j * P . Q u i Q d S + S γ i k n k Θ Q U i j * P , Q d S Ω γ i k Θ , k Q U i j * P , q d Ω ,
where c i j , u i , t i , γ i k , U i j * , and T i j * are the free term, displacements, tractions, thermal moduli, displacement fundamental solutions, and traction fundamental solutions, respectively.
The volume integral is represented by
V j = Ω γ i k Θ , k Q U i j * P , q d Ω
According to Lekhnitskii [12], the fundamental solutions are
U i j * P , Q q = 2 Re β i k A j k ln z k ,
T 1 j * P , Q = 2 n 1 Re μ k 2 A j k z k 2 n 2 Re μ k A j k z k
T 2 j * P , Q = 2 n 1 R e μ k A j k z k + 2 n 2 R e A j k z k .
where β i k and A j k are complex constants, and Re is the real value of the variable between the brackets.
The volume integral in Equation (7) should be redefined in the transformed domain Ω ˆ as follows:
V j = Ω γ i k Θ , k Q U i j * P , q d Ω ˆ
in which
style   γ i k = γ 11 γ 11 K 12 + γ 12 K 11 Δ γ 21 γ 21 K 12 + γ 22 K 11 Δ ,
where Δ = k 11 k 22 k 12 > 0
Then, the volume integral in Equation (11) can be written as follows [13]:
V j = s γ i k Q i j k , t * Θ γ i k Q i j k * Θ , t n t γ i k U i j * Θ n k d S ˆ
where
Q i j k * = 2 Re β i m A j m μ k m μ 1 m 2 + μ 2 m 2 z m ln z m ,
Q i j k , t * = 2 Re β i m A j m μ k m μ t m μ 1 2 + μ 2 2 1 + ln z m ,
and
z m = μ j m x ˆ j x ˆ p j ,
μ m n = K 11 + μ 1 K 12 Δ K 11 + μ 2 K 12 Δ μ 1 μ 2 .
The above transformation can be valid upon adding the line integral as follows:
V j = s [ ( γ i k Q i j k , t * Θ γ i k Q i j k * Θ , t ) n t γ i k U i j * Θ n k ] d S ^ + n = 1 m L 2 n 1 L 2 n 2 L j ( ζ 1 ) d ζ 1 ,
where ζ 1 = x 1 x p 1 is the local coordinate of P , and L 1 , L 2 ,   L 3 , L 4 , ߪ ,   L 2 n 1 , L 2 n 2 are integral intervals of all source points where branch-cut lines intersect the domain.
Furthermore, L j is expressed as
L j = 4 π Θ γ i k K 12 K 11 Im β i p A j p μ 1 p μ k p μ 1 p 2 + μ 2 p 2 + Δ K 11 Im β i p A j p μ 2 p μ k p μ 1 p 2 + μ 2 p 2   4 π γ i k K 12 K 11 Θ , 1 + Δ K 11 Θ , 2 ζ 1 Im β i p A j p μ k p μ 1 p 2 + μ 2 p 2   4 π Θ K 12 K 11 γ i 1 + Δ K 11 γ i 2 Im β i p A j p
where Im is the imaginary value of the variable between the brackets.
All source points whose branch-cut lines cross the domain require the additional line integral to be integrated with intervals L 1 , L 2 ,   L 3 , L 4 , ,   L 2 n 1 , L 2 n 2 . To do so, a strong and robust code is needed to determine all possible intersection points of all source points. Consequently, not only is this inefficient computationally, but it is also extremely difficult to build a strong and robust code to account for all conceivable occurrences of a complex shape. Another significant disadvantage is the requirement to compute internal thermal data along the branch-cut lines involved for evaluation of the extra line integral. However, these internal thermal data can only be obtained after solving the boundary integral equation for thermal analysis, and positions of the internal points are unknown in advance. To do so, the BEM code for computing the associated thermal field must be included into the mechanical analysis code, allowing the computations to be interconnected. This necessitates validating the exact transformation without using the extra line integrals. This procedure is covered in detail in below.
The additional line integral must be integrated over all source points with branch-cut lines that cross the domain.
According to Shiah and Wang [13], we can write
Q j j k , t t * = U i j , k * ,
where
Ω γ i k U i j , k * Θ d Ω ˆ = s γ i k Q i j k , t * Θ γ i k Q i j k * Θ , t n t d S ˆ .  
The derivative of U i j * in Equation (8) yields
U i j , k * = 2 Re β i m A j m μ k m z m .
The generalized variable z m in the polar coordinates can be expressed as
z m = r ˆ μ 11 cos θ ˆ + μ 21 sin θ ˆ .
By substituting Equation (23) into Equation (22), we have
U i j , k = 2 r ˆ Re β i 1 A j 1 μ k 1 μ 11 cos θ ˆ + μ 21 sin θ ˆ + β i 2 A j 2 μ k 2 μ 12 cos θ ˆ + μ 22 sin θ ˆ .
Then, we can re-express U i j , k * as
U i j , k * = 2 r ˆ n = b b D i j k ( n ) e in θ ˆ ,
where b is an integer big enough to ensure series convergence, r ˆ , θ ˆ are the polar coordinates, and D i j k n is the material complex constant. In the current work we considered b = 16 .
From Equation (23) and using the theory of Fourier series, we can write D i j k n as
D n jjk = 1 2 π π π Re β i 1 A j 1 μ k 1 μ 11 cos θ ˆ + μ 21 sin θ ˆ + β i 2 A j 2 μ 2 μ 12 cos θ ˆ + μ 22 sin θ ˆ e in θ ˆ d θ ˆ .
Then, Q i j k * which satisfies Equation (20) can be expressed as
Q i j k * r ˆ , θ ˆ = r ˆ ρ i j k θ ˆ .
Substituting Equation (23) into Equation (20), we get
d 2 ρ i j k ( θ ˆ ) d θ ˆ 2 + ρ i j k ( θ ˆ ) = 2 n = b b D i j k ( n ) e in θ ˆ .
From Equation (28), ρ i j k θ ˆ is determined to be
ρ i j k ( θ ˆ ) = n = b b 2 ( 1 n 2 ) D i j k ( n ) e in θ ˆ ,
when n = ± 1 , ρ i j k θ ˆ is not properly defined in Equation (29). As a result, the Fourier series of U i j , k * is split into two halves—one for n   ± 1 and one for n = ± 1 .
Thus, we can write U i j , k * as
U i j , k = 2 r ˆ D i j k 1 e i θ ˆ + D i j k 1 e i θ ˆ ,
where Q i j k * can be expressed as
Q i j k * r ˆ , θ ˆ = r ˆ ln r ˆ λ i j k θ ˆ .
Substituting Equation (31) into Equation (20), we obtain
λ i j k θ ˆ + ln r ˆ d 2 λ i j k θ ˆ d θ ˆ 2 + λ i j k θ ˆ = 2 D i j k 1 e i θ ˆ + D i j k 1 e i θ ˆ ,
where Equation (32) can be satisfied under the following condition:
λ i j k θ ^ = D i j k 1 e i θ ^ + D i j k 1 e i θ ^ .
Then, by using Equations (27), (29), (31), and (33), we can write Q i j k * as
Q i j k * ( r ^ , θ ^ ) = r ^ n = b n ± 1 b 2 D i j k ( n ) ( 1 n 2 ) e i n θ ^ + r ^ ln r ^ ( D i j k ( 1 ) e i θ ^ + D i j k ( 1 ) e i θ ^ ) ,
where D i j k n is computed using Equation (26).
The explicit expression of Q i j k , t * may be written as follows [13]:
Q i j k , 1 * = n = b n ± 1 b 2 D j j k ( n ) ( 1 n 2 ) ( cos θ insin θ ) e in θ ^ + ln r ^ ( D i j k ( 1 ) + D i j k ( 1 ) ) + cos θ ^ ( D i j k ( 1 ) e i θ ^ + D i j k ( 1 ) e i θ ^ ) ,  
Q i j k , 2 * = n = b n ± 1 b 2 D i j k ( n ) ( 1 n 2 ) ( sin θ ^ + in   cos θ ^ ) e in θ ^ + ln r ^ ( D i j k ( 1 ) D i j k ( 1 ) ) + sin θ ( D i j k ( 1 ) e i θ ^ + D i j k ( 1 ) e i θ ^ ) .

4. Numerical Results and Discussion

To validate the analysis, we consider a thick hollow disc with four inside holes [13]. The current problem was analyzed using FlexPDE 7 which is based on the finite element method (FEM). For the elastic boundary conditions, the exterior surface was fully constrained, although the interior was traction-free in all directions. Moreover, insulation was required on all surfaces of the four inside holes.
The thermoelastic constants of monoclinic graphite–epoxy can be written as follows:
Elasticity tensor
C p j k l = 430.1 130.4 18.2 130.4 116.7 21.0 18.2 21.0 73.6 0   0   201.3 0   0   70.1 0   0   2.4 0 0 0 0 0 0 201.3 70.1   2.4 19.8 8.0     0   8.0 29.1 0     0 0   147.3   GPa ;
Mechanical temperature coefficient
β p j = 1.01 2.00 0 2.00 1.48 0 0 0 7.52 · 10 6   N / km 2 ;
Thermal conductivity tensor
k p j = 5.2 0 0 0 7.6 0 0 0 38.3   W / km ;
Mass density ρ = 7820   kg / m 3 ;
Heat capacity c = 461   J / kg   K .
Figure 1, Figure 2 and Figure 3 show the effect of transient and steady-state heat conduction on the quasi-static thermal stresses σ 11 , σ 12 , and σ 22 along the x -axis in the transient, as well as the steady-state heat conduction for isotropic, orthotropic, and anisotropic cases.
Figure 1 depicts the distribution of the quasi-static thermal stress component σ 11 in the context of the isotropic, orthotropic, and anisotropic materials. It demonstrates that, in the context of the three considered materials, σ 11 increases at first to a maximum value, before decreasing to a minimum value in the transient heat conduction (THC) and steady-state heat conduction (SHC). It shows that, in the context of the three considered materials, σ 11 converges to zero with increasing distance x.
Figure 2 depicts the distribution of the quasi-static thermal stress component σ 12 , demonstrating that the stress component σ 12 reaches a zero value. It shows that, in the context of isotropic and orthotropic materials, σ 12 increases to a maximum value, before sharply decreasing in the transient and steady-state heat conduction. However, in the context of the anisotropic material,   σ 12 increases to a maximum value, before decreasing to a minimum value in the transient and steady-state heat conduction. In the context of the three materials, σ 12 converges to zero with increasing distance x , while the values of the stress component σ 12 in the transient heat conduction are higher than those in the steady-state heat conduction.
Figure 3 shows the distribution of the quasi-static thermal stress component σ 22 in the context of the three materials; it begins with a negative value decrease in the transient and steady-state heat conduction. The values of σ 22 in the context of the isotropic and orthotropic materials increase initially in the range 0 x 0.7 , before decreasing in the range 0.7 x 1 , and again increasing to a maximum value in the transient and steady-state heat conduction. However, in the context of the anisotropic material, σ 22 increases at first, before decreasing to a minimum value in the transient and steady-state heat conduction The values of the stress component σ 22 in the transient heat conduction are higher than those in the steady-state heat conduction and converge to zero with increasing distance x in the transient and steady-state heat conduction.
Quasi-static thermal stresses σ 11 and σ 12 of anisotropic thermoelastic metal and alloy discs with holes were calculated on the inner surface of the disc, as well as the other four holes inside. Figure 4 and Figure 5 show the sensitivity variations of the thermal stresses σ 11 and σ 12 along the x -axis in the transient and steady-state heat conduction for isotropic, orthotropic, and anisotropic cases.
The numerical results of the quasi-static thermal stresses σ 11 and σ 12 for the four holes are presented in Figure 6 and Figure 7 for BEM and FEM analysis of anisotropic thermoelastic metal and alloy discs with holes. It can be seen from these figures that the BEM results are in very good agreement with the FEM results.
The advantages of spark plasma sintering (SPS) over conventional hot pressing (HP) or hot isostatic pressing (HIP), such as reduced sintering time and temperatures, minimize grain growth and frequently result in improved mechanical, physical, or optical properties. Therefore, SPS has been used successfully to manufacture a variety of metals. The thermoelastic problem studied in Wang et al. [14] can be treated as a special case of our study of analysis of anisotropic thermoelastic metal and alloy discs with holes. To analyze the thermal stress sensitivity distribution in the SPS process, we used the same model as in [14] with the boundary element method (BEM). In this special case under consideration, the numerical results for the quasi-static thermal stresses σ 11 and σ 12 are shown in Figure 8 and Figure 9. It can be seen from these figures that our BEM results are in excellent agreement with the FEM results and experimental (Exp) results of Wang et al. [14]. We refer the interested readers to [15,16,17,18,19,20,21,22] for further references.
Table 1 shows a comparison of the computer resources needed in the analysis of anisotropic thermoelastic metal and alloy discs, for BEM with additional line integrals (Case 1) versus BEM without additional line integrals (Case 2). It can be seen from this table that the proposed BEM without additional line integrals is more accurate and efficient than the BEM with additional line integrals.

5. Conclusions

The examination of the numerical results and figures enables us to make some concluding remarks:
  • The current research has received a lot of attention because of its practical applications in fields such as astronautics, geomechanics, earthquake engineering, nuclear reactors, material science, and other industrial applications.
  • Because the proposed boundary element approach only needs to solve the boundary unknowns, it solves problems faster and more accurately than domain approaches while also minimizing the solver’s processing costs.
  • Avoiding the use of additional line integrals by using branch-cut redefinitions in the current study plays a significant role in all the physical quantities and their design sensitivities.
  • The current results were validated against the numerical and experimental results obtained through other methods previously validated. It should be noted that the BEM results are in excellent agreement with the FEM and experimental results, confirming the accuracy of the BEM technique.
  • Current numerical results for our complex and general problem may be of interest to engineers and material science researchers, as well as those working on the development of anisotropic thermoelastic metal and alloy discs with holes.
  • It can be concluded from analysis results that the proposed technique is more efficient than other techniques in the literature for analyzing anisotropic thermoelastic metal and alloy discs with holes.
  • The numerical results show that the proposed BEM is ideal for analyzing anisotropic thermoelastic metal and alloy discs with holes.

Author Contributions

Conceptualization, M.A.F. and M.O.A.; methodology, M.A.F.; software, M.A.F.; validation, M.A.F. and M.O.A.; formal analysis, M.A.F.; investigation, M.A.F.; resources, M.A.F.; data curation, M.A.F.; writing—original draft preparation, M.A.F. and M.O.A.; writing—review and editing, M.A.F.; visualization, M.A.F.; supervision, M.A.F.; project administration, M.A.F.; funding acquisition, M.A.F. and M.O.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Deanship of Scientific Research at Umm Al-Qura University] grant number [22UQU4340548DSR01] And the APC was funded by [Deanship of Scientific Research at Umm Al-Qura University].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analysed during this study are included in this published article.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: (22UQU4340548DSR01).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sládek, V.; Sládek, J. Boundary integral equation method in two-dimensional thermoelasticity. Eng. Anal. 1984, 1, 135–148. [Google Scholar] [CrossRef]
  2. Sládek, V.; Sládek, J. Boundary integral equation method in thermoelasticity part III:—Uncoupled thermoelasticity. Appl. Math. Model. 1984, 8, 413–418. [Google Scholar] [CrossRef]
  3. Sládek, J.; Sládek, V.; Markechová, I. Boundary element method analysis of stationary thermoelasticity problems in non-homogeneous media. Int. J. Numer. Methods Eng. 1990, 30, 505–516. [Google Scholar] [CrossRef]
  4. Nardini, D.; Brebbia, C.A. A new approach to free vibration analysis using boundary elements. Appl. Math. Model. 1983, 7, 157–162. [Google Scholar] [CrossRef]
  5. Nowak, A.J.; Brebbia, C.A. The multiple-reciprocity method. A new approach for transforming BEM domain integrals to the boundary. Eng. Anal. Boundary Elem. 1989, 6, 164–167. [Google Scholar] [CrossRef]
  6. Deb, A.; Banerjee, P.K. BEM for general anisotropic 2D elasticity using particular integrals. Commun. Appl. Numer. Methods 1990, 6, 111–119. [Google Scholar] [CrossRef]
  7. Rizzo, F.J.; Shippy, D.J. An advanced boundary integral equation method for three-dimensional thermoelasticity. Int. J. Numer. Methods Eng. 1977, 11, 1753–1768. [Google Scholar] [CrossRef]
  8. Danson, D. Linear isotropic elasticity with body forces. In Progress in Boundary Element Methods; Brrebbia, C.A., Ed.; Springer: New York, NY, USA, 1983; Volume 2, pp. 101–135. [Google Scholar]
  9. Shiah, Y.C.; Tan, C.L. Exact boundary integral transformation of the thermoelastic domain integral in BEM for general 2D anisotropic elasticity. Comput. Mech. 1999, 23, 87–96. [Google Scholar] [CrossRef]
  10. Fahmy, M.A. A time-stepping DRBEM for the transient magneto-thermo-visco-elastic stresses in a rotating non-homogeneous anisotropic solid. Eng. Anal. Bound. Elem. 2012, 36, 335–345. [Google Scholar] [CrossRef]
  11. Fahmy, M.A. Transient Magneto-Thermo-Elastic Stresses in an Anisotropic Viscoelastic Solid with and Without a Moving Heat Source. Numer. Heat Transf. Part A Appl. 2012, 61, 633–650. [Google Scholar] [CrossRef]
  12. Lekhnitskii, S.G. Anisotropic Plates; Gordon & Breach Science Publisher: New York, NY, USA, 1968. [Google Scholar]
  13. Shiah, Y.C.; Wang, S.-H. New Domain Integral Transformation in Boundary Element Analysis for 2D Anisotropic Thermoelasticity. J. Eng. Mech. 2016, 142, 04016065. [Google Scholar] [CrossRef]
  14. Wang, C.; Cheng, L.; Zhao, Z. FEM analysis of the temperature and stress distribution in spark plasma sintering: Modelling and experimental validation. Comput. Mater. Sci. 2010, 49, 351–362. [Google Scholar] [CrossRef]
  15. Fahmy, M.A. A novel BEM for modeling and simulation of 3T nonlinear generalized anisotropic micropolar-thermoelasticity theory with memory dependent derivative. CMES-Comput. Model. Eng. Sci. 2021, 126, 175–199. [Google Scholar]
  16. Fahmy, M.A. A new boundary element formulation for modeling and simulation of three-temperature distributions in carbon nanotube fiber reinforced composites with inclusions. Math. Methods Appl. Sci. 2021. [Google Scholar] [CrossRef]
  17. Fahmy, M.A. Boundary element modeling of 3 T nonlinear transient magneto-thermoviscoelastic wave propagation problems in anisotropic circular cylindrical shells. Compos. Struct. 2021, 277, 114655. [Google Scholar] [CrossRef]
  18. Fahmy, M.A. A new boundary element algorithm for a general solution of nonlinear space-time fractional dual-phase-lag bio-heat transfer problems during electromagnetic radiation. Case Stud. Therm. Eng. 2021, 25, 100918. [Google Scholar] [CrossRef]
  19. Fahmy, M.A. A new boundary element algorithm for modeling and simulation of nonlinear thermal stresses in micropolar FGA composites with temperature-dependent properties. Adv. Model. Simul. Eng. Sci. 2021, 8, 6. [Google Scholar] [CrossRef]
  20. Fahmy, M.A. A New BEM for Fractional Nonlinear Generalized Porothermoelastic Wave Propagation Problems. Comput. Mater. Contin. 2021, 68, 59–76. [Google Scholar] [CrossRef]
  21. Fahmy, M.A. A New BEM Modeling Algorithm for Size-Dependent Thermopiezoelectric Problems in Smart Nanostructures. Comput. Mater. Contin. 2021, 69, 931–944. [Google Scholar] [CrossRef]
  22. Fahmy, M.A. Boundary element algorithm for nonlinear modeling and simulation of three temperature anisotropic generalized micropolar piezothermoelasticity with memory-dependent derivative. Int. J. Appl. Mech. 2020, 12, 2050027. [Google Scholar] [CrossRef]
Figure 1. Variation of the thermal stress σ11 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Figure 1. Variation of the thermal stress σ11 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Materials 15 01828 g001
Figure 2. Variation of the thermal stress σ12 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Figure 2. Variation of the thermal stress σ12 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Materials 15 01828 g002
Figure 3. Variation of the thermal stress σ22 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Figure 3. Variation of the thermal stress σ22 along x-axis in the transient and steady-state heat con-duction for isotropic, orthotropic, and anisotropic cases.
Materials 15 01828 g003
Figure 4. Variation of the thermal stress σ11 sensitivity along x-axis of the disc inner surface for BEM and FEM.
Figure 4. Variation of the thermal stress σ11 sensitivity along x-axis of the disc inner surface for BEM and FEM.
Materials 15 01828 g004
Figure 5. Variation of the thermal stress σ12 sensitivity along x-axis of the disc inner surface for BEM and FEM.
Figure 5. Variation of the thermal stress σ12 sensitivity along x-axis of the disc inner surface for BEM and FEM.
Materials 15 01828 g005
Figure 6. Variation of the thermal stress σ11 sensitivity along x-axis on the holes inside for BEM and FEM.
Figure 6. Variation of the thermal stress σ11 sensitivity along x-axis on the holes inside for BEM and FEM.
Materials 15 01828 g006
Figure 7. Variation of the thermal stress σ12 sensitivity along x-axis on the holes inside for BEM and FEM.
Figure 7. Variation of the thermal stress σ12 sensitivity along x-axis on the holes inside for BEM and FEM.
Materials 15 01828 g007
Figure 8. Variation of the thermal stress σ11 sensitivity along x-axis on the holes inside for BEM, FEM, and Exp.
Figure 8. Variation of the thermal stress σ11 sensitivity along x-axis on the holes inside for BEM, FEM, and Exp.
Materials 15 01828 g008
Figure 9. Variation of the thermal stress σ12 sensitivity along x-axis on the holes inside for BEM, FEM, and Exp.
Figure 9. Variation of the thermal stress σ12 sensitivity along x-axis on the holes inside for BEM, FEM, and Exp.
Materials 15 01828 g009
Table 1. Comparison of computer resources required for BEM with additional line integrals (Case 1) and BEM without additional line integrals (Case 2).
Table 1. Comparison of computer resources required for BEM with additional line integrals (Case 1) and BEM without additional line integrals (Case 2).
FEMBEM (Case 1)BEM (Case 2)
CPU time (min)28244
Memory (MB)26221
Disc space (MB)38320
Accuracy of results (%)2.22.11.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fahmy, M.A.; Alsulami, M.O. Boundary Element and Sensitivity Analysis of Anisotropic Thermoelastic Metal and Alloy Discs with Holes. Materials 2022, 15, 1828. https://doi.org/10.3390/ma15051828

AMA Style

Fahmy MA, Alsulami MO. Boundary Element and Sensitivity Analysis of Anisotropic Thermoelastic Metal and Alloy Discs with Holes. Materials. 2022; 15(5):1828. https://doi.org/10.3390/ma15051828

Chicago/Turabian Style

Fahmy, Mohamed Abdelsabour, and Mohammed Owaidh Alsulami. 2022. "Boundary Element and Sensitivity Analysis of Anisotropic Thermoelastic Metal and Alloy Discs with Holes" Materials 15, no. 5: 1828. https://doi.org/10.3390/ma15051828

APA Style

Fahmy, M. A., & Alsulami, M. O. (2022). Boundary Element and Sensitivity Analysis of Anisotropic Thermoelastic Metal and Alloy Discs with Holes. Materials, 15(5), 1828. https://doi.org/10.3390/ma15051828

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