Investigation of temperature modes in thermosensitive non-uniform elements of radioelectronic devices Текст научной статьи по специальности «Физика»

PA^IOE^EKTPOHIKA TA TE^EKOMYfflKA^Ï

PA^HOMEKTPOHHKA H TЕ^ЕKOММУНHKAЦHH

RADIO ELECTRONICS AND TELECOMMUNICATIONS

UDC 536.24

INVESTIGATION OF TEMPERATURE MODES IN THERMOSENSITIVE NON-UNIFORM ELEMENTS OF RADIOELECTRONIC DEVICES

Havrysh V. I. - Dr. Sc., Professor of Department of Software, Lviv Polytechnic National University, Lviv, Ukraine.

Baranetskij Ya. O. - PhD, Assosiate Professor of Department of Mathematics, Lviv Polytechnic National University, Lviv, Ukraine.

Kolyasa L. I. - Doctor of Philosophy, Senior Lecture of Department of Mathematics, Lviv Polytechnic National University, Lviv, Ukraine.

ABSTRACT

Context. The non-linear boundary value problem of heat conduction for a thermosensitive non-homogeneous strip-shaped element of a radio-electronic system with a through inclusion has been solved whose analytical-numerical solution enables us to analyze temperature regimes in the element.

Objective. Is to develop such a method of linearization of mathematical model of heat conduction which enables us to obtain analytical numerical solution of the corresponding non-linear boundary value problem for determination of temperature field in elements of radio electronic devices, which are geometrically represented by a thermosensitive plate with a through inclusion.

Method. A linearizing function which enables us to partially linearize the initial non-linear mathematical model of heat conduction for a thermosensitive non-homogeneous element of a radio electronic system in the form of "plate-inclusion" structure has been suggested. The introduced piece-wise linear approximation of temperature on plate-inclusion interfaces has enabled us to completely linearize the corresponding partially linearized boundary value problem relative to the linearizing function. After this, it became possible to apply Fourier's integral transformation to the obtained linear problem with respect to one of the spatial coordinates, as well as to determine the linearizing function. The linear dependence of the coefficient of heat conductivity on temperature for structure materials with the use of the linearizing function has been considered. By solving the boundary value problem, the formulae for determination of temperature field in the "plate-inclusion" thermosensetive structure have been obtained.

Results. The obtained formulae for determination of temperature field in a thermosensitive non-homogeneous element of radio electronic system were used to create the software which enables us to obtain distribution of value of temperature and to analyze temperature regimes.

Conclusions. A mathematical model for the calculation for the temperature field in a "plate-inclusion" thermosensitive structure is adequate to the actual physical process, because no jump of temperature at "plate-inclusion" interfaces is observed. The numerical results for the chosen materials under linear dependence of the coefficient of thermoconductivity on temperature differ by 7% from the results which are obtained for constant coefficient of heat conductivity. Prospect of further investigation will consider more complicated geometric representation of elements of radio electronic systems.

KEYWORDS: temperature, heat conduction, nonlinear boundary-value problem, isotropic infinite thermosensitive plate with insulated faces, through inclusion, perfect thermal contact, heat flow.

NOMENCLATURE

25 is the thickness of a heat-sensitive isotropic with respect to its thermophysical parameters plate ;

2h is the length of a foreign through inclusion;

Q0 is the domain of the inclusion;

q1 is the domain of the plate (beyond the inclusion);

K is the boundary surface of the plate, the system is heated by a heat stream; heat stream, whose surface density is q0 = const ;

K+ is the part of the surface of the plate beyond the inclusion and is thermo-insulated;

K+h is the boundary surfaces of the inclusion;

t(x,y) is the distribution of temperature with respect to the coordinates of space;

X(x, t) is the coefficient of heat conductivity of the et-erogeneous heat-sensitive plate;

Xo(t ) is the coefficient of het conductivity of the material of the plate;

X1 (t) is the coefficient of het conductivity of the material of the inclusion;

S± (Z) are the asymmetric unit functions;

x, y) is the linearizing function; A is the Laplace operator in the Cartesian system of coordinates;

tj are the unknown approximated values;

s± (z) is the asymmetric Dirac delta-function;

y) is the transform formula s(x, y); \ is the parameter of Fourier integral transformation; n is the number of intervals in the partition ]_i; i [ ;

X °m are basic coefficients of heat conductivity of the materials of the inclusion (m = 0); and the plate (m = 1); km is the temperature coefficients of heat conductivity

of the materials of the inclusion (m = 0); and the plate (m = 1);

t* is the dimensionleess temperature; x , y* are the dimensionless coordinates of space;

X( x) is the coefficient of heat conductivity of the

heterogeneous plate for the corresponding linear model.

INTRODUCTION Composite materials, development of which is a liding problem of modern materials science, have become especially valuable in radioelectronic device's. The emergence of new composite materials with improved physic-mechanical performance will contribute to creation of new technologies in air-space, naval, power-generating, and electronic branches, in mechanical engineering, and transport. Among the composite materials an important place occupy structures with foreign inclusions, which are widely used structures of microelectronic devices, in particular, in integral sensors for monitoring of temperature and humidity, in light emitting elements for dynamic LED-backlighting etc. Because the aforesaid structures operate in a wide range of temperatures, their high performances cause the necessity of considering and solving of the problems which are non-linear due to the dependence of thermophysical characteristics of the materials on the temperature of the structure and on the conditions of heat exchange, on the temperature of their surfaces; because the calculations of the temperature fields which are performed on the basis of the linear mathematical models of heat conduction processes do not always give us satisfactory results [1]. Therefore, for creation of the most adequate to actual process mathematical model, it is necessary to take into account the dependencies of thermo-physical characteristics of materials on temperature, on density of surface heat fluxes, and on intensity of internal heat sources; change of the shape of the body and possible phase and structural transformations are also to be taken into account [2, 3].

The object of study is the process of heat conduction for an element of radio electronic devices presented by a thermosensitive plate with a foreign inclusion.

© Havrysh V. I., Baranetskij Ya. O., Kolyasa L. I., 2018 DOI 10.15588/1607-3274-2018-3-1

In order to create a non-linear mathematical model of the process of heat-conduction in a thermosensitive non-homogeneous element of radio electronic systems, theory of generalized function has been used; and in order to construct a numerical-analytical solution of the corresponding boundary value problem, integral Fourier transformation and methods of solving partial differential equations were used.

The subject of study consists of non-linear mathematical models of heat-conduction process and methods of construction of numerical analytical solution of the corresponding boundary value problems for thermosensi-tive elements of radio electronic devices of piecewise-homogeneous structure.

The purpose of the work is the development of such a method of linearization of mathematical model of heat conduction which enables us to obtain the analytical-numerical solution of the corresponding non-linear boundary value problem for determination of temperature field in elements of radio electronic devices, which are geometrically represented by a thermosensitive plate with a through inclusion.

1 PROBLEM STATEMENT Let us consider a heat-sensitive isotropic with respect to its thermophysical parameters plate of 28 thickness with thermo-insulated faces |z| = ±8. The plate contains a

foreign through inclusion of 2h length. The plate is referenced to the Cartesian system of coordinates (Oxyz) with whose origin is in the centre of the inclusion. In the domain of the inclusion q0 = {(x, -l, z): |x| < h,|y| < l,\z| < 8}

of the boundary surface K- ={(x, -l, z ):|x| <<»,| z| <8} of

the plate, the system is heated by a heat stream, whose surface density is q0 = const, and the other part of the surface of the plate beyond the inclusion and the surface K+ ={(x,l,z): |x| <<»,|z| <8} are thermo-insulated. At the

boundary surfaces of the inclusion

K±h = {(±h, y, z): |y| < l, |z| <8} perfect heat contact takes place

|x| = h (0 - for inclusion,

/0 t, .

10 (t >£■

X. (t>d/1 for

u ' dx

1 - for plane) (fig. 1).

K

21

K _

%

A'

X

y

Figure 1 - Section of isotropic thermosensitive plate with foreign through inclusion by the plane z = 0

In the given structure, it is necessary to determine the distribution of temperature t(x, y) with respect to the coordinates of space.

2 REVIEW OF THE LITERATURE

The determination of temperature conditions, both of homogeneous and non-homogeneous structures, draws attention of many researchers [4, 10].

In the work [11], a mathematical model of quasi-stationary temperature field in a continuous cylinder of rotation made of a composite material with non-linear boundary conditions is developed; the dependence of thermophysical parameters of material on temperature being taken into account. The obtained analytical expressions for determination of temperature fields enable us to select the composition of composite materials for cylinder-shaped parts in order to increase there are service life.

One-dimentional (plane, cylindrically symmetric, and spherically symmetric) non-linear problems of heat-conduction are considered for a given flux of heat at the origin of coordinate systems in the form of exponential function depending on time. Approximation solutions of the aforesaid problems have been obtained, and there are convergence has been analyzed [12].

An analytical solution of the non-linear problem of heat-conduction has been constructed on the basis of integral method of heat balance [13]. In order to improve the accuracy of the solution, the temperature function is approximated by high-degree polynomials. To determine the coefficients of the polynomials, additional boundary conditions have been introduced. It is shown that the introduction of additional boundary conditions, as early as in the second approximation, leads to considerable improvement of the accuracy of the solution of the problem.

In the work [14], numerical-analytical solution of the non-stationary problem of heat-conduction for a hollow sphere has been constructed; thermophysical parameters of materials of the sphere being variable with temperature. In the particular case, a solution for a continuous sphere has been obtained. On the basis of variational approach, a non-linear mathematical model of heat-conduction process have been constructed for 2D-space with a thin inclusion, this anable us to take into account small thickness of a thick inclusion. For linearization of the stated problem, the Newton-Raphson's methods has been applied. Discreatization with respect to temporal variable has been carried out in accordance with intermediate point scheme. The Discreatization statement of the problem is presented in the form of minimization of a functional [15].

Investigations of temperature regimes for structural thermosensitive elements of a piece-wise homogeneous structure have been carried out in the works [16-19].

Article [20] solved a non-stationary problem on thermal conductivity and thermoelasticity for functional-gradient thick-wall spheres. Thermal-physical and ther-moelastic parameters of materials, except for Poisson

coefficient, are arbitrary functions of the radial coordinate.

Axisymmetric stationary problem on thermal conductivity and thermoelasticity of the hollow functionally gradient areas relative to the heat source was considered. The solutions are obtained as functions from spatial coordinates for temperature, the displacement component vector and stress tensor by using boundary conditions for radial and angular coordinates [21].

In paper [22], a thermoelasticity solution for steady state response of thick cylinders which are subjected to pressure and external heat flux in inner surface is presented.

A general solution for the one-dimensional steady-state thermal and mechanical stresses in a hollow thick sphere made of porous functionally graded material (FGPM) is presented in [23].

The boundary value non-linear problem of heat conduction is determined below, the technique of its linearization and calculation formulae for determination of the temperature field in a heat-sensitive plate with a through inclusion which is heated by a heat stream concentrated at a surface of the inclusion is given. Numerical analysis under the condition of linear dependence of the coefficient heat conductivity on temperature is performed.

3 MATERIALS AND METHODS

In the given structure (in problem statement), it is necessary to determine the distribution of temperature t (x, y)

with respect to the coordinates of space which is obtained from solving the nonlinear equation of heat conduction [24, 25]

д r(W ч dt д —[x( x, t )—] +—

дх дх dy

x( x, t )—

ду.

= 0

(1)

with the boundary condition

4 i = 0,

3t_

дх

= 0,

ду

= 0,

^ (t >f

y=l

= —% s- (h — x| ),

(2)

y=-i

where x(x,t) = x 1(t) + [x 0(t)_ x1 (t)] S_(h_|x|) is the

coefficient of heat conductivity of the heterogeneous heat-sensitive plate; x0 (t ), X1 (t ) _ are the coefficients of hat

conductivity of the materials of the plate and the inclusion, respectively;

1, z > 0,

S±(Z) = 0,5 + 0,5, Z = 0, is the asymmetric unit func-

0, Z < 0

tions [26].

Let us introduce the linearizing function

t (x y) t{ x y)

S = j X^Ç)dÇ + S-(x + h) j [xo(Z)-

0 t (—h, y)

t ( x, y )

-X1(Z)]+s+ (x - h) j [x 0(Z) -^1(0] dZ,

t (h, y )

(3)

having differentiated which with respect to the variables x and y, we obtain

x(t, x )» =d», V dx dx

X",x )| = f + {[x 0("-x ■(' )]|

S_ ( x + h) -

-{[X 0(t )-x 1(t )]f }

x=-h

S+ (x - h).

(4)

Taking into account the expression (4), the original equation (1), assumes the following form:

■» + !; {[X 0«>-X 1<t ) ]|}

i {[X 0(t )-X.(' ) ]|

S- (x + h)-

=-h

(5)

S+ ( x - h) = 0,

where A =

d 2 d 2

dx2 dy2

is the Laplace operator in the

as

¥

= 0, s

5»

dy

y=l

= -1 q0S-(h - lxl )+

0 d»

|x|= 0, &

= 0,

(x 0(t) -X 1(t)))

x=-h

y=-i

S- ( x+h) -

(X 0(t ) -X 1(t )))

(6)

(7)

x=h

y=-i

S+ ( x - h)\

t (±h, y ) = t1 +X(t,+1 - tj )S- ( y - y. ),

(8)

j=1

Having substituted the expression (8) into the relations (6), (7), we obtain the following partial linear differential equation relating the linearizing function

A» = -51(tj+1 -1j )[X0(tj+1) - X1 (tj+1 )] S- (h - |x|)5-(y - yj ) (9)

j=1

with the boundary condition

= -<?0 S- (h -| x| ).

9y

y=-l

(10)

Here s± (ç) = dS± (ç) is the asymmetric Dirac delta-

d ç

function [26].

Having applied Fourier integral transformation with respect to the coordinate x to the equation (9) and to the boundary condition (6), (10), we obtain the following ordinary differential equation with constant coefficients

Cartesian system of coordinates.

The boundary conditions with the use of the relation (3), are written as follows:

-1 y 1JV"

0 )'„-, I

Figure 2 - Approximation of the function t (±h, y )

d2s

-^HïT -1))0(0+1)-

(11)

The linearizing function (3) has enabled us to reduce the nonlinear boundary value problem (1), (2) to the partially linearized equation with discontinuous coefficients (5) and to the boundary conditions (6), (7).

Let us approximate the function t(±h, y) in the form of Fig. 2 to the expression

-^j+,)]8-(y -y). with the boundary conditions

d S

dy

= 0,

y=l dy

2 — sin h\,

n i

(12)

y=-l

x is the transform formula

where yj e ]_l; l [; y1 < y2 < — < yn-1; n is the number of intervals in the partition ]_i;i [ ; t (j = 1, n) are the unknown approximated values of temperature.

© Havrysh V. I., Baranetskij Ya. O., Kolyasa L. I., 2018 DOI 10.15588/1607-3274-2018-3-1

where 9 = _L j Se^d

x, y); £ is the parameter of Fourier integral transformation; i2 =-1.

Having solved the problem (11), (12), then having applied the inverse to its solution Fourier integral transformation, we obtain the expression of the function 9

ô = n J 7 sin h7 cos x7 ] Z (tj+1 - tj ) P-0 (tj+1 ) -n 0 7 [j=1

-'•1('j+1) ] [ ^ - yj)-

- yj - y- )]+f ^} d7.

(13)

Having substituted the temperature dependence of the coefficient of heat conductivity of the materials of the plate and that of the inclusion into the relations (3), (13), after some transformations, we obtain the system of nonlinear equations for determining the unknown approximated values of the temperature tj (j = 1, n).

For the given system, the unknown temperature field we determine by means of the obtained nonlinear equation with the use of the relations (3), (13) after substitution into them the concrete expressions of the dependence of coefficient of heat conductivity of the structural materials on temperature.

4 EXPERIMENTS

To solve many practical problems, the following dependence of the coefficient of thermal conductivity on temperature is used [27, 28]:

+P>-1' )(^ *» - y- )+

+(1 - ch7(y-y- ))S- (y- y j ))] [sin 7( x + h) - sin 7( x - h)] + 2q0ch7(y -1)

7sh27

-cos 7x sin 7h}d 7.

x = xm (1- kmt),

(14)

The given expression enables us to obtain numerical results represented in Fig. 3b for the considered structure materials with the constant coefficients of their heat conductivities

= 13,4W/ (Km), X0 = 419W/(Km).

The dependence of the dimensionless temperature t

* II

on the special dimensionleess coordinate x = x / h and y* = y / h is shown in Fig. 3. Note, that the maximal value of temperature is reached in the domain of the action of the concentrated heat stream, and at the edge K±h of the inclusion, the satisfaction of the conditions of perfect heat contact (no temperature jump) is observed; this is in accordance to the considerate mathematical model.

The change in dimensionleess temperature t' depending on the dimensionless coordinates y* for x* = 0 (Fig. 4a) and on x* for y* = 0 (Fig. 4b) is shown in Fig. 4. The behavior of the curves indicates the adequacy of the mathematical model to the actual physical process, since at the surface of the inclusion the perfect contact is observed (no jump in temperature).

where Xm, km are basic and temperature coefficients of heat conductivity of the materials of the inclusion (m=0) and the plate (m=1).

Numerical analysis of the dimensionless temperature

t* = X0t / (q0h) is carried out for the following initial

data: material of the layer is VK94-I ceramics, the material of the inclusion is silver, n=10 is the number of intervals in the partition ]-/; I [ ; I / h = 1. In the temperature interval [20°C;1230°C] , the temperature dependencies of the the coefficient of heat conductivity for the given materials are expressed as follows:

I", Ci

X1(t ) = 13,67—(1 - 0,00064—t ), Km K

X0 (t ) = 422,54—(1 - 0,00031—t ), Km K

(15)

which is a partial case of the relation (14).

Let us express the temperature field for the corresponding linear model from the following relation

x) t (x, y)=1 f !{(x0+ n 0 %

Figure 3 - Dependence of dimensionless temperature t on dimensionless coordinates x and y for linearly variable (a) and constant (b) coefficient of heat conductivity of materials of structure

Figure 4 - Dependence of dimensionless temperature t on di* * _ n

mensionless coordinates y for x — 0 (a)

* * „

and x for y — 0 (b). 5 RESULTS

Having taken into account the relationship (14), from the expressions (3), (13), we obtain the formulae for determining the temperature t in the domain of the inclusion

t = -

( ( \ 2S „

1 - \ 1 - k0 ^ K 0 V u /

I /

(16)

in the domain Q — {(x,y, z): |x| > h, |y| < l, |z| <8} of the plate (beyond the inclusion)

t = -

1 - 1 --

2k, S

-0

(17)

Here

2 - V--1 (2 - ¥)

x=h

x=h

1 -

1 - 21S

K01 x=h

The formulae (16), (17) completely describe the temperature field in a heat-sensitive layer with a foreign-through inclusion (fig. 1).

The number of intervals «=10 in the partition ]-/; l[

for the given thermophysical (basic and temperature coefficients of heat conductivity for the materials of the plate and that of the inclusion) and geometric (length and width of the inclusion, width of the plate) parameters of the structure enables us to perform calculations accurate to s = 10-6.

6 DISCUSSION

In the course of development and investigation of linear and non-linear mathematical models of heat conduction process for structures which are geometrically de© Havrysh V. I., Baranetskij Ya. O., Kolyasa L. I., 2018 DOI 10.15588/1607-3274-2018-3-1

scribed by means of the presented piece-wise homogeneous structures, it is detected that the numerical results of temperature field for the considered materials in the case of constant coefficient of thermal conductivity differ by 7% from those for linear variable. This indicates that the consideration of the dependence of thermophysical parameters on temperature of materials of structural elements of complex systems are important, since results which are obtained with the use of nonlinear models are more accurate. The conscideration of piece-wise homogeneous structure of elements of the object is also important in the investigations presented; this makes the solving of corresponding linear and non-linear boundary value problems considerably more complicated; but solution of such problems describes temperature distribution in more adequate way to the actual process.

CONCLUSIONS

The introduced linearizing function enables us to partially linearize the original nonlinear boundary value problem, and the suggested piecewise-linear approximation of temperature at the boundary surfaces of the inclusion to completely linearize the equation and the boundary condition. Due to this it became possible to apply Fourier integral transformation to the obtained linear problem relating the suggested function and to obtain its analytical-numerical solution. The linear temperature dependence of coefficient of heat conductivity of the materials of the inclusion and that of the plate is considered. On the basis of this, the formulae for the calculation of the values of the temperature in the considered "plate-inclusion" structure are given. The obtained results for the chosen materials under linear dependence of coefficient of heat conductivity on temperature (Fig. 3 a) differ from the results obtained under the condition of constant coefficient (Fig. 3b) by 7%. Their inconsiderable difference is accounted for by the fact that the actual values of the temperature coefficient of the heat conductivity for the considered materials are small.

REFERENCES

1. Savatorova V. L. Reshenie uravnenija teploprovodnosti v neodnorodnoj srede s uchetom temperaturnoj zavisimosti kojefficienta teploprovodnosti, Obozrenie prikladnoj i promyshlennoj matematiki, 2009, 17, No. 1, pp. 135-137.

2. Kudinov V. A. Analiz nelinejnoj teploprovodnosti na osnove opredelenija fronta temperaturnogo vozmushhenija, Teplofizika vysokih temperatur, 2006, 44, No. 3, pp. 577585.

3. Kudrjashov N. A. Priblizhennye reshenija odnomernyh zadach nelinejnoj teploprovodnosti pri zadannom potoke, Zhurnal vychislitel'noj matematiki i matematicheskoj fiziki, 2007, 47, No. 1, pp. 110-120.

4. Carpinteri A., Paggi M. Thermoelastic mismatch in nonhomogeneous beams, J. Eng. Math, 2008, 61, No. 2-4. pp. 371-384.

5. Noda N. Thermal stresses in materials with temperature-dependent properties, Appl. Mech. Rev,, 1991, 44, pp. 383397.

0

t

6. Noda N. Thermal stresses in materials with temperature-dependent properties, Appl. Mech. Rev., 1991, 44, pp. 383— 397.

7. Otao Y. Tanigawa O., Ishimaru O. Optimization of material composition of functionality graded plate for thermal stress relaxation using a genetic algorithm, J. Therm. Stresses, 2000, 23, pp. 257-271.

8. Tanigawa Y., Akai T., Kawamura R. Transient heat conduction and thermal stress problems of a nonhomogeneous plate with temperature-dependent material properties, J. Therm. Stresses, 1996, 19, No. 1, pp. 77-102.

9. Tanigawa Y., Otao Y. Transient thermoelastic analysis of functionally graded plate with temperature-dependent material properties taking into account the thermal radiation, Nihon Kikai Gakkai Nenji Taikai Koen Ronbunshu, 2002, 2, pp. 133-134.

10. Yangian Xu., Daihui Tu. Analysis of steady thermal stress in a ZrO2/FGM/Ti-6Al-4V composite ECBF plate with temperature-dependent material properties by NFEM, 2009-WASE Int. Conf. on Informa. Eng., Vol. 2-2, pp. 433-436.

11. Golicyna E. V., Osipov Ju. R. Kvazistacionarnaja trehmernaja zadacha teploprovodnosti vo vrashhajushhemsja sploshnom cilindre iz kompozicionnogo materiala s nelinejnymi granichnymi uslovijami, Konstrukcii iz kompozicionnyh materialov, 2007, No. 4, pp. 47-58.

12. Kudrjashov N. A., Chmyhov M. A. Priblizhennye reshenija pervoj i vtoroj kraevyh zadach nelinejnoj teploprovodnosti na polubeskonechnoj prjamoj, Inzhenernaja fizika, 2007, No. 3, pp. 12-15.

13. Kudinov V. A., Averin B. V., Stefanyuk E. V., Nazarenko S. A. Analysis of nonlinear heat conduction based on determining the front of temperature perturbation, High Temperature, 2006, 44, № 4, pp. 574-583.

14. Popovich V. S., Ivankiv K. S. Nelinijna zadacha teploprovidnosti dlja kuli z teploobminom, VisnikL'viv. un-tu. Ser. Prikl. matematika ta informatika, 2002, No. 5, pp. 136-144.

15. Savula Ja. G., Djakonjuk L. M. Doslidzhennja variacijnoï zadachi teploprovidnosti u bagatosharovih seredovishhah z tonkimi vkljuchennjami, Visnik LNU im. I. Franka, Ser. prikl. matem. ta informat, 2000, No. 3, pp. 125-130.

16. Gavrish V. I., Fedasjuk D. V. Modeljuvannja temperaturnih rezhimiv u kuskovo-odnoridnih strukturah. L'viv: V-vo Nac. un-tu «L'vivs'ka politehnika», 2012, 176 p.

17. Gavrysh V. I. Chislenno-analiticheskoe reshenie nelinejnoj stacionarnoj zadachi teploprovodnosti dlja beskonechnoj termochuvstvitel'noj mnogoslojnoj plastiny, Jelektronnoe modelirovanie, 2014, 36, No. 3, pp. 59-70.

18. Gavrish V. I. Doslidzhennja temperaturnih rezhimiv u termochutlivij plastini z chuzhoridnim naskriznim vkljuchennjam, Fiziko-matematichne modeljuvannja ta informacijni tehnologiï. Naukovij zbirnik, 2013, Vipusk 18, pp. 43-50.

19. Gavrish V. I. Nelinijna krajova zadacha teploprovidnosti dlja sharuvatoï plastini z vkljuchennjam, Fiziko-himichna mehanika, 2015, 51, No. 3, pp. 32-38.

20. Bayat A., Moosavi H., Bayat Y. Thermo-mechanical analysis of functionally graded thick spheres with linearly time-dependent temperature [Text], Scientia Iranica, 2015, Vol. 22, Issue 5, pp. 1801-1812.

21. Mohazzab A. H., Jabbari M. Two-Dimensional Stresses in a Hollow FG Sphere with Heat Source [Text], Advanced Materials Research, 2011, Vol. 264-265, pp. 700-705. DOI: 10.4028/www.scientific.net/amr.264-265.700

22. Ghannad M., Yaghoobi M. P. A thermoelasticity solution for thick cylinders subjected to thermo-mechanical loads under various boundary conditions, International Journal of Advanced Design & Manufacturing Technology, 2015, Vol. 8, No. 4, pp. 1-12.

23. Jabbari M., Karampour S., Eslami M. R. Radially symmetric steady state thermal and mechanical stresses of a poro FGM hollow sphere, International Scholarly Research Network ISRN Mechanical Engineering, 2011, Vol. 2011, pp. 1-7. Article ID 305402. DOI: 10.5402/2011/305402 24.

24. Podstrigach Ja. S., Lomakin V. A., Koljano Ju. M. Termouprugost' tel neodnorodnoj struktury. Moscow, Nauka, 1984, 368 p.

25. Koljano Ju. M. Metody teploprovodnosti i termouprugosti neodnorodnogo tela. Kiev, Naukova dumka, 1992, 280 p.

26. Korn G., Korn T. Spravochnik po matematike dlja nauchnyh rabotnikov i inzhenerov. Moscow, Nauka, 1977, 720 p.

27. Lomakin V. A. Teorija uprugosti neodnorodnyh tel. Moscow, Izd-vo Mosk. un-ta, 1976, 376 p.

28. Berman R. Teploprovodnost' tverdyh tel. Moscow, Mir, 1979, 288 p.

Received 11.06.2017.

Accepted 18.01.2018.

y,qK 536.24

flOC^ia^EHHa TEMnEPATyPHHX PE^HMIB y HEOflHOPiaHHX TEPMOWT.HBHX E.EMEHTAX

PAajOE.nEKTPOHHHX nPHCTPOÏB

TaBpnm B. I. - g-p TexH. HayK, npo^ecop Ka^egpu nporpaMHoro 3a6e3neneHHfl, Ha^oHamtHHH ymBepcmeT «JltBÎBCtKa nomiTexHiKa», JtBiB, yKpaïHa.

Eapaiieui.Kiiii a O. - KaHg. $i3.-MaT. HayK, go^HT Ka^egpu bh^oï MareMaTHKH, Ha^oHamtHHH yHÏBepcHTeT «JtBiBctKa nomiTexHiKa», Jlbïb, yKpaïHa.

Koaaca .. I. - KaHg. $i3.-MaT. HayK, CTapmuft BHKmagan Ka^egpu bh^oï MaTeMaTHKH, Ha^oHamtHHH yHiBepcHTeT « JtBÎBCtKa nomiTexHiKa», JtBiB, yKpaïHa.

AHOTAUja

AKTyaabHÏCTb. Po3B'a3aHo HemimftHy KpaftoBy 3agany TenmonpoBigHocri gna: TepMonyTmHBoro HeogHopigHoro emeMerna pagioemeKTpoHHoï CHCTeMH y BHrmagi CMyru 3 HacKpi3HHM BKmroneHHAM, aHamiTHHHo-HHcmoBHH po3B'a3oK aKoï gae 3Mory aHami3yBaTH y HtoMy TeMnepaTypHi pernMH. MeTa po6oTH - po3po6Ka MeTogy miHeaproa^ï HerniHiftHoï MaTeMamnHoï Mogemi TenmonpoBigHocTi, akhh gae 3Mory oTpHMam aHamiTHHHo-HHcmoBHH po3B'a3oK BignoBigHoï HerniHiftHoï KpaftoBoï 3agani gma BH3Ha-

чення температурного поля в елемент радюелектронних пристро!в, яю геометрично зображено термочутливою пластиною з наскр1зним включенням.

Метод. Запропоновано лшеаризуючу функщю, яка дае змогу частково лшеаризувати вихщну нелшшну математичну модель теплопров1дносл для термочутливого неоднорщного елемента радюелектронно! системи у вигляд1 конструкцп «пластина-включення». Введена кусково-лшшна апроксимащя температури на поверхнях спряження пластини з включенням дозволила повшстю лшеаризувати в1дпов1дну частково лшеаризовану крайову задачу в1дносно лшеаризуючо! функцп. Шсля цього стало можливим застосувати штегральне перетворення Фур'е за одтею з просторових координат до отримано! лшшно! задач! та визначити лшеаризуючу функщю. Розглянуто лшшну залежшсть коефщента теплопровщност вщ температури для конструкцшних матер1атв з використанням лшеаризуючо! функцп. Шляхом розв'язання крайово! задач1 отримано формули для визначення температурного поля в термочутливш конструкцп «пластина-включення».

Результати. 1з використанням отриманих формул для визначення температурного поля у термочутливому неоднор1дно-му елеменп радюелектронно! системи створено обчислювальш програми, яю дають змогу отримати числов1 значення розподшу температури та анал1зувати температурш режими.

Висновки. Розроблена математична модель розрахунку температурного поля в термочутливш конструкцп «пластина-включення» е адекватною до реального ф1зичного процесу так як не спостерггаеться стрибка температури на поверхнях спряження пластини з включенням. Числов1 результати для вибраних матер1ал1в за лшшно! залежноси коефщента тепло-провщноста в1д температури в1др1зняються вщ результата, отриманих для сталого коефщента теплопровщносп, на 7%. Перспективи подальших дослщжень полягатимуть у розгляд1 складнших геометричних зображень елеменпв радюелек-тронних систем.

КЛЮЧОВ1 СЛОВА: температура, теплопровщтсть, нелшшна крайова задача, 1зотропна безмежна термочутлива пластина з тепло1зольованими лицевими поверхнями, чужорщне наскр1зне включення, щеальний тепловий контакт, тепловий потж.

УДК 536.24

ИССЛЕДОВАНИЕ ТЕМПЕРАТУРНЫХ РЕЖИМОВ В НЕОДНОРОДНЫХ ТЕРМОЧУВСТВИТЕЛЬНЫХ ЭЛЕМЕНТАХ РАДИОЭЛЕКТРОННЫХ УСТРОЙСТВ

Гаврыш В. И. - д-р техн. наук, профессор кафедры программного обеспечения, Национальный университет «Львовская Политехника», Львов, Украина.

Баранецкий Я. О. - канд. физ.-мат. наук, доцент кафедры высшей математики, Национальный университет «Львовская Политехника», Львов, Украина.

Коляса Л. И. - канд. физ.-мат. наук, старший преподаватель кафедры высшей математики, Национальный университет «Львовская Политехника», Львов, Украина.

АННОТАЦИЯ

Актуальность. Получено решение нелинейной краевой задачи теплопроводности для термочувствительного неоднородного элемента радиоэлектронной системы в виде полосы со сквозным включением, численно-аналитическое решение которой позволяет анализировать в нем температурные режимы. Цель работы - разработка метода линеаризации нелинейной математической модели теплопроводности, который позволяет получить численно-аналитическое решение соответствующей нелинейной краевой задачи для определения температурного поля в элементе радиоэлектронных устройств, которые геометрически изображены термочувствительной пластиной со сквозным включением.

Метод. Предложено линеаризирующую функцию, которая позволяет частично линеаризировать исходную нелинейную математическую модель теплопроводности для термочувствительного неоднородного элемента радиоэлектронной системы в виде конструкции «пластина-включение». Введеная кусочно-линейная аппроксимация температуры на поверхностях сопряжения пластины с включением позволила полностью линеаризировать соответствующую частично линеаризованную краевую задачу относительно линеаризирующей функции. После этого стало возможным применить интегральное преобразование Фурье по одной из пространственных координат к полученной линейной задаче и определить линеаризирующую функцию. Рассмотрено линейную зависимость коэффициента теплопроводности от температуры для конструкционных материалов с использованием линеаризирующей функции. Путем решения краевой задачи получены формулы для определения температурного поля в термочувствительной конструкции «пластина-включение».

Результаты. С использованием полученных формул для определения температурного поля в термочувствительном неоднородном элементе радиоэлектронной системы созданы вычислительные программы, которые позволяют получить числовые значения распределения температуры и анализировать температурные режимы.

Выводы. Разработана математическая модель расчета температурного поля в термочувствительной конструкции «пластина-включение» является адекватной к реальному физическому процессу так как не наблюдается скачка температуры на поверхностях сопряжения пластины с включением. Числовые результаты для выбранных материалов по линейной зависимости коэффициента теплопроводности от температуры отличаются от результатов, полученных для постоянного коэффициента теплопроводности, на 7%. Перспективы дальнейших исследований будут заключаться в рассмотрении сложных геометрических изображений элементов радиоэлектронных систем.

КЛЮЧЕВЫЕ СЛОВА: температура, теплопроводность, нелинейная краевая задача, изотропная термочувствительная пластина с теплоизолированными лицевыми поверхностями, инородное сквозное включения, идеальный тепловой контакт, тепловой поток.

Л1ТЕРАТУРА / ЛИТЕРАТУРА

1. Саваторова В. Л. Решение уравнения тепло проводнос-ти в неоднородной среде с учетом температурной зависимости коэффициента теплопроводности /

B. Л. Саваторова // Обозрение прикладной и промышленной математики. - 2009. - 17, № 1. -

C. 135-137.

2. Кудинов В. А. Анализ нелинейной теплопроводности на основе определения фронта температурного возмущения / В. А. Кудинов // Теплофизика высоких температур. - 2006. - 44, № 3. - С. 577-585.

3. Кудряшов Н. А. Приближенные решения одномерных задач нелинейной теплопроводности при заданном потоке / Н. А. Кудряшов // Журнал вычислительной математики и математической физики. - 2007. - 47, № 1. - С. 110-120.

4. Carpinteri A. Thermoelastic mismatch in nonhomogeneous beams / A. Carpinteri, M. Paggi // J. Eng. Math. - 2008. -61, No. 2-4. - P. 371-384.

5. Noda N. Thermal stresses in materials with temperature-dependent properties / N. Noda // Appl. Mech. Rev. -1991. - 44. - P. 383-397.

6. Noda N. Thermal stresses in materials with temperature-dependent properties / N. Noda // Appl. Mech. Rev. -1991. - 44. - P. 383-397.

7. Otao Y. Optimization of material composition of functionality graded plate for thermal stress relaxation using a genetic algorithm / Y. Otao, O. Tanigawa, O. Ishimaru // J. Therm. Stresses. - 2000. - 23. - P. 257-271.

8. Tanigawa Y. Transient heat conduction and thermal stress problems of a nonhomogeneous plate with temperature-dependent material properties / Y. Tanigawa, T. Akai, R. Kawamura // J. Therm. Stresses. - 1996. - 19, No. 1. -P. 77-102.

9. Tanigawa Y. Transient thermoelastic analysis of functionally graded plate with temperature-dependent material properties taking into account the thermal radiation / Y. Tanigawa, Y. Otao // Nihon Kikai Gakkai Nenji Taikai Koen Ronbunshu. - 2002. - 2. - P. 133-134.

10. Yangian Xu. Analysis of steady thermal stress in a ZrO2/FGM/Ti-6Al-4V composite ECBF plate with temperature-dependent material properties by NFEM / Xu. Yangian, Tu. Daihui //2009-WASE Int. Conf. on Informa. Eng. - Vol. 2-2. - P. 433-436.

11. Голицына Е. В. Квазистационарная трехмерная задача теплопроводности во вращающемся сплошном цилиндре из композиционного материала с нелинейными граничными условиями / Е. В. Голицына, Ю. Р. Осипов // Конструкции из композиционных материалов. -2007. - № 4. - C. 47-58.

12. Кудряшов Н. А. Приближенные решения первой и второй краевых задач нелинейной теплопроводности на полубесконечной прямой / Н. А. Кудряшов, М. А. Чмыхов // Инженерная физика. - 2007. - № 3. -С. 12-15.

13. Kudinov V. A. Analysis of nonlinear heat conduction based on determining the front of temperature perturbation / V. A. Kudinov, B. V. Averin, E. V. Stefanyuk, S. A. Nazarenko // High Temperature. - 2006. - 44, № 4. -P. 574-583.

14. Попович В. С. Нелшшна задача теплопровщносп для кулi з теплообмшом / В. С. Попович, К. С. 1ванкв //

Вгсник Львiв. ун-ту. Сер. Прикл. математика та шформатика. - 2002. - № 5. - С. 136-144.

15. Савула Я. Г. Дослвдження варiацiйноl задачi теплопровщносп у багатошарових середовищах з тонкими включеннями / Я. Г. Савула, Л. М. Дяконюк // Вюник ЛНУ iм. I. Франка, Сер. прикл. матем. та шформат. - 2000. - № 3. - С. 125-130.

16. Гавриш В. I. Моделювання температурних режимiв у кусково-однорщних структурах / В. I. Гавриш, Д. В. Федасюк. - Львiв : В-во Нац. ун-ту «Льв1вська полггехнжа», 2012. - 176 с.

17. Гаврыш В. И. Численно-аналитическое решение нелинейной стационарной задачи теплопроводности для бесконечной термочувствительной многослойной пластины / В. И. Гаврыш // Электронное моделирование. - 2014. - 36, № 3.- С. 59-70.

18. Гавриш В. I. Дослщження температурних режимiв у термочутливш пластиш з чужорiдним нас^зним включенням / В. I. Гавриш // Фiзико-математичне моделювання та шформацшш технологи. Науковий збiрник. - 2013. - Випуск 18. - С. 43-50.

19. Гавриш В. I. Нелшшна крайова задача теплопровщност для шарувато! пластиш з включенням / В. I. Гавриш // Фiзико-хiмiчна мехашка матерiалiв. -2015. - 51, № 3 - С. 32-38.

20. Bayat A. Thermo-mechanical analysis of functionally graded thick spheres with linearly time-dependent temperature [Text] / A. Bayat, H. Moosavi, Y. Bayat // Scientia Iranica. - 2015. - Vol. 22, Issue 5. - P. 18011812.

21. Mohazzab A. H. Two-Dimensional Stresses in a Hollow FG Sphere with Heat Source [Text] / A. H. Mohazzab, M. Jabbari // Advanced Materials Research. - 2011. -Vol. 264-265. - P. 700-705. DOI: 10.4028/www.scientific.net/amr.264-265.700

22. Ghannad M. A thermoelasticity solution for thick cylinders subjected to thermo-mechanical loads under various boundary conditions / M. Ghannad, M. P. Yaghoobi // International Journal of Advanced Design & Manufacturing Technology. - 2015. - Vol. 8, No. 4. -P. 1-12.

23. Jabbari M. Radially symmetric steady state thermal and mechanical stresses of a poro FGM hollow sphere / M. Jabbari, S. Karampour, M. R. Eslami // International Scholarly Research Network ISRN Mechanical Engineering. - 2011. - Vol. 2011. - P. 1-7. Article ID 305402, DOI: 10.5402/2011/305402

24. Подстригач Я. С. Термоупругость тел неоднородной структуры / Я. С. Подстригач, В. А. Ломакин, Ю. М. Коляно. - М. : Наука, 1984. - 368 с.

25. Коляно Ю. М. Методы теплопроводности и термоупругости неоднородного тела / Ю. М. Коляно. -К. : Наукова думка, 1992. - 280 с.

26. Корн Г. Справочник по математике для научных работников и инженеров / Г. Корн, Т. Корн. - М. : Наука, 1977. - 720 с.

27. Ломакин В. А. Теория упругости неоднородных тел / В. А. Ломакин. - М. : Изд-во Моск. ун-та, 1976. -376 с.

28. Берман Р. Теплопроводность твердых тел / Р. Берман. - М. : Мир, 1979. - 288 с.