No. 5 2010
ON THE TOLERANCE AVERAGING OF HEAT CONDUCTION FOR PERIODIC HEXAGONAL-TYPE
COMPOSITES
Ewaryst WIERZBICKI Warsaw University of Life Sciences Faculty of Applied Informatics and Mathematics
Departament of Applied Mathematics ul. Nowoursynowska 159, 02-776 Warsaw, Poland
ewaryst_wierzbicki@sggw.pl
The paper deals periodic composites which two-dimensional material structure is determined by the anisotropic conductivity matrix invariant with respect to the pair of translations such that the periodic cell determined by them coincide with the shapely hexagon. Moreover, it will be assumed that the material structure is invariant over the π/3-rotations with the centers of every hexagon as the origin of the rotation. The main aim of this paper is to prove that the tolerance model of such composite, obtained by the orthogonalization method, is isotropic.
Keywords: Tolerance averaging, periodicity cell, shapely hexagon, orthogonalization method
1. ITRODUCTION
Throughout this contribution we shall deal with the microheterogeneous periodic hexagonal-type composite solid made of the perfectly bonded constituents. The behaviour of these solids will be restricted to the heat conduction problem based on the Fourier heat conduction law and will be inveatigated in the framework of the well known two-dimensional parabolic equation which, under denotations
3 2
1 3
2
1, ,0], [0,0, ], [ , ,0],
[∂ ∂ ∂= ∂ x= x x z=x
=
∇ , can be written in the form
0 )
( )
( ∇ −∂⋅ ∂ + =
⋅
∇
−
∂ w w a w f
c t A (1.1)
Symbol w=w(⋅) stand here for a temperature field defined in Ω⊂R2, (⋅)
= f
f be the known density of heat sources. In the above equation c=c(⋅) is the heat flux and
=
) )
) ,
( ] 0 , 0 [
] 0 , 0 [ , (
( a z
z x x
z T
K A
is the conductivity martrix. We shall also assume the heat flux continuity conditions in normal directions on the interfaces Γ, Γ⊂Ω, between the constituents of the considered composite.
The well-known fact is that due to the discontinuous and highly oscillating form of functional coefficients c(·), A(·), the direct application of (1.1) to the analysis of special problems in most cases is difficult. That is why the mentioned above heat conduction problem is usually replaced, under some additional assumption, by some other problems with more regular coefficients. The characteristing feature of these assumptions is that microstructure of considered material is characterized by a certain scalar parameter λ>0. In this case conductivity matrix A in (1.1) depends on
λ
, K=Kλ. In many cases, in which the local averaged temperature u can be introduced, the tolerace averaging technique approach is very useful to the averaging of the heat conduction problem (1.1).The tolerance averaged model obtained on this way consists of the system of differential equations with more regular coefiicients for averaged temperature
) ( )
(⋅ ∈ 01 Ω
=u SV
u and fluctuation amplitudes wA =wA(⋅,t)∈SV01(Ω), which are new basic unknowns. Just introduced functional space SV01(Ω) is a certain new space consisting of slowly-varying functions. For particulars the reader is reffered to [1,2].
2. FORMULATION OF THE MODELLING EQUATIONS
There are known two methods via which the tolerace averaged model can be obtained. The first one is the method based of a new concept of extended stationary action principle introduced in [1]. This method have been resulted in many applications dealing functionally graded materials. To the periodic problems the orthogonalization metod explained in [2] usually have been applied. Since the subject of this paper is a certain periodic material structures we shall restrict ourselves to the tolerance model equations obtained by orthogonalization metod. To thos end following procedure explained in [2] we look for the temperature foield in the form
) , , ( ) ( ) , , ( ) , ,
(x z t u x z t g xW x z t
w = + A A (2.1)
Where u(x,t)=〈c〉−1〈cw〉(x,t) is the averaged temperature field and WA(x,t) are extra unknowns which are usually referred to as fluctuation amplitudes.
Here and in the sequel 〈⋅〉 stand for the averaged operator, [2], superscripts
denoted by capital latins A,B,... run over 1,2,...,N, where N is a number of fluctuation amplitudes. Shape functions gA(x), caused by the periodic structure of the composite, should be periodic and satisfy some additional conditions like
=0
〉
〈gA and gA∈O(λ), cf. [2]. Foolowing [2] the system of tolerance averaged equations will be rewritten in the form
−〈 〉
=
〉∇
⋅
〈∇ +
+
∇ 〉
⋅
⋅
〈∇ +
〉∂
−〈
〉∂
〈
〉
=〈
〉∇
∇
⋅ +〈
〉∂
−〈
〉∂
+〈
〉∇
〈
A A
B B A
B B A B t B A
A A t
fg u g
W g g
W g ag W
g cg
f W g u
c u a u
A
A A
A
2 2
2
(2.2)
Under the additional assumption that unknowns depends exclusively on x the above system ofequations reduces to the form
−〈 〉
=
〉∇
⋅
〈∇ +
∇ 〉
⋅
⋅
〈∇ +
〉∂
〈
〉
=〈
〉∇
∇
⋅ +〈
〉∂
−〈
〉∇
〈
A A
B B A
B t B A
A A t
fg u g
W g g
W g cg
f W g u
c u
A A
A A 2
(2.3)
investigated in [3]. Following [3] we shall restrict considerations to the case of honeycomb-type composites satisfying two following assumptions:
Assumption 1. The material structure of the anisotropic conductor is invariant under 2π /3 – rotations with respect to the center of a regular periodicity cell.
Assumption 2. The sequence g ,...,1 gN of the shape functions is invariant under 2π /3 – rotations with respect to the center of a regular periodicity cell.
We are to show that under the above assumptions tolerance equations (2.3) can be transformed to the equivalent system of equations with the isotopic coefficients.
3. ISOTROPIC PROPERTIES OF MODEL EQUATIONS
To transform tolerance equations (2.3) to the isotropic form we are to outline the line of approach similar to this presented in [3]. Firstly, we shall rewrite decomposition (2.1) in the form
) , , ( ) ( ) , , ( ) , ,
(x z t u x z t g xW x z t
w = + ra ra (3.1)
where indices a, r run over the sequences 1,2,,,,n and 1,2,3, respectively. In the decomposition mentioned above decomposition
) ( )
1(x g Qx
gar+ = ra (3.2)
where Q denote the matrix of rotation over π/3 and Wra are fluctuation amlitudes W also renumerated on the same way. Hence N=3n. A
Now, instead of fluctuation amplitud Wra we shall introduce new amplitudes
a r r a a a a
a W W W W
U = 1 + 2 + 3 , v =t (3.3)
where tr are three vectors determines the basic hexagon. Formula (1.7) represents an invertible linear transformation between W and ra Ua,va. Let us introduce the following averaged coefficients
, ,
) (
, ) (
) (
, ) )(
(
, ,
, ,
3 2 1
3 2 1 3
2 1 2
3 2 1 3 2 1 2
2 2
2
r a r a a a a a
b b b a
a a ab
b b b a a a ab
s r b s a r ab s r b s a r ab
r a r a
s r b s a
r ab
fg g
g g f f
g g g g
g g a
g g g g g g c c
g ag g
cg
g g
g
t f
A
t t K
t t C
t A
B t t A A
〉
=〈 + 〉
+
=〈
+ 〉 +
∇
⋅
⋅ + +
〈∇
=
+ 〉 + +
+
=〈
⊗
〉
=〈
⊗
〉
=〈
⊗
∇ 〉
⋅
=〈
⊗
∇ 〉
⋅
⋅
〈∇
=
(3.4)
where t2 =Qt1, t3 =Qt2, t1=Qt3, are three vectors determines the basic hexagon. Rather simple manipulations yield to
a b ab b ab
a T
a b ab b ab b ab
t a a
f U a U c
u f
u c u
= +
−
=
∇
⋅ +
⋅ +
∂
⋅
−
⋅
〉
=〈
〉∂
−〈
∇ +
∇
⋅
〉
⋅〈
∇
2 2
2 2
2
2 ( )
:
&
& K v A v B f
V C
v B A
(3.5)
It can be proved that
, ,
, ,
,
2 2 2 2
2 2
2 2 2
∈ +
=
∈ +
=
∈ +
=
∈ +
=
=
〉
〈
ab ab ab ab
ab ab
a a a ab ab ab
k k c
c
b b a
a k
( ( )
)
( ( )
)
1 K
1 C
1 B
1 A
1 A
(3.6)
for
, ,
,
~ , )
8 ( , 3
) (
~ , )
4 ( , 3 )
4 ( , 3
5 . 0
2 2
2 2
2 2
rs b s a r ab rs b s a r ab rs b s a r ab rs b s a r ab
s b s a
r r ab s b s a
r r ab
r a r a
r a r a
g ag k
g ag k
g cg c
g cg c
g tr g a
g tr g a
g tr b
g tr b
tr k
∈
〉
=〈
〉
=〈
∈
〉
=〈
〉
=〈
⋅
∇ 〉
〈∇
⋅
=
⋅
∇ 〉
〈∇
⋅
=
⊗
∇ 〉
= 〈
⊗
∇ 〉
= 〈
〉
= 〈
( ( )
)
( )
( )
δ δ
t A
t t
A t
t A
t A
A
(3.7)
Hence system of equations (3.5) is isotropic.
4. CONCLUDING REMARKS
Characteristic features of equations (3.5) can be listed as follows.
1o Model equations (3.5) are isotropic and hence we have proved the main thesis of this note that the averaged heat transfer response of the honeycomb- type rigid conductor is isotrpic even material properties of the conductor are anisotropic.
2o The obtained equations (3.5) include a special case in which every cell of the honeycomb-type considered conductor has a threefold axis of symmetry.
Similar result in elastodynamic was obtained In [4].
3o As it has been mentioned at the end of Section 2, in the special case in which basic unknowns do not depend on z results obtained In this note reduces to results obtained in [3].
REFERENCES
1. Thermomechanics of microheterogeneous solids and structures. Tolerance averaging approach, Eds. Woźniak Cz., Michalak B., Jędrysiak J., Wydawnictwo Politechniki Łódzkiej, Łódź 2009
2. Woźniak Cz., Wierzbicki E.: Averaging techniques in thermomechanics of composite solids, Wydawnictwa Politechniki Częstochowskiej, Częstochowa 2000
3. Wierzbicki E., Siedlecka U.: Isotropic models for a heat transfer in periodic composites, PAMM, 4, 1 (2004), 502-503
4. Nagórko W., Wągrowska M.: A contribution to modeling of composite solids, J. Theor. Appl. Mech., 1, 40 (2002), 149-158
UŚREDNIANIE TOLERANCYJNE PRZEWODNICTWA CIEPŁA W KOMPOZYTACH O PERIODYCZNEJ STRUKTURZE HEKSAGONALNEJ
S t r e s z c z e n i e
W pracy zaproponowano równoważną postać tolerancyjnie uśrednionego modelu przewodnictwa ciepła kompozytu o mikroperiodycznej strukturze materialnej typu plastra miodu. Założono niezmienniczość tej mikrostruktury względem obrotów o kąt 2π/3 a także taką samą niezmienniczość stosowanych funkcji kształtu. W otrzymanej nowej postaci model ten ma współczynniki izotropowe.