Published: 21 November 2020

Numerical analysis of the damage mechanics variable and vibration of a viscothermoelastic microbeam with variable thermal conductivity

Hamdy M. Youssef1
Alaa A. El-Bary2
Haitham M. Atef3
A. H. El-Sharif4
1Mathematics Department, Faculty of Education, Alexandria University, Alexandria, Egypt
1Mechanical Engineering Department, College of Engineering and Islamic Architecture, Umm Al-Qura University, Makkah, Saudi Arabia
2Basic and Applied Science Institute, Arab Academy for Science, Technology and Maritime Transport, P.O. Box 1029, Alexandria, Egypt
3Faculty of Science, Department of Mathematics, Damanhur University, Damanhur, Egypt
4Mathematics Department, Faculty of Science, Benghazi University, Benghazi, Libya
Corresponding Author:
Hamdy M. Youssef
Views 149
Reads 52
Downloads 1062

Abstract

In the present paper, the analysis for thermoelastic homogeneous isotropic microbeams has been constructed. A generalized viscothermoelasticity theory of one relaxation time with variable thermal conductivity in the context of damage mechanics definition has been applied based on simply supported conditions for aspect ratios. Laplace transform has been applied for the governing differential equations, and its inverse has been carried out by using the Tzou method. Microbeam of silicon nitride has been considered when it is subjected to ramp-type heating and simply supported. The results have been illustrated in figures to stand on the impacts of the viscothermoelastic parameter, the thermal conductivity parameter, the value of the beam thickness, and the ramp-type heating parameter. The influences of the mentioned parameters are significant on all the studied functions, and the ramp-type heating parameter plays a vital role in thermodynamically damping of the energy which has been generated on the beam.

Numerical analysis of the damage mechanics variable and vibration of a viscothermoelastic microbeam with variable thermal conductivity

Highlights

  • An analysis for thermoelastic homogeneous isotropic microbeams has been constructed.
  • A generalized viscothermoelasticity theory of one relaxation time with variable thermal conductivity in the context of damage mechanics definition has been applied.
  • Microbeam of silicon nitride has been considered when it is subjected to ramp-type heating and simply supported.
  • The thermal conductivity parameter, the value of the beam thickness, and the ramp-type heating parameter have significant effects on all the studied functions.
  • The ramp-type heating parameter plays a vital role in thermodynamically damping of the energy which has been generated on the beam.

1. Introduction

Tzou is the first one who studied the heat conduction by solving mathematical models based on dual-phase lag (DPL) [1, 2]. Many authors and researchers applied that model in the thermal transfer applications [3-8]. The model of coupled thermoelasticity theory is one of the first types of heat conduction, which consists the equation of motion and the equation of energy conservation using the classical Fourier’s law of heat conduction [9-12]. Lord and Shulman (L-S) modified the classical Fourier’s law by inserting the relaxation time (lag time) for an isotropic body [13]. Within this model, the heat conduction law has been modified to include both heat flux and its time derivative which is called non-Fourier’s or Cattaneo’s law of heat conduction. In the context of this model, the heat equation law is a hyperbolic type of differential equation and then eliminates the paradox which coming from infinite speeds of propagation of the thermal wave [14-18].

The vibration of microbeam resonators is the most important of the micro-/nano-beam subject. Alghamdi studied the damping vibration of beam resonator with voids by dual-phase-lag generalized thermoelasticity model [9]. Sharma and Grover discussed the transverse vibrations of a homogenous isotropic and thermoelastic microscale and nanoscale beam resonators with voids [19]. Sun and Saka induced the thermal damping vibration for thermoelastic of micro-circular plate resonators [20]. They inserted a new factor in their formula of thermal damping based on Poisson's ratio, which is not the same of Lifshitz and Roukes [21]. Al-Lehaibi and Youssef studied the vibration of gold nano-beam due to thermal shock [22]. Kidawa studied the effects of internal and external damping on transverse vibrations of a microbeam due to moving heat source by using the properties of the Green functions [23]. Boley discussed the thermal vibrations of a simply supported rectangular nanobeam due to a thermal shock distributed through its span [24]. Manolis and Beskos studied the heat transfer induced vibration of microbeams structures; he used a numerical method of analysis to the thermal of the elastic dynamic response of beam structure to thermal loading [25]. Al-Huniti et al. studied the displacements and stresses of heated a bar due to a thermal wave generated by a moving laser beam with high-power [26].

Studying of viscoelastic materials or thermal and mechanical relaxations properties became vital in mechanics. The theory of viscothermoelasticity and vibrational principles have been studied in thermodynamics by Biot [27, 28]. Drozdov introduced a constitutive thermal model for the viscothermoelasticity behavior with finite strains of polymers [29]. Ezzat and El-Karmany applied a different model of viscothermoelasticity for isotropic and homogenous medium to study the influence of the thermal and mechanical relaxation of volume properties of viscothermoelastic materials [30]. Carcione et al. applied a computational algorithm for wave simulations in a thermoelastic materials by using the Kelvin-Voigt model [31]. Grover investigated the transverse vibrations in viscothermoelastic microbeams [32-34]. Sharma and Grover discussed the closed-form formulations for the transverse vibrations of a homogenous thermoelastic thin micro/nano-beam with voids [19]. Grover and Seth solved a problem of viscothermoelastic micro-scale beam resonators based on the model of dual-phase-lag (DPL) [35].

Usually, the properties of any material are constant during any process. In engineering ceramics applications, significant variations do, however, occur over the operating temperature range, in the thermal conductivity coefficient particularly. Godfrey found out that up to 45 % of decreases in the value of the thermal conductivity of different samples of silicon nitride in the range of increment (1.0 °C-400 °C). So, what are the impacts of these variations on the temperature increment, lateral deflection, stress, and energy in metal components [36]. Thus, the temperature-dependent material properties must be proposed in the analysis of thermal stress. Youssef, with many co-authors, solved many applications for thermoelastic materials with variable thermal conductivity [37-39]. Zenkour and Abbas solved a model of an infinite annular cylinder based on generalized thermoelasticity with one relaxation time for temperature-dependent properties [40].

Products typically contain many flaws in the original condition, such as microcracks or voids. Such internal holes or voids may expand and converge during a deformation cycle, whereas other content is totally removed by the development of new micro defects in stress concentrators called mechanical damage. At the latter stage, the credibility of the materials has been totally lost and macro-cracks have been created. Methane damage was graded as brittle, crawl, ductile, and fatigue according to the growing macroscopic phenomena [41].

Throughout the continuum mechanics system there is also a definition of the macros actions of the mechanically affected material. In certain cases, the principle of mechanical harm may relate the dynamics of fracture to conventional continuum mechanics. Damage variables can be introduced in different ways.

In a cross-section of the mechanically damaged medium, we thus assume an area of small element dA with the n vector of the unit normal. The area of the defects is dAυ and the quantity of mechanical damage can be measured by the area fraction [41, 42]:

1
υn=dAυdA, 0υn1,

where υn=0 classifies the undamaged material, while υn=1 theoretically comply with the totally damaged with a complete loss of stress hold ability (fracture case). In any natural material, with values of υn0.2,...,0.5 processes taking place, which leads to a total failure. If the mechanical damage is constant through a finite area, under uniaxial tension, the relation Eq. (1) reduces to:

2
υn=AυA.

The effects of microcracks that are inclined to the cross-section of the materials cannot be described correctly in the same way. Correspondingly, in case of isotropic damage υ independent of n, hence, the effective stresses are given by [41, 42]:

3
σij=1-υσ~ij,

where σ~ij are the average stresses in the undamaged material. Many applications and problems have been published under this definition of damage mechanics [42-47].

2. Basic Equations

A homogenous isotropic thermally conducting, Kelvin-Voigt type viscothermoelastic solid in Cartesian coordinate has been considered and initially undeformed at a uniform temperature T0. The displacement vector is defined as Ux,y,z,t=u,v,w, and the absolute temperature is given by Tx,y,z,t. The governing differential equations in the context of generalized thermoelasticity based on the non-Fourier heat conduction law will be constructed in the absence of any external body forces and heat sources as follows [34]:

The equations of motion as [2, 37]:

4
σij,j=ρU¨i.

The heat conduction law in the form [48]:

5
Kφφ,i,i=t+τ02t2ρCυφ+βT01-υδijeij.

The stress-strain constitutive equations are in the form [48]:

6
σij=1-υλδijekk+2μeij-βδijφ.

The deformation-displacement relations are in the form [48]:

7
eij=12Ui,j+Uj,i.

For the viscothermoelastic materials, Lame’s parameters are in the following forms [32, 35, 49]:

8
λ=λ01+λ1 t, μ=μ01+μ1 t, β=3λ+2μαT,

where the indices i,j=x,y,z, ρ is the density, αT is the coefficient of linear thermal expansion, λ1, μ1 are the viscoelastic relaxation times, τ0 is the lag time or thermal relaxation time, λ0, μ0 Lame’s parameter in the usual case, Kφ is the thermal conductivity, and φ=T-T0 is the temperature increment.

The specific heat at constant strain satisfies the following relation [48]:

9
Cυ=Kρ κ,

where κ is the diffusivity.

Consider the following mapping [48]:

10
θ=1K00 φKξdξ,

where K0 is the usual thermal conductivity.

Differentiating Eq. (10) for the coordinates xi, we get:

11
K0θ,i=Kφ φ,i.

Hence, we have:

12
K0θ,ii=Kφ φ,i,i.

Differentiating Eq. (10) for time, we get:

13
K0θ˙=Kφφ˙.

Thus, the heat Eq. (5) takes the form:

14
θ,ii= t+τ02 t2θκ+βT0K01-υδijeij.

3. Problem formulation

We assume small flexural deflections of a thin visco-thermoelastic nanobeam of length l(0xl), with b(-b/2yb/2), and thickness h(-h/2zh/2), for which the x, y, and z axes are defined along the longitudinal, with, and thickness directions of the beam, respectively. In equilibrium, the beam is unstrained, unstressed, without any mechanical damping, bending [6].

In this work, the well-known Euler-Bernoulli equation in one-dimensional form will be used. Thus, the displacements components are given by [32, 50, 51]:

15
Ux,z,t=-z wx,t x ,0 ,wx,t.

The flexural moment of the cross-section is given by [1, 2, 9, 10, 19]:

16
Mx,t=λ+2μ1-υI2wx,tx2+1-υβMTx,t,

where I=bh3/12 is the moment of inertia of the cross-section about the x-axis, and MT is the thermal moment of the beam about the x-axis which is given by [42, 50-52]:

17
MTx,t=b-h/2h/2φx,z,tzdz.

Thus, the equation of the induced lateral deflection may be expressed in the form [32, 49]:

18
1-υλ+2μI4 wx,t x4 +ρA2 wx,t t2+1-υ3λ+2μαT2 MTx,t x2=0,

where wx,t is the lateral deflection, A=hb is the cross-section area.

The non-Fourier heat conduction law in Eq. (14) takes the form [32]:

19
2θx,z,t x2+2θx,z,t z2= t+τ02 t2θx,z,tκ +αTT01-υ3λ+2μK0ex,z,t,

and:

20
ex,z,t=ux,z,tx+vx,z,ty+wx,z,tz,

where ex,z,t is the volumetric strain which gives from Eq. (15) that:

21
ex,z,t=-z2wx,tx2.

From the relation in Eq. (8), we have the following:

22
λ+2μ=λ0+2μ01+β1t, 3λ+2μ=3λ0+2μ01+β2t.

β1=λ0λ1+2μ0μ1λ0+2μ0 and β2=3λ0λ1+2μ0μ13λ0+2μ0 are the aggregation of the viscothermoelastic relaxation times parameters.

We consider that the thermal conductivity is a function on the temperature increment as follows [48]:

23
Kφ=K01+K1φ.

Thus, the mapping in Eq. (10) takes the form:

24
θ=φ+K12φ2,

where K1 is a small constant (K1=0gives the usual thermal conductivity).

Because there is no heat-flow across the upper and lower surfaces of the beam, so that φx,y,zzz=±h/2=0. For a thin beam and assuming the temperature varies in terms of a sin(pz) function through the thickness direction, where p=π/h , gives:

25
φx,z,t=ϑ1x,tsinpz.

From Eq. (24) and (25), we have:

26
θx,z,tzz=±h/2= φx,z,tz+K1φx,z,tφx,z,tzz=±h/2=0.

Hence, we consider that:

27
θx,z,t=ϑ2x,tsinpz.

Then, form Eqs. (17), (18), and (25), we obtain:

28
1-υλ+2μ4 wx,t x4 +12ρh22 wx,t t2 +12αTh31-υ3λ+2μ2ϑ1x,t x2-h/2h/2sinpzzdz=0.

Furthermore, from Eqs. (19) and (27), we get:

29
2ϑ2x,t x2-p2ϑ2x,tsinpz= t+τ02 t2ϑ2x,tκsinpz
-T0αTK01-υ3λ+2μ t+τ02 t2z2wx,tx2.

Thus, Eqs. (8) and (28) gives:

30
1-υλ0+2μ01+β1t4 wx,t x4 +12ρh22 wx,t t2
+12αT1-υ3λ0+2μ0h31+β2t2ϑ1x,t x2-h/2h/2 z sinpzdz=0 .

By carrying out the integrations, the Eq. (30) takes the form:

31
1-υλ0+2μ01+β1t4 wx,t x4 +12ρh22 wx,t t2
+24αT1-υ3λ0+2μ0hπ21+β2t2ϑ1x,t x2=0 .

Multiply both sides of Eq. (29) by z and integrating for z from -h/2 to h/2 , then, we obtain:

32
2ϑ2x,t x2-p2ϑ2x,t= ε t+τ02 t2ϑ2x,t
-T0hπ2αT1-υ3λ0+2μ024K0 t+τ02 t21+β2t2wx,t x2,

where ε=1κ=ρCυK0.

From Eqs. (24), (25), and (27), we have:

33
ϑ2x,z,t= ϑ1x,t+K12ϑ12x,tsinpz,

which gives:

34
2ϑ2x,tx2=2ϑ1x,tx2+K1sinpzϑ1x,t2ϑ1x,tx2+ϑ1x,tx2.

Since:

35
K1sinpzϑ1x,t2ϑ1x,tx2+ϑ1x,tx22ϑ1x,tx2.

Then, we can consider the following approximation:

36
2ϑ2x,tx22ϑ1x,tx2.

Hence, the Eq. (31) takes the form:

37
λ0+2μ01+β1t4 wx,t x4 +12ρh21-υ2 wx,t t2
+24αT3λ0+2μ0hπ21+β2t2ϑ2x,t x2=0 .

For simplicity, we use the non-dimensional variables as follows [12, 37]:

38
x',w',h',l'=ε c0x,w,h,l , t',τ'0,β'1,β'2=ε c02t,τ0,β1,β2,
σ'=σλ0+2μ0, ϑ'2=ϑ2T0,

where c0=λ0+2μ0/ρ.

Hence, we have:

39
1+β1t4 wx,t x4 +ε12 wx,t t2+ε21+β2t2ϑ2x,t x2=0,
40
2ϑ2x,t x2-ε3ϑ2x,t= t+τo2 t2ϑ2x,t-ε41+β2t2wx,t x2,
41
σxxx,z,t=1+β1tex,z,t-γ1+β2tφx,z,t,

where:

ε1=12h21-υ, ε2=24γπ2h , ε3=p2,
ε4=π2h αT1-υ3λ0+2μ024K0ε , γ=T0αT3λ0+2μ0λ0+2μ0.

(We dropped the prime for simplicity).

4. Formulation in the Laplace transform domain

The Laplace transform for Eq. (39) and (40), which is defined by the following formula will be applied:

42
f- s=0fte-stdt,

where the inverse of the Laplace transform takes the form:

43
L-1f-s=fteκtt12f-κ+Ren=1N-1nf-κ+inπt,

where κt4.7 Tzou [2], Re is the real part and i is the imaginary unit, for faster convergence, many numerical applications and experiments have shown that the value of κ satisfies the relation.

Hence, we obtain the following system of differential equations:

44
1+β1s4 w-x,s x4 +ε1s2w-+ε21+β2s2ϑ-2x,s x2=0,
45
2ϑ-2x,s x2-ε3ϑ-2x,s=s+τos2ϑ-2x,s-ε4s+τos21+β2s2w-x,s x2,
46
σ-xxx,z,s=1-υ1+β1se-x,z,s-γ1+β2sφ-x,z,s,
47
e-x,z,s=-z2w-x,sx2.

During applying the Laplace transform, the following initial conditions have been used:

48
ϑ2x,0=wx,0=ϑ2x,0t=wx,0t=0.

We can re-write the above system in Eqs. (44) and (45) to be in the forms:

49
D4 +ε5w-x,s=-ε6D2ϑ-2x,s,
50
D2-ε7ϑ-2x,s=-ε8D2w-x,s,

where:

Dr=r xr, ε5=ε1s21+β1s, ε6=ε21+β2s1+β1s,
ε7=ε3+s+τos2, ε8= ε4s+τos21+β2s.

The system in Eqs. (49) and (50) generates the following characteristic equation:

51
D6-LD4+MD2-Nϑ-2x,s,w-2x,s =0,

where L=ε7 +ε6ε8, M=ε5, N=ε5ε7.

The solutions of the Eqs. (51) take the forms:

52
ϑ-2x,s=-ε8i=13Aiki2sinhkil-x,
53
w-x,s=i=13Aiki2-ε7sinhkil-x,

where ±k1, ±k2, ±k3 are the roots of the characteristic equation:

54
k6 -L k4+M k2-N =0.

To calculate the constants Ai=Ais, i= 1, 2, 3, we must apply any set of boundary conditions, so we consider that the beam is thermally loaded and simply supported as following:

55
w0,t=2w0,tx2=0, φ0,t=ϑ0gt,
56
wl,t=2wl,tx2= φl,t=0,

where ϑ0 is constant and gt is the function of the thermal loading.

Eqs. (25) and (55) gives that:

57
ϑ10,t=ϑ0gt.

From Eqs. (33) and (57), we obtain:

58
ϑ20,t=ϑ0gt+K1sinpzϑ022g2t.

Apply the Laplace transform, we have:

59
w-0,s=2w-0,sx2=0,
60
ϑ-20,s= ϑ0g-t+K1sinpzϑ022g2t¯=Gs.
61
w-l,s=2w-l,sx2=ϑ-2l,s=0.

Then, we obtain the following system of linear equations:

62
i=13Aiki2sinhkil=-Gsε8,
63
i=13Aiki2-ε7sinhkil=0,
64
i=13Aiki2-ε7ki2sinhkil=0.

After solving the above system, then, we get the solutions in the Laplace transform domain as follows:

65
ϑ-2x,s=Gs ε7ε7-k22ε7-k32k12k12-k22k12-k32sinhk1lsinhk1l-x
+ε7-k12ε7-k32k22k22-k12k22-k32sinhk2lsinhk2l-x
+ε7-k12ε7-k22k32k32-k12k32-k22sinhk3lsinhk3l-x,

which gives:

66
θ-x,s=Gs ε7sinpzε7-k22ε7-k32k12k12-k22k12-k32sinhk1lsinhk1l-x
+ε7-k12ε7-k32k22k22-k12k22-k32sinhk2lsinhk2l-x
+ε7-k12ε7-k22k32k32-k12k32-k22sinhk3lsinhk3l-x.

We get the temperature increment φx,z,t from Eqs. (24) and (66) as follows:

67
φx,z,t= -1+1+2K1 L-1θ-x,z,sK1.

The lateral deflection is:

68
w-x,s=Gsε6ε71k12-k22k12-k32sinhk1lsinhk1l-x
+1k22-k12k22-k32sinhk2lsinhk2l-x
+1k32-k12k32-k22sinhk3lsinhk3l-x,

and the deformation takes the form:

69
e-x,s=-z Gsε6ε71k12-k22k12-k32sinhk1lsinhk1l-x
+1k22-k12k22-k32sinhk2lsinhk2l-x
+1k32-k12k32-k22sinhk3lsinhk3l-x.

The stress-strain energy through the beam is given by [12, 37]:

70
ϖx,z,t=i,j=1312σijeij=12σxxx,z,tex,z,t.

Hence, we have:

71
ϖx,z,t=12 L-1σ-xxx,z,sL-1e-x,s,

where L-1 is defined as the inversion of Laplace transform.

To get the complete solutions, we have to determine the function of the thermal loading gt, so, we will consider the thermal loading is ramp-type heat as follows [38]:

72
gt=tt0,0<t<t0,1,tt0.

Hence, we have:

73
Gs= ϑ01-e-s t0s2t0+K1sinpzϑ0221-e-st0-s t0e-st0s3t02.

5. Numerical results and discussion

Now, we will consider a numerical example for which computational results are given. For this purpose, silicon nitride is taken as the thermoelastic material for which we take the following values of the different physical constants [34]: k= 43.5 W/(m K), αT= 2.71×10-6 K-1, ρ= 3200 kg/m3, T0= 293 K, Cυ= 630 J/(kg K), λ0= 217×109 N/m2, μ0= 108×109 N/m2, τ0= 4.32×10-13 s, λ1=μ1= 6.89×10-13s.

The aspect ratios of the beam are fixed asb=h/2 and l/h=6. For the microscale beam, we will assume the length of the beam is in range l (1.0-100)×10-9 m , and the original time t and the relaxation time τ0 are of order 10-12 sec and 10-14 sec, respectively. The figures were carried by using the dimensionless variables for a beam with length l=1.0θ0=1.0, τ0=0.01z=h/2 and t=1.0.

Figs. 1(a-d) represent the temperature increment, the lateral deflection, the stress, and the stress-strain energy distributions, respectively, with various values of the thermal conductivity (K1=0.0, K10.0) when tt0. It has been noted that the parameter K1 has significant effects on the temperature increment distribution, the lateral deflection, the stress, and the stress-strain energy distributions. An increase in the value of the parameter K1 leads to an increase in the values of the temperature increment, lateral deflection, stress, and stress-strain energy.

Fig. 1The state-functions distributions based on the thermal conductivity when t≥t0

The state-functions distributions based on the thermal conductivity when t≥t0

a) The temperature increment distribution

The state-functions distributions based on the thermal conductivity when t≥t0

b) The lateral deflection distribution

The state-functions distributions based on the thermal conductivity when t≥t0

c) The stress distribution

The state-functions distributions based on the thermal conductivity when t≥t0

d) The stress-strain energy distribution

Figs. 2(a-d) show the temperature increment, lateral deflection, stress, and stress-strain energy distributions, respectively, with various values of the thermal conductivity parameter (K1= 0.0, K1 0.0) when t<t0. It has been noted that the parameter K1 has a limited effect on the temperature increment distribution, while it has significant effects on the lateral deflection, the stress, and the stress-strain energy distributions. The values of the peak points of the lateral deflection, stress, and stress-strain energy increase when the thermal conductivity becomes variable.

Fig. 2The state-functions distributions based on the thermal conductivity when t<t0

The state-functions distributions based on the thermal conductivity when t<t0

a) The temperature increment distribution

The state-functions distributions based on the thermal conductivity when t<t0

b) The lateral deflection distribution

The state-functions distributions based on the thermal conductivity when t<t0

c) The stress distribution

The state-functions distributions based on the thermal conductivity when t<t0

d) The stress-strain energy distribution

Figs. 3(a-d) and 4(a-d) show the temperature increment, lateral deflection, stress, and stress-strain energy distributions, respectively, for various values of the mechanical damage variable (υ= 0.0, υ= 0.3) when tt0 and t<t0, respectively. It has been noted that the parameter υ has a limited effect on the temperature increment distribution. The mechanical damage variable has significant impacts on the lateral deflection, stress, and stress-strain energy distributions in the two cases tto and t>to. An increase in the value of the mechanical damage variable leads to a decrease in the values of the peak points of the lateral deflection, stress, and stress-strain energy.

Figs. 5(a-d) and 6(a-d) show the temperature increment, lateral deflection, stress, and stress-strain energy distributions, respectively, for various values of the viscothermoelastic parameter (β1=β2= 0.0,…, β1β20.0) when tt0 and t<t0, respectively. It is observed that the mechanical relaxation times parameters β1 and β2 have a limited effect on the temperature increment distribution, while they have significant effects on the lateral deflection, stress, and stress-strain energy distributions. An increase in the values of the mechanical relaxation times parameters β1 and β2 leads to a decrease in the values of the peak points of the lateral deflection, stress, and stress-strain energy.

Fig. 3The state-functions distributions based on the damage mechanics variable when t≥t0

The state-functions distributions based on the damage mechanics variable when t≥t0

a) The temperature increment distribution

The state-functions distributions based on the damage mechanics variable when t≥t0

b) The lateral deflection distribution

The state-functions distributions based on the damage mechanics variable when t≥t0

c) The stress distribution

The state-functions distributions based on the damage mechanics variable when t≥t0

d) The stress-strain energy distribution

Fig. 4The state-functions distributions based on the damage mechanics variable when t<t0

The state-functions distributions based on the damage mechanics variable when t<t0

a) The temperature increment distribution

The state-functions distributions based on the damage mechanics variable when t<t0

b) The lateral deflection distribution

The state-functions distributions based on the damage mechanics variable when t<t0

c) The stress distribution

The state-functions distributions based on the damage mechanics variable when t<t0

d) The stress-strain energy distribution

Fig. 5The state-functions distributions based on the viscothermoelastic parameters when t≥t0

The state-functions distributions based on the viscothermoelastic parameters when t≥t0

a) The temperature increment distribution

The state-functions distributions based on the viscothermoelastic parameters when t≥t0

b) The lateral deflection distribution

The state-functions distributions based on the viscothermoelastic parameters when t≥t0

c) The stress distribution

The state-functions distributions based on the viscothermoelastic parameters when t≥t0

d) The stress-strain energy distribution

Fig. 6The state-functions distributions based on the viscothermoelastic parameters when t<t0

The state-functions distributions based on the viscothermoelastic parameters when t<t0

a) The temperature increment distribution

The state-functions distributions based on the viscothermoelastic parameters when t<t0

b) The lateral deflection distribution

The state-functions distributions based on the viscothermoelastic parameters when t<t0

c) The stress distribution

The state-functions distributions based on the viscothermoelastic parameters when t<t0

d) The stress-strain energy distribution

Fig. 7The state-functions distributions based on the thickness of the nanobeam when t≥t0

The state-functions distributions based on the thickness of the nanobeam when t≥t0

a) The temperature increment distribution

The state-functions distributions based on the thickness of the nanobeam when t≥t0

b) The lateral deflection distribution

The state-functions distributions based on the thickness of the nanobeam when t≥t0

c) The stress distribution

The state-functions distributions based on the thickness of the nanobeam when t≥t0

d) The stress-strain energy distribution

Fig. 8The state-functions distributions based on the thickness of the nanobeam when t<t0

The state-functions distributions based on the thickness of the nanobeam when t<t0

a) The temperature increment distribution

The state-functions distributions based on the thickness of the nanobeam when t<t0

b) The lateral deflection distribution

The state-functions distributions based on the thickness of the nanobeam when t<t0

c) The stress distribution

The state-functions distributions based on the thickness of the nanobeam when t<t0

d) The stress-strain energy distribution

Figs. 7(a-d) and 8(a-d) show the temperature increment, lateral deflection, stress, and stress-strain energy distributions, respectively, for various values of beam’s thickness zz=h/4,z=h/2 when t0=1 and t0=2, respectively. It has been noted that the beam’s thickness has significant effects on the temperature increment, lateral deflection, stress, and stress-strain energy distributions.

Figs. 9(a-d) and 10(a-d) show the temperature increment, lateral deflection, stress, and stress-strain energy distributions, respectively, for different values of damage mechanics parameter υ= (0.0, 0.3) when t=t0=1.0, for constant and variable thermal conductivity, respectively. In those figures, we discuss the effect of the damage variable to the thermal conductivity.

Figs. 9(a) and 10(a) show that the values of the temperature increment based on variable thermal conductivity are higher than the values of the temperature increment based on constant thermal conductivity for the two cases of damage mechanics variable.

Fig. 9The state-functions distributions when the thermal conductivity is constant

The state-functions distributions when the thermal conductivity is constant

a) The temperature increment distribution

The state-functions distributions when the thermal conductivity is constant

b) The lateral deflection distribution

The state-functions distributions when the thermal conductivity is constant

c) The stress distribution

The state-functions distributions when the thermal conductivity is constant

d) The stress-strain energy distribution

Figs. 9(b) and 10(b) state that the values of the difference between the peak points of the lateral deflection based on the damage and non-damage situations with variable thermal conductivity are higher than the values of the difference between the peak points of the lateral deflection based on the damage and non-damage situations with constant thermal conductivity.

Fig. 9(c) and 10(c) state that the value of the difference between the peak points of the stress based on the damage and non-damage situations with variable thermal conductivity is higher than the values of the difference between the peak points of the stress based on the damage and non-damage situations with constant thermal conductivity.

Figs. 9(d) and 10(d) state that the value of the difference between the peak points of the stress-strain energy based on the damage and non-damage situations with variable thermal conductivity are higher than the values of the difference between the peak points of the stress-strain energy based on the damage and non-damage situations with constant thermal conductivity. Thus, thermal conductivity has a significant effect on the damage mechanics variable and vis versa.

Fig. 10The state-functions distributions when the thermal conductivity is variable

The state-functions distributions when the thermal conductivity is variable

a) The temperature increment distribution

The state-functions distributions when the thermal conductivity is variable

b) The lateral deflection distribution

The state-functions distributions when the thermal conductivity is variable

c) The stress distribution

The state-functions distributions when the thermal conductivity is variable

d) The stress-strain energy distribution

Fig. 1-10 show that the ramp-type heating parameter to has significant effects on the temperature increment, lateral deflection, stress, and stress-strain energy. An increase the value of the ramp-time heat parameter t0 leads to a decrease the values of all the studied functions. Thus, the ramp-time heat parameter t0 plays a vital role in the propagation of the thermal and mechanical waves and it could be used to control the vibration of the microbeam and damping the energy generated in the beam.

Fig. 11(a-d) represents only the lateral deflection for wide range of time t 0.0t2.0 and at distance x=0.2 when t0=1.0 with different cases.

Fig. 11(a) shows the later deflection with respect to the thermal conductivity case, where the solid line represents the case of constant thermal conductivity while the dotted line represents the variable thermal conductivity case. It is noted that the considering of variable thermal conductivity has a significant effect on the lateral deflection distribution. Considering variable thermal conductivity leads to an increase in the value of the lateral deflection.

Fig. 11(b) shows the later deflection with various values of the mechanical damage variable, where the solid line represents the undamaged case, while the dotted line represents the damaged case. It is noted that the mechanical damage variable has a significant effect on the lateral deflection distribution. An increase in the value of mechanical damage variable leads to a decrease in the value of the lateral deflection.

Fig. 11(c) shows the later deflection with various values of the mechanical relaxation time parameter, where the solid line represents the thermoelastic case (non-visco), while the dotted line represents the viscothermoelastic case. It is noted that the mechanical relaxation time parameter has a limited effect on the lateral deflection distribution.

Fig. 11(d) shows the later deflection with various values of the beam’s width z, where the solid line represents the case of beam’s width z=h/4, while the dotted line represents the case of beam’s width z=h/2. It is noted that the beam’s width z has a significant effect on the lateral deflection distribution. An increase in the value of z leads to an increase in the value of the lateral deflection.

Fig. 11The lateral deflection for wide range of time t 0.0≤t≤2.0 and at distance x=0.2 when t0=1.0

The lateral deflection for wide range of time t 0.0≤t≤2.0  and at distance x=0.2 when t0=1.0

a) The lateral deflection when υ=0.3, z=h/2, and β1β20.0

The lateral deflection for wide range of time t 0.0≤t≤2.0  and at distance x=0.2 when t0=1.0

b) The lateral deflection when K10.0, z=h/2, and β1β20.0

The lateral deflection for wide range of time t 0.0≤t≤2.0  and at distance x=0.2 when t0=1.0

c) The lateral deflection when υ=0.3, z=h/2, and K10.0

The lateral deflection for wide range of time t 0.0≤t≤2.0  and at distance x=0.2 when t0=1.0

d) The lateral deflection when K10.0 , υ=0.3, and β1β20.0

Overall, the results in this work agree with the results in other papers, mainly in the cases K1=υ=β1=β2=0 [9, 10, 12, 19, 22, 25, 32, 37].

6. Conclusions

A simply supported viscothermoelastic microbeam has been thermally loaded by ramp-type heating considering variable thermal conductivity and damage mechanics variable.

The ramp-type heating parameter, thermal conductivity, and thickness of the microbeam have significant effects, while the damage mechanics variable and the mechanical relaxation time have a limited impact on the temperature increment. The ramp-type heating parameter, thermal conductivity, thickness of the microbeam, damage mechanics variable, and mechanical relaxation time have significant effects on the lateral deflection, deformation, stress, and stress-strain energy distributions. Thus, the considering of the variability of the thermal conductivity and damage mechanics variable are more significant among studying the microbeams resonators. The thermal conductivity and damage mechanics variable between them have a significant mutual effect.

The ramp-time heat parameter can be used as a tuner to the vibration and the generated entire energy of the microbeam.

References

  • Tzou D. On the thermal shock wave induced by a moving heat source. Journal of Heat Transfer, Vol. 111, Issue 2, 1989, p. 232-238.
  • Tzou D., Transfer M.-T.-M.-H. The Lagging Behavior. Taylor Francis, Washington, 1997.
  • Xu M., Guo J., Wang L., Cheng L. Thermal wave interference as the origin of the overshooting phenomenon in dual-phase-lagging heat conduction, International Journal of Thermal Sciences, Vol. 50, Issue 5, 2011, p. 825-830.
  • Al Huniti N.-S., Al Nimr M. Thermoelastic behavior of a composite slab under a rapid dual-phase-lag heating. Journal of Thermal Stresses, Vol. 27, Issue 7, 2004, p. 607-623.
  • Ho J.-R., Kuo C.-P., Jiaung W.-S. Study of heat transfer in multilayered structure within the framework of dual-phase-lag heat conduction model using lattice Boltzmann method. International Journal of Heat and Mass Transfer, Vol. 46, Issue 1, 2003, p. 55-69.
  • Lee Y.-M., Tsai T.-W. Ultra-fast pulse-laser heating on a two-layered semi-infinite material with interfacial contact conductance. International Communications in Heat and Mass Transfer, Vol. 34, Issue 1, 2007, p. 45-51.
  • Liu K.-C. Numerical analysis of dual-phase-lag heat transfer in a layered cylinder with nonlinear interface boundary conditions. Computer Physics Communications, Vol. 177, Issue 3, 2007, p. 307-314.
  • Ramadan K. Semi-analytical solutions for the dual phase lag heat conduction in multilayered media. International Journal of Thermal Sciences, Vol. 48, Issue 1, 2009, p. 14-25.
  • Alghamdi N. Dual-phase-lagging thermoelastic damping vibration in micro-nano scale beam resonators with voids. International Journal of Multidisciplinary and Current research, Vol. 5, 2017, p. 71-78.
  • Alghamdi N. A., Youssef H. M. Dual-phase-lagging thermoelastic damping in-extensional vibration of rotating nano-ring. Microsystem Technologies, Vol. 23, Issue 10, 2017, p. 4333-4343.
  • Biot M. A. Thermoelasticity and irreversible thermodynamics. Journal of Applied Physics, Vol. 27, Issue 3, 1956, p. 240-253.
  • Youssef H. M., Alghamdi N. Thermoelastic damping in nanomechanical resonators based on two-temperature generalized thermoelasticity theory. Journal of Thermal Stresses, Vol. 38, Issue 12, 2015, p. 1345-1359.
  • Lord H. W., Shulman Y. A generalized dynamical theory of thermoelasticity. Journal of the Mechanics and Physics of Solids, Vol. 15, Issue 5, 1967, p. 299-309.
  • Dhaliwal R. S., Sherief H. H. Generalized thermoelasticity for anisotropic media. Quarterly of Applied Mathematics, Vol. 38, Issue 1, 1980, p. 1-8.
  • Hoang C. M. Thermoelastic damping depending on vibration modes of nano beam resonator. Communications in Physics, Vol. 25, Issue 4, 2015, p. 317.
  • Naik A. K., Hanay M., Hiebert W., Feng X., Roukes M. L. Towards single-molecule nanomechanical mass spectrometry. Nature Nanotechnology, Vol. 4, Issue 7, 2009, p. 445.
  • O’Connell A. D., Hofheinz M., Ansmann M., Bialczak R. C., Lenander M., Lucero E., Neeley M., Sank D., Wang H., Weides M. Quantum ground state and single-phonon control of a mechanical resonator. Nature, Vol. 464, Issue 7289, 2010, p. 697.
  • Van Beek J., Puers R. A review of MEMS oscillators for frequency reference and timing applications. Journal of Micromechanics and Microengineering, Vol. 22, Issue 1, 2011, p. 013001.
  • Sharma J., Grover D. Thermoelastic vibrations in micro-/nano-scale beam resonators with voids. Journal of Sound and Vibration, Vol. 330, Issue 12, 2011, p. 2964-2977.
  • Sun Y., Saka M. Thermoelastic damping in micro-scale circular plate resonators. Journal of Sound and Vibration, Vol. 329, Issue 3, 2010, p. 328-337.
  • Lifshitz R., Roukes M. L. Thermoelastic damping in micro-and nanomechanical systems. Physical review B, Vol. 61, Issue 8, 2000, p. 5600.
  • Al Lehaibi E.-A., Youssef H. M. Vibration of gold nano-beam with variable young’s modulus due to thermal shock. World Journal of Nano Science and Engineering, Vol. 5, Issue 4, 2015, p. 194.
  • Kidawa Kukla J. Application of the Green functions to the problem of the thermally induced vibration of a beam. Journal of Sound and Vibration, Vol. 262, Issue 4, 2003, p. 865-876.
  • Boley B. A. Approximate analyses of thermally induced vibrations of beams and plates. Journal of Applied Mechanics, Vol. 39, Issue 1, 1972, p. 212-216.
  • Manolis G., Beskos D. Thermally induced vibrations of beam structures. Computer Methods in Applied Mechanics and Engineering, Vol. 21, Issue 3, 1980, p. 337-355.
  • Al Huniti N.-S., Al Nimr M., Naji M. Dynamic response of a rod due to a moving heat source under the hyperbolic heat conduction model. Journal of Sound and Vibration, Vol. 242, Issue 4, 2001, p. 629-640.
  • Biot M. A. Theory of stress‐strain relations in anisotropic viscoelasticity and relaxation phenomena. Journal of Applied Physics, Vol. 25, Issue 11, 1954, p. 1385-1391.
  • Biot M. A. Variational principles in irreversible thermodynamics with application to viscoelasticity. Physical Review, Vol. 97, Issue 6, 1955, p. 1463.
  • Drozdov A. A constitutive model in finite thermoviscoelasticity based on the concept of transient networks. Acta Mechanica, Vol. 133, Issues 1-4, 1999, p. 13-37.
  • Ezzat M. A., El-Karamany A.-S. The relaxation effects of the volume properties of viscoelastic material in generalized thermoelasticity. International Journal of Engineering Science, Vol. 41, Issue 19, 2003, p. 2281-2298.
  • Carcione J. M., Poletto F., Gei D. 3-D wave simulation in anelastic media using the Kelvin-Voigt constitutive equation. Journal of Computational Physics, Vol. 196, Issue 1, 2004, p. 282-297.
  • Grover D. Viscothermoelastic vibrations in micro-scale beam resonators with linearly varying thickness. Canadian Journal of Physics, Vol. 90, Issue 5, 2012, p. 487-496.
  • Grover D. Damping in thin circular viscothermoelastic plate resonators. Canadian Journal of Physics, Vol. 93, Issue 12, 2015, p. 1597-1605.
  • Grover D. Transverse vibrations in micro-scale viscothermoelastic beam resonators. Archive of Applied Mechanics, Vol. 83, Issue 2, 2013, p. 303-314.
  • Grover D., Seth R. Viscothermoelastic micro-scale beam resonators based on dual-phase lagging model. Microsystem Technologies, Vol. 24, Issue 3, 2018, p. 1667-1672.
  • Claussen N., Hasselman D., Heller R. Thermal Stresses in Severe Environments. Plenum Press, New York, 1980, p. 381-95.
  • Youssef H. Vibration of gold nanobeam with variable thermal conductivity: state-space approach. Applied Nanoscience, Vol. 3, Issue 5, 2013, p. 397-407.
  • Youssef H., El Bary A. Thermal shock problem of a generalized thermoelastic layered composite material with variable thermal conductivity. Mathematical Problems in Engineering, Vol. 2006, 2006, p. 87940.
  • Youssef H. M., El-Bary A. Mathematical model for thermal shock problem of a generalized thermoelastic layered composite material with variable thermal conductivity. Computational Methods in Science and Technology, Vol. 12, Issue 2, 2006, p. 165-171.
  • Zenkour A. M., Abbas I. A. A generalized thermoelasticity problem of an annular cylinder with temperature-dependent density and material properties. International Journal of Mechanical Sciences, Vol. 84, 2014, p. 54-60.
  • Gross D., Seelig T. Fracture Mechanics: with an Introduction to Micromechanics. Springer, 2017.
  • Rezaiee Pajand M., Kazemiyan M. S., Aftabi A. S. Static damage identification of 3D and 2D frames. Mechanics Based Design of Structures and Machines, Vol. 42, Issue 1, 2014, p. 70-96.
  • Öchsner A. Continuum damage mechanics, Continuum Damage and Fracture Mechanics. Springer, 2016.
  • Voyiadjis G. Z. Handbook of Damage Mechanics: Nano to Macro Scale for Materials and Structures. Springer, 2015.
  • Yao Y., He X., Keer L. M., Fine M. E. A continuum damage mechanics-based unified creep and plasticity model for solder materials. Acta Materialia, Vol. 83, 2015, p. 160-168.
  • Voyiadjis G. Z., Kattan P. I. Introducing damage mechanics templates for the systematic and consistent formulation of holistic material damage models. Acta Mechanica, Vol. 228, Issue 3, 2017, p. 951-990.
  • Khatir A., Tehami M., Khatir S., Wahab Abdel M. Multiple damage detection and localization in beam-like and complex structures using co-ordinate modal assurance criterion combined with firefly and genetic algorithms. Journal of Vibroengineering, Vol. 18, Issue 8, 2016, p. 5063-5073.
  • Hetnarski R. B., Eslami M. R., Gladwell G. Thermal Stresses: Advanced Theory and Applications. Springer, 2009.
  • Kwok K. Shape recovery of viscoelastic beams after stowage. Mechanics Based Design of Structures and Machines, Vol. 43, Issue 1, 2015, p. 95-111.
  • Xu L., Yang Q. Multi-field coupled dynamics for a micro beam. Mechanics Based Design of Structures and Machines, Vol. 43, Issue 1, 2015, p. 57-73.
  • Sahmani S., Madyira D. M. Nonlocal strain gradient nonlinear primary resonance of micro/nano-beams made of GPL reinforced FG porous nanocomposite materials. Mechanics Based Design of Structures and Machines, 2019, https://doi.org/10.1080/15397734.2019.1695627.
  • Haskul M. Elastic state of functionally graded curved beam on the plane stress state subject to thermal load. Mechanics Based Design of Structures and Machines, 2019, https://doi.org/10.1080/15397734.2019.1660890.

About this article

Received
05 May 2020
Accepted
06 August 2020
Published
21 November 2020
SUBJECTS
Fault diagnosis based on vibration signal analysis
Keywords
vibration
viscothermoelastic
microbeam
silicon nitride
ramp-type heat
thermal conductivity
damage