Next Article in Journal
Hybrid Model Development for HVAC System in Transportation
Next Article in Special Issue
A Model of Damage for Brittle and Ductile Adhesives in Glued Butt Joints
Previous Article in Journal
RoboEye, an Efficient, Reliable and Safe Semi-Autonomous Gaze Driven Wheelchair for Domestic Use
Previous Article in Special Issue
Effective Complex Properties for Three-Phase Elastic Fiber-Reinforced Composites with Different Unit Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interface Models in Coupled Thermoelasticity

1
Department of Civil and Building Engineering, and Architecture, Università Politecnica delle Marche, 60121 Ancona, Italy
2
IMAG CNRS UMR 5149, University of Nîmes, 30000 Nîmes, France
3
Department of Engineering, University of Ferrara, 44122 Ferrara, Italy
4
CNRS, Centrale Marseille, Laboratoire de Mécanique et d’Acoustique, Aix-Marseille University, 13453 Marseille, France
*
Author to whom correspondence should be addressed.
Submission received: 12 February 2021 / Revised: 1 March 2021 / Accepted: 2 March 2021 / Published: 4 March 2021
(This article belongs to the Special Issue Advances in Multiscale and Multifield Solid Material Interfaces)

Abstract

:
This work proposes new interface conditions between the layers of a three-dimensional composite structure in the framework of coupled thermoelasticity. More precisely, the mechanical behavior of two linear isotropic thermoelastic solids, bonded together by a thin layer, constituted of a linear isotropic thermoelastic material, is studied by means of an asymptotic analysis. After defining a small parameter ε , which tends to zero, associated with the thickness and constitutive coefficients of the intermediate layer, two different limit models and their associated limit problems, the so-called soft and hard thermoelastic interface models, are characterized. The asymptotic expansion method is reviewed by taking into account the effect of higher-order terms and defining a generalized thermoelastic interface law which comprises the above aforementioned models, as presented previously. A numerical example is presented to show the efficiency of the proposed methodology, based on a finite element approach developed previously.

1. Introduction

The use of composite structures, obtained by bonding together simpler structural members, has spread in all fields of engineering in the last decades. On the one hand, the structural assembly presents a significant improvement of the mechanical properties and an enhancement of its performances. On the other hand, the bonded joints among the composite constituents may cause a jump of the physical fields at the interface level and radically modify the global mechanical response. Thus, the correct modeling of composite interfaces is crucial in the understanding and design of complex structures.
From a theoretical point of view, the bonded region is considered as a thin interphase between two adjacent parts. By letting the thickness of this layer tend to zero, the interphase is reduced into a two-dimensional surface, called imperfect interface, where ad-hoc transmission conditions in terms of the representative physical fields are prescribed. The contact laws can be derived by means of classical variational tools and more refined mathematical techniques, in different physical frameworks, involving uncoupled (thermal conduction and elasticity) and coupled (piezoelectricity and multiphysics) phenomena.
Concerning the thermal (or electrical) conduction case, two main interface laws have been formulated: the lowly-conducting (LC) or Kapitza’s model and highly-conducting (HC) model. The LC model provides a discontinuity of the temperature field (electric potential) and a continuity of the normal heat flow (electric displacement) across the interface (see, e.g., [1,2,3]). The HC model gives rise to two-dimensional Young–Laplace equation, defined on the interface, depending on the jump of the normal heat flow (electric displacement) and maintaining the temperature (electric potential) continuous (see, e.g., [4,5]). A unifying approach of a general imperfect interface model, involving the concurrent jump of both the temperature field and the normal heat flow, recovering both the LC and HC models, was proposed by [6,7]. Concerning the linear elastic case, three types of imperfect interfaces have been proposed: the spring-layer interface model (SL) (soft interface), the coherent interface (CI) (rigid interface), and the general imperfect interface. The SL models considers that the traction vector is continuous across the interface, while the displacement presents a jump linearly proportional to the traction vector (see, e.g., [8,9]). The CI model has been developed for continuum theories with surface effects and nano-sized materials (see, e.g., [10,11,12]): the traction vector suffers a jump, while the displacement field is continuous across the interface. Finally, in the general imperfect model, both the displacement and normal traction fields are discontinuous across the interface [13,14].
The asymptotic expansions method and convergence approaches represent mathematical tools, usually employed in the derivation and justification of classical thin structures and layered plates [15,16,17,18,19]. These methodologies are based on the behavior of the problem solution, when a small parameter ε , related to the thickness of the interphase, tends to zero. Considering that the material properties of the intermediate layer depend on ε p , different limit behaviors can be derived by means of the asymptotic analysis: for p = 1 , an SL interface model can be recovered (see, e.g., [20,21]); and, for p = 1 , the CI interface model is mathematically justified by means of strong convergence arguments in [22,23]. Within the framework of a higher-order theory, assuming the interphase elastic constants are independent of the small thickness ( p = 0 ), the asymptotic analysis yields to a general stiff imperfect interface condition, prescribing both the jumps of the displacement and traction vector fields and recovering as a particular case the perfect contact conditions at the zeroth-order [24,25,26,27,28].The above transmission conditions have been generalized by considering some multiphysics and multifield couplings, such as in piezoelectricity and magneto-electro-thermo-elasticity [29,30], poroelasticity [31], and micropolar elasticity [32].
The goal of the present work is to identify the interface limit models of a composite constituted by a thin thermoelastic layer surrounded by two thermoelastic bodies in the framework of dynamic coupled thermoelasticity. Different situations are analyzed by varying the stiffnesses ratios between the middle layer and the adherents: namely, the soft thermoelastic lowly conducting interface, where the intermediate material coefficients have the order of magnitude ε with respect to those of the surrounding bodies, and the hard thermoelastic moderately conducting interface, where the constitutive parameters have the same order of magnitude. Following the ideas of [33], a generalized interface law is derived, comprising the aforementioned behaviors. A numerical investigation was performed in the framework of the finite element method (FEM), employing the approach developed in [34] for multiphysics problems, in order to assess the validity of the asymptotic models. Convergence results and a comparison between the full 3D model and the generalized interface problem are given.

2. Position of the Problem

In the sequel, Greek indices range in the set { 1 , 2 } , Latin indices range in the set { 1 , 2 , 3 } , and the Einstein’s summation convention with respect to the repeated indices is adopted. Let us consider a three-dimensional Euclidian space identified by R 3 and such that the three vectors e i form an orthonormal basis. Let M n be the space of n × n square matrices. We introduce the following notations for the inner products: a · b : = a i b i , for all vectors a = ( a i ) and b = ( b i ) in R 3 and A : B : = A i j B i j , for all A = ( A i j ) and B = ( B i j ) in M 3 .
Let us define a small parameter 0 < ε < 1 . We consider the assembly constituted of two solids Ω ± ε R 3 , called the adherents, bonded together by an intermediate thin layer B ε : = S × ( ε 2 , ε 2 ) of thickness ε , called the adhesive, with cross-section S R 2 . In the following, B ε and S are called interphase and interface, respectively. Let S ± ε be the plane interfaces between the interphase and the adherents and let Ω ε : = Ω + ε B ε Ω ε denote the composite system comprising the interphase and the adherents (cf. Figure 1a).
Let us assume that Ω ± ε and B ε are constituted by three homogeneous linear isotropic thermoelastic materials, whose constitutive laws are defined as follows:
σ i j ε = λ ε e p p ε δ i j + 2 μ ε e i j ε β ε θ ε δ i j , S ε = c v ε θ ε + β ε e p p ε , q i ε = k ε θ , i ε ,
where σ ε = ( σ i j ε ) is the Cauchy stress tensor, associated with the linearized strain tensor e ε = ( e i j ε ) : = 1 2 ( u i , j ε + u j , i ε ) , S ε represents the thermodynamic entropy and q ε = ( q i ε ) is the heat flow field. Constants λ ε , μ ε , β ε , c v ε , and k ε represent the Lamé’s constants, the thermal stress coefficient, the calorific capacity, and the thermal conductivity, respectively.
The thermoelastic state is defined by the couple s ε : = ( u ε , θ ε ) , where u ε and θ ε are the displacement field and variation of temperature, respectively. The thermoelastic composite is subject to body forces f ε = ( f i ε ) : Ω ± ε × ( 0 , T ) R 3 and heat source h ε : Ω ± ε × ( 0 , T ) R , applied on the top and bottom bodies, while all thermo-mechanical loadings and inertia forces are neglected in the intermediate layer B ε . The thermoelastic state s ε verifies the following coupled thermoelasticity system:
ρ ε u ¨ ε div ε σ ε = f ε in Ω ± ε × ( 0 , T ) , S ˙ ε + 1 T 0 div ε q ε = h ε in Ω ± ε × ( 0 , T ) , div ε σ ε = 0 in B ε × ( 0 , T ) , S ˙ ε + 1 T 0 div ε q ε = 0 in B ε × ( 0 , T ) ,
where f ˙ = t f denotes the time derivative of f and T 0 represents a reference temperature. The transmission conditions across the interfaces S + , ε and S , ε implies the continuity of the state s ε and of its normal dual counterpart with respect to S ± , ε , meaning that [ u ε ] = 0 , [ θ ε ] = 0 , σ ε e 3 = 0 , q ε · e 3 = 0 on S ± , ε × ( 0 , T ) , where [ f ] stands for the jump function evaluated at the interface S ± , ε . The boundary conditions are posed on Γ ε × ( 0 , T ) , with Γ ε : = Γ + , ε Γ , ε ; we recall that Γ ε = Γ g ε Γ u ε . For simplicity, we assume homogeneous boundary conditions on Γ u ε × ( 0 , T ) , concerning displacements and temperature, and non-homogeneous boundary conditions on Γ g ε × ( 0 , T ) , concerning surface forces g ε = ( g i ε ) and surface heat flow q ε . Hence, one has: σ ε n ε = g ε and q ε · n ε = q ε on Γ g ε × ( 0 , T ) , and u ε = 0 and θ ε = 0 on Γ u ε × ( 0 , T ) , where n ε = ( n i ε ) is the outer unit normal vector to Ω ε . The initial conditions are posed in Ω ε . Let θ i n ε , u i n ε , and u ˙ i n ε be, respectively, the variation of temperature, the displacement, and velocity fields at time t = 0 ; one has θ ε ( x ε , 0 ) = θ ε ( 0 ) = θ i n ε u ε ( x ε , 0 ) = u ε ( 0 ) = u i n ε and u ˙ ε ( x ε , 0 ) = u ˙ ε ( 0 ) = u ˙ i n ε in Ω ε .
Let us introduce the functional spaces V ( Ω ε ) : = { v ε H 1 ( Ω ε ) ; v ε = 0 on Γ u ε } and V ( Ω ε ) : = [ V ( Ω ε ) ] 3 . Given a certain state s ε : = ( u ε , θ ε ) V ( Ω ε ) : = V ( Ω ε ) × V ( Ω ε ) , for all test functions r ε = ( v ε , ξ ε ) V ( Ω ε ) and for any fixed t ( 0 , T ) , we introduce the following bilinear and linear forms:
A ε ( s ε , r ε ) : = Ω ε ρ ε u ¨ ε · v ε + σ ε : e ε ( v ε ) + S ˙ ε ξ ε 1 T 0 q ε · ε ξ ε d x ε ,
L ε ( r ε ) : = Ω ± ε f ε · v ε + h ε ξ ε d x ε + Γ g ε g ε · v ε + q ε ξ ε d Γ ε .
The variational form of the coupled thermoelastic system defined on the variable domain Ω ε reads as follows:
Find s ε ( t ) V ( Ω ε ) , t ( 0 , T ) , such that A ¯ ε ( s ε , r ε ) + A ¯ + ε ( s ε , r ε ) + A ^ ε ( s ε , r ε ) = L ε ( r ε ) ,
for all r ε V ( Ω ε ) , with initial condition θ i n ε , u i n ε , and u ˙ i n ε . The coupled hyperbolic–parabolic equations associated with variation problem (5) imply a degenerate system. Hence, the standard existence theorems are not applicable. For instance, in [35,36], by applying the pseudo-monotone theory, a weak solution is provided for a sufficiently small thermal stress coupling coefficient. In [37,38], a solution to the implicit evolution equation is derived after time-differentiation of the equilibrium equation provided by sufficiently smooth data of the problem. Under suitable regularity properties of the initial data, source and boundary values, and constitutive parameters, the well-posedness of thermo-electro-elastic evolution problem is extensively discussed in [39]: the proof of existence, uniqueness, and regularity of the solution has been obtained through the Faedo–Galerkin method. The existence and uniqueness theorems have also been extended to the thermo-electro-magneto-elastic case [40] and can be easily adapted to the present coupled thermoelastic problem.

Rescaling

To study the asymptotic behavior of the solution of problem (5) when ε tends to zero, we rewrite the problem on a fixed domain Ω independent of ε . By using the approach of [15], we consider the bijection π ε : x Ω ¯ x ε Ω ¯ ε given by
π ε : π ¯ ε ( x 1 , x 2 , x 3 ) = ( x 1 , x 2 , x 3 1 2 ( 1 ε ) ) , for all x Ω ¯ ± , π ^ ε ( x 1 , x 2 , x 3 ) = ( x 1 , x 2 , ε x 3 ) , for all x B ¯ ,
where, after the change of variables, the adherents occupy Ω ± : = Ω ± ε ± 1 2 ( 1 ε ) e 3 and the interphase B = { x R 3 : ( x 1 , x 2 ) S , | x 3 | < 1 2 } . The sets S ± = { x R 3 : ( x 1 , x 2 ) S , x 3 = ± 1 2 } denote the interfaces between B and Ω ± and Ω = Ω + Ω B is the rescaled configuration of the composite. Lastly, Γ u and Γ g indicate the images through π ε of Γ u ε and Γ g ε (cf. Figure 1b). Consequently, x α ε = x α and x 3 ε = x 3 in Ω ± , and x α ε = x α and x 3 ε = 1 ε x 3 in B . In the sequel, only if necessary, s ¯ ε = ( u ¯ ε , θ ¯ ε ) and s ^ ε = ( u ^ ε , θ ^ ε ) denote the restrictions of functions s ε = ( u ε , θ ε ) to Ω ± and B.
The constitutive coefficients of Ω ± ε are assumed to be independent of ε , while the constitutive coefficients of B ε present the following dependences on ε : λ ^ ε = ε p λ ^ , μ ^ ε = ε p μ ^ , β ^ ε = ε p β ^ , c ^ v ε = ε p c ^ v , and k ^ m , ε = ε p k ^ , with p { 0 , 1 } . Two different limit behaviors are characterized according to the choice of the exponent p: by choosing p = 1 , a model for a soft thermoelastic interface with low conductivity is deduced; and, when p = 0 , a model for a hard thermoelastic interface with moderate conductivity is obtained. Finally, the data, unknowns, and test functions verify the following scaling assumptions: s ε ( x ε ) = s ε ( x ) , r ε ( x ε ) = r ( x ) x Ω , f ε ( x ε ) = f ( x ) , h ε ( x ε ) = h ( x ) x Ω ± , g ε ( x ε ) = g ( x ) , q ε ( x ε ) = q ( x ) , x Γ g . Thus, L ε ( r ε ) = L ( r ) .
According to the previous hypothesis, problem (5) can be reformulated on a fixed domain Ω independent of ε . Thus, the following rescaled problem (in the sequel, we omit the explicit dependences on time t of the unknowns and data) is obtained:
Find s ε V ( Ω ) , t ( 0 , T ) , such that A ¯ ( s ε , r ) + A ¯ + ( s ε , r ) + ε p + 1 A ^ ( s ε , r ) = L ( r ) ,
for all r V ( Ω ) , p { 0 , 1 } , with initial condition θ i n , u i n , and u ˙ i n , where
A ¯ ± ( s ε , r ) : = Ω ± ρ ε u ¨ ε · v ε + σ ε : e ( v ) + S ˙ ε ξ 1 T 0 q ε · ξ d x ,
A ^ ( s ε , r ) : = 1 ε 2 a 0 ( s ε , r ) + 1 ε a 1 ( s ε , r ) + a 2 ( s ε , r ) ,
where
a 0 ( s ε , r ) : = B K ^ u , 3 ε · v , 3 + k ^ T 0 θ , 3 ε ξ , 3 d x ,
a 1 ( s ε , r ) : = B ( K ^ α ) T u , α ε · v , 3 + K ^ α u , 3 ε · v , α β ^ θ ε v 3 , 3 + β ^ u ˙ 3 , 3 ε ξ d x ,
a 2 ( s ε , r ) : = B K ^ α β u , β ε · v , α β ^ ε θ ε v τ , τ + k ^ T 0 θ , α ε ξ , α + ( c ^ v θ ˙ ε + β ^ u ˙ α , α ε ) ξ d x
and
K ^ : = μ ^ 0 0 0 μ ^ 0 0 0 2 μ ^ + λ ^ , K ^ 1 : = 0 0 λ ^ 0 0 0 μ ^ 0 0 , K ^ 2 : = 0 0 0 0 0 λ ^ 0 μ ^ 0 ,
K ^ 11 : = 2 μ ^ + λ ^ 0 0 0 μ ^ 0 0 0 μ ^ , K ^ 22 : = μ ^ 0 0 0 2 μ ^ + λ ^ 0 0 0 μ ^ ,
K ^ 12 : = 0 λ ^ 0 μ ^ 0 0 0 0 0 , K ^ 21 = ( K ^ 12 ) T .
Now, an asymptotic analysis of the rescaled problem (7) can be performed. Since the rescaled problem (7) has a polynomial structure with respect to the small parameter ε , we can look for the solution s ε of the problem as a series of powers of ε :
s ε = s 0 + ε s 1 + ε 2 s 2 + , s ¯ ε = s ¯ 0 + ε s ¯ 1 + ε 2 s ¯ 2 + , s ^ ε = s ^ 0 + ε s ^ 1 + ε 2 s ^ 2 + .
where s ¯ ε = s ε π ¯ ε and s ^ ε = s ε π ^ ε . By substituting (16) into the rescaled problem (7), and by identifying the terms with identical power of ε , as customary, a set of variational problems is obtained to be solved in order to characterize the limit thermoelastic state s 0 , the first-order corrector term s 1 and their associated limit problem, for p { 0 , 1 } . The order 1 can be considered as a corrector term of the order 0, giving a better approximation of the initial model.

3. The Soft Thermoelastic Interface Model

In this section, the limit model for a soft thermoelastic interface model, corresponding to an adhesive which is weaker with respect to the adherents, is derived. By choosing p = 1 and injecting (16) into (7), the following set of variational problems P q is obtained:
P 0 : A ¯ ( s 0 , r ) + A ¯ + ( s 0 , r ) + a 0 ( s 0 , r ) = L ( r ) , P 1 : A ¯ ( s 1 , r ) + A ¯ + ( s 1 , r ) + a 0 ( s 1 , r ) + a 1 ( s 0 , r ) = 0 , P q : A ¯ ( s q , r ) + A ¯ + ( s q , r ) + a 0 ( s q , r ) + a 1 ( s q 1 , r ) + a 2 ( s q 2 , r ) = 0 , q 2
In the sequel, the limit problems at order 0 and order 1 are presented, by skipping all the mathematical technicalities involved in the solution of problems P q (see [33] for a detailed description of the asymptotic analysis).
  • Order 0 model
    Governing equations Transmission conditions on S ± ρ u ¯ ¨ 0 div σ ¯ 0 = f in Ω ± , S ¯ ˙ 0 + 1 T 0 div q ¯ 0 = h in Ω ± , σ ¯ 0 n = g on Γ g q ¯ 0 · n = q on Γ g , s 0 = 0 on Γ u , [ u ¯ 0 ] = K ^ 1 σ ¯ 0 e 3 , [ θ ¯ 0 ] = T 0 k ^ q ¯ 0 · e 3 , [ σ ¯ 0 e 3 ] = 0 , [ q ¯ 0 · e 3 ] = 0 .
  • Order 1 model
    Governing equations Transmission conditions on S ± ρ u ¯ ¨ 1 div σ ¯ 1 = 0 in Ω ± , S ¯ ˙ 1 + 1 T 0 div q ¯ 1 = 0 in Ω ± , σ ¯ 1 n = 0 on Γ g q ¯ 1 · n = 0 on Γ g , s 1 = 0 on Γ u , [ u ¯ 1 ] = K ^ 1 σ ¯ 1 e 3 ( K ^ α ) T u ¯ 0 , α + β ^ θ ¯ 0 e 3 , [ θ ¯ 1 ] = T 0 k ^ q ¯ 1 · e 3 , [ σ ¯ 1 e 3 ] = K α [ u ¯ 0 ] , α , [ q ¯ 1 · e 3 ] = β ^ [ u ¯ ˙ 3 0 ] ,
where f : = 1 2 ( f ( x ˜ , + ( 1 / 2 ) + ) + f ( x ˜ , ( 1 / 2 ) ) and f : = f ( x ˜ , + ( 1 / 2 ) + ) f ( x ˜ , ( 1 / 2 ) ) , x ˜ : = ( x α ) S denote, respectively, the mean value and the jump functions at the interfaces. The soft thermoelastic interface models at order 0 and order 1 present various similarities, compared with the linear elastic case [27]. At order 0, from a mechanical point of view, the interface behaves as linear springs reacting to the jump between the top and bottom displacements and temperature, while the traction vector and normal heat flow are remains continuous. The order 1 model provides a mixed contact law, expressed by a concurrent discontinuity in terms of thermoelastic state, traction vector, and normal heat flow. The order 1 transmission conditions can be also rewritten in terms of q ¯ 1 · e 3 and σ ¯ 1 e 3 , as follows:
σ ¯ 1 e 3 = K ^ [ u ¯ 1 ] + ( K ^ α ) T u ¯ 0 , α β ^ θ ¯ 0 e 3 , q ¯ 1 · e 3 = k ^ T 0 [ θ ¯ 1 ] .
The jump and mean values of the traction vector and normal heat flow at the interface depend on s 0 and are analogous to those obtained for the soft elastic case in [26]. It is interesting to notice that, at order 1, the jump of the heat flow at the interface inside the intermediate layer depend on the variation in time of the normal displacement u 3 .

4. The Hard Thermoelastic Interface Model

In this section, the limit model for a hard thermoelastic interface, corresponding to an intermediate layer having the same rigidities of the top and bottom bodies, is derived. Let p = 0 , the asymptotic expansion (16) is inserted in (7), and the following set of variational problems P q is obtained:
P 1 : a 0 ( s 0 , r ) = 0 , P 0 : A ¯ ( s 0 , r ) + A ¯ + ( s 0 , r ) + a 0 ( s 1 , r ) + a 1 ( s 0 , r ) = L ( r ) , P 1 : A ¯ ( s 1 , r ) + A ¯ + ( s 1 , r ) + a 0 ( s 2 , r ) + a 1 ( s 1 , r ) + a 2 ( s 0 , r ) = 0 , P q : A ¯ ( s q , r ) + A ¯ + ( s q , r ) + a 0 ( s q + 1 , r ) + a 1 ( s q , r ) + a 2 ( s q 1 , r ) = 0 , q 2
A detailed equivalent analysis on the solution of the variational problems P q can be found in [33]. In the sequel, the limit problems at order 0 and order 1 are presented.
  • Order 0 model
    Governing equations Transmission conditions on S ± ρ u ¯ ¨ 0 div σ ¯ 0 = f in Ω ± , S ¯ ˙ 0 + 1 T 0 div q ¯ 0 = h in Ω ± , σ ¯ 0 n = g on Γ g , q ¯ 0 · n = q on Γ g , s 0 = 0 on Γ u , [ u ¯ 0 ] = 0 , [ θ ¯ 0 ] = 0 , [ σ ¯ 0 e 3 ] = 0 , [ q ¯ 0 · e 3 ] = 0 .
  • Order 1 model
    Governing equations Transmission conditions on S ± ρ u ¯ ¨ 1 div σ ¯ 1 = 0 in Ω ± , S ¯ ˙ 1 + 1 T 0 div q ¯ 1 = 0 in Ω ± , σ ¯ 1 n = 0 on Γ g , q ¯ 1 · n = 0 on Γ g , s 1 = 0 on Γ u , [ u ¯ 1 ] = K ^ 1 σ ¯ 0 e 3 ( K ^ α ) T u ¯ , α 0 + β ^ θ ¯ 0 e 3 , [ θ ¯ 1 ] = T 0 k ^ q ¯ 0 · e 3 , [ σ ¯ 1 e 3 ] = K ^ α K ^ 1 σ ¯ 0 e 3 , α + L ^ α β u ¯ , α β 0 + + β ^ ( K ^ α K ^ 1 θ ¯ , α 0 e 3 θ ¯ , τ 0 e τ ) , [ q ¯ 1 · e 3 ] = β ^ λ ^ + 2 μ ^ σ ¯ ˙ 33 0 + Σ ˜ ˙ 0 k ^ T 0 Δ s θ ¯ 0 ,
where L ^ α β : = K ^ α β K ^ β K ^ 1 ( K ^ α ) T , Σ ˜ 0 : = β ˜ u ¯ α , α 0 + c ˜ v θ ¯ 0 , with β ˜ : = 2 μ ^ β ^ λ ^ + 2 μ ^ and c ˜ v : = c ^ v + β ^ 2 λ ^ + 2 μ ^ , Δ s denotes the two-dimensional Laplacian operator. Note that, in this case, θ ¯ 0 = θ ¯ 0 and u ¯ 0 = u ¯ 0 .
The hard thermoelastic interface problems above present the same structures of the analogous linear elastic hard interface models [24,25,26]. Concerning the order 0, the transmission conditions provide a continuity of the thermoelastic state and of its conjugated counterpart, which is typical for adhesives having the same rigidity properties of the adherents. In this case, the upper and lower bodies are perfectly bonded together. At order 1, a mixed interface model is obtained, characterized by a jump of the state and traction vector depending on the values of the thermoelastic state and traction vector at order 0. These order 0 terms are known since they have been determined in the previous problem and they appear in the formulation as source terms. The interface conditions at order 1 can be interpreted as the two-dimensional coupled thermoelastic problem defined on the plane of the interface.

5. Generalized Interface Transmission Conditions

In [26,33], it has been shown that it is possible to obtain a condensed form of transmission conditions summarizing both the orders 0 and 1 of the soft and hard cases in only one couple of equations in terms of the jump of the displacement field and tractions at the interface. Equivalently, it is possible to define an implicit general thermoelastic interface law starting from the hard case, comprising the order 0 and order 1 soft and hard thermoelastic interface models.
To this end, by denoting by s ˜ ε : = s ¯ 0 + ε s ¯ 1 , σ ˜ ε : = σ ¯ 0 + ε σ ¯ 1 and q ˜ ε : = q ¯ 0 + ε q ¯ 1 , suitable approximations of s ε , σ ¯ ε and q ¯ ε , respectively, and following the approach developed in [26,33], one can obtain the implicit form of the transmission conditions:
σ ˜ ε e 3 = 1 ε K ^ [ u ˜ ε ] + ( K ^ α ) T u ˜ ε , α β ^ θ ˜ ε e 3 , q ˜ ε · e 3 = k ^ ε T 0 [ θ ˜ ε ] , σ ˜ ε e 3 = K ^ α [ u ˜ ε ] , α ε K ^ α β u ˜ ε , α β + ε β ^ θ ˜ ε , α e α , q ˜ ε · e 3 = β ^ [ u ˜ ˙ 3 ε ] + ε Σ ˙ ε ε k ^ T 0 Δ s θ ˜ ε ,
with Σ ε : = c ^ v θ ˜ ε + β ^ u ˜ α ε , α .
To write the variational formulation of the general coupled thermoelastic interface problem, the expression of the general transmission conditions presented in (24) is employed. In what follows, for the sake of simplicity, the indices ε and symbol ( · ) ˜ are omitted. Let us write the variational form of the equilibrium equations on each sub-domain Ω + and Ω . The sum of the two equations leads to
Ω ± ρ u ¨ · v + σ : e ( v ) + S ˙ ξ 1 T 0 q · ξ d x S σ ( x ˜ , 0 + ) n ( x ˜ , 0 + ) · v + σ ( x ˜ , 0 ) n ( x ˜ , 0 ) · v d Γ + + S q ( x ˜ , 0 + ) · n ( x ˜ , 0 + ) ξ + q ( x ˜ , 0 ) · n ( x ˜ , 0 ) ξ d Γ = L ( r ) ,
which can be written
Ω ± ρ u ¨ · v + σ : e ( v ) + S ˙ ξ 1 T 0 q · ξ d x + S [ σ e 3 · v ] [ q · e 3 ξ ] d x ˜ = L ( r ) ,
letting e 3 = n ( x ˜ , 0 ) = n ( x ˜ , 0 + ) and d Γ = d x ˜ . Then, using the property [ a b ] = a [ b ] + [ a ] b and relations (24), and after an integration by parts, one has
Find s W ( Ω ˜ ) , Ω ˜ : = Ω + S Ω , t ( 0 , T ) , such that A ¯ ( s , r ) + A ¯ + ( s , r ) + A ( s , r ) = L ( r ) ,
for all r W ( Ω ˜ ) , where W ( Ω ˜ ) : = W ( Ω ˜ ) × W ( Ω ˜ ) , with W ( Ω ˜ ) : = { r H 1 ( Ω ˜ ) , r | S H 1 ( S ) , r = 0 on Γ u } , W ( Ω ˜ ) : = [ W ( Ω ˜ ) ] 3 , and
A ( s , r ) : = S 1 ε K ^ [ u ] · [ v ] + ( K ^ α ) T u , α · [ v ] + K ^ α [ u ] · v , α + ε K ^ α β u , β · v , α β ^ θ [ v 3 ] ε β ^ θ v α , α + 1 ε k ^ T 0 [ θ ] [ ξ ] + ε k ^ T 0 θ , α ξ , α + β ^ [ u ˙ 3 ] ξ + ε Σ ˙ ξ d x ˜ ,
L ( r ) : = Ω ± f · v + h ξ d x + Γ g g · v + q ξ d Γ + S F · v + H ξ d γ ,
where Σ : = c ^ v θ + β ^ u α , α , K ^ α [ u ] + ε K ^ α β u , β ε β ^ θ e α ν α : = F and ε k ^ T 0 θ , α ν α : = H denote the loads on the lateral boundary of the interface S , with outer unit normal vector ( ν α ) (see [33]).

6. FEM Implementation

The numerical simulations were carried out by means of the finite element method, discretizing the variational problem (27). This helped validate the proposed asymptotic approach. The FEM analysis was performed considering the coupled dynamic thermoelastic problem and comparing the solution of the three-phase model (two adherents and adhesive) with the generalized interface (two adherents + interface). The problem was solved employing the software GetFem++ (see [41,42] for more details), with a standard linear solver (conjugate gradient). For that purpose, standard piecewise linear finite elements were considered.
Let us consider a thermoelastic laminated plate occupying a 3D domain defined by Ω ε = [ 0 , L 1 ] × [ 0 , L 2 ] × [ 0 , 2 h + ε ] , with h = 1 cm, L 1 / h = 10 , L 2 / h = 5 . (see Figure 2). Clearly, with self-explanatory notation, x 1 = x , x 2 = y , and x 3 = z . The adherents are made of Material 1, while the adhesive is constituted by Material 2.
Simply supported boundary conditions are considered on the bottom edges of the composite plate. The plate is subject to a thermal shock q · >n | Γ t o p = q ( t ) = a t e b t on the top face, with a = 30 and b = 0.8 , whose plot is illustrated in Figure 3. The bottom face is thermally insulated q · n | Γ b o t t o m = 0 . No volume or surface mechanical loads were applied f = g = 0 .
The FEM discretization was carried out using piecewise linear finite elements on hexahedrons, with 7280 nodes (29,203 degrees of freedom) for the three-phase problem and 5824 nodes (23,635 degrees of freedom) for the problem with the generalized interface law. The time discretization was realized using a Newmark-beta scheme with β = 0.25 and γ = 0.5 for the second time derivative of the displacements, and a Crank–Nicholson scheme for the first derivative of the thermodynamic entropy. It is worth noting that both methods are unconditionally stable and of order two in time. The time step of discretization is equal to δ t = 0.1 s.
The numerical example considers a composite plate, in which the adherents and the adhesive have very different thermo-mechanical properties. Material 1 is aluminum (Al), while Material 2 is a polyvinyl chloride (PVC) foam. The constitutive parameters are listed in Table 1.
To evaluate the accuracy of the asymptotic analysis, the influence of the relative thickness ε / h , for fixed time instants, on the L 2 -relative error was investigated. The L 2 -relative errors u ε u L 2 u L 2 and θ ε θ L 2 θ L 2 was computed taking into account the solution ( u ε , θ ε ) of the initial three-phase problem, discretized with a FE mesh, and the solution ( u , θ ) of the interface problem (27). Table 2 and Table 3 report the relative error values for increasing time and vanishing relative thickness.
The convergence diagrams of the the relative L 2 -norms of the displacements and temperatures, obtained with the three-phase problem and the reduced interface problem, are plotted in Figure 4, as the thickness ratio ε / h tends to zero, at time t = 1 s and t = 10 s, respectively. Moreover, the evolution in time of the L 2 -relative is reported in Figure 5.
From the results in Table 2 and Table 3 and, especially, Figure 4, it can be noticed that, by decreasing the thickness ratio ε / h , the relative errors present an immediate reduction for fixed times. The convergence rate is of the order ( ε / h ) 2 and remains constant for increasing time instants. As illustrated in Figure 5, the evolution in time of the L 2 -relative error, for fixed ε / h , becomes approximately steady after t = 10 s for the displacement field, while it presents a decreasing trend concerning the temperature field. Besides, even for a relative thickness ε / h = 0.1 , at time t = 10 s, the relative error is close to about 1.89 × 10 3 , for the displacement field, and about 5.16 × 10 4 , for the temperature field. Hence, the proposed general thermoelastic interface model provides an acceptable solution and it is able to correctly approximate the solution ( u ε , θ ε ) of the three-phase problem. Moreover, the reduced model can also be employed for moderately thick adhesives.
In the sequel, the numerical results obtained by solving the general interface model are presented, considering a relative thickness of ε / h = 0.01 . Following the approach by [43], hereinafter, the results are provided using dimensionless variables:
  • U ( X 1 , X 2 , X 3 , t ) : = 1 ν 1 ( 1 + ν 1 ) α 1 T 0 u ( x 1 , x 2 , x 3 , t ) ,
  • Θ ( X 1 , X 2 , X 3 , t ) : = θ ( x 1 , x 2 , x 3 , t ) T 0 T 0 ,
  • Σ i j ( X 1 , X 2 , X 3 , t ) : = 1 ρ 1 V σ i j ( x 1 , x 2 , x 3 , t )
where X i = x i / , t = V , and and V are defined by
V = e 1 ( 1 ν 1 ) ( 1 + ν 1 ) ( 1 2 ν 1 ) ρ 1 , = k 1 ρ 1 c ν 1 V .
Let us notice that the domain Ω ε is chosen such that X 1 [ 0 ; 10 ] , X 2 [ 0 ; 10 ] and X 3 [ 0 ; 1 ] .
Figure 6 represents the trend of the displacement U 3 and temperature Θ , evaluated along X 3 on the orthogonal fiber to the mid-plane of the interface at point ( X ¯ 1 = 6 , X ¯ 2 = 6 ) , for given times. The plot shows that, after the thermal shock, the displacement U 3 evolves in opposite directions within the adherents: the composite laminated plate tends to expand and contract itself along the through-the-thickness axis. On the other hand, the temperature field Θ remains constant along the X 3 -axis within the adherents, for given times, reaching a steady value after a certain time interval. As expected, the plots also report a jump of the state fields ( U 3 , Θ ) in correspondence of the intermediate layer, and, thus, the adhesive behaves as a soft thermoelastic interface. This is mainly due to the material properties of the adhesive, which are smaller with respect to the those of the adherents.
Figure 7 illustrates the evolution of the displacement field U and temperature Θ with respect to the time t, at a given point X ¯ = ( 6.5 , 3.5 , 0 ) , placed on the bottom face of the composite plate. As expected, the thermal shock induces an oscillatory trend concerning the displacements. Conversely, the temperature evolves to a steady state, corresponding to a constant value, after a sudden increase related to the thermal shock application.
Figure 8 represents the trend of the stresses Σ 33 and Σ 13 , evaluated along X 3 on the orthogonal fiber to the mid-plane of the interface at point ( X ¯ 1 = 6 , X ¯ 2 = 6 ) , for given times. The plot shows that, after the thermal shock, the stress Σ 33 remains constant along the X 3 -axis within the adherents. In this particular case, the thermal contribution to Σ 33 is predominant with respect to the elastic one, i.e., Σ 33 B ( Θ + 1 ) : indeed, their diagrams present analogous trends and differ for a constant of proportionality B : = β T 0 ρ 1 V (see Figure 6 and Figure 8). The stress Σ 13 presents an oscillating behavior along X 3 inside the adherents, but its contribution is negligible compared with Σ 33 . Moreover, the normal Σ 33 and shear stresses Σ 13 , evaluated at the top and bottom faces of the intermediate layer, are very similar and, thus, their jump almost vanishes. This is typical of soft interface models, in which the thermoelastic state presents a discontinuity, while its conjugated quantities (traction vector and normal heat flow) are continuous across the interface.
Figure 9 shows a comparison between the evolution in time of U 3 ( X ¯ , t ) , at a given point X ¯ = ( 6.5 , 3.5 , 0 . ) , of a homogeneous three-layer plate, made of aluminum, and the Al/PVC composite plate. Although the thickness ratio is small ( ε / h = 0.01 ), the effect of the adhesive becomes relevant concerning the response of the plate to the thermal shock. The homogeneous plate appears to be stiffer with respect to the composite one, which manifests a significant amplitude and period increase of the U 3 motion.

7. Concluding Remarks

General imperfect interface conditions are proposed in the framework of coupled thermoelasticity, simulating the thermomechanical behavior of a thin-bonded joint. The approach used to obtain the transmission conditions is based on the asymptotic expansions method. Zero- and higher-order interface models are derived for soft and hard interphases. Following [33], a general transmission law, comprising the two regimes (soft and hard) at the various order, is derived. To assess the validity of the previous asymptotic approach, numerical simulations were developed using a finite element method, which generalizes an analogous methodology to dynamical coupled thermoelasticity, already proposed in [34] in the framework of piezoelectricity. The numerical example consisted of a thermoelastic composite three-layer aluminum plate, with a PVC adhesive, subject to a thermal shock. Two different configurations were considered: the first one consisted of an initial three-phase problem, while the second one took into account the FE discretized form of interface problem (27). The most significant fields (displacement and temperature) and their L 2 -relative errors were then computed and compared to test the validity of the proposed interface laws and the accuracy of the asymptotic model. The proposed general thermoelastic interface model provides an acceptable solution and it is able to correctly approximate the solution of the three-phase problem. These findings clearly indicate that the approach of substituting the interphase with the proposed interface law provides a robust modeling for the composite.

Author Contributions

The following statements should be used “Conceptualization, M.S., S.D., R.R. and F.L.; methodology, M.S., S.D., R.R. and F.L.; software, S.D.; validation, M.S. and S.D.; formal analysis, M.S., S.D. and R.R.; investigation, M.S.; data curation, M.S., S.D. and R.R.; writing—original draft preparation, M.S.; writing—review and editing, M.S., S.D., R.R. and F.L.; funding acquisition, R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Benveniste, Y. The effective conductivity of composites with imperfect thermal contact at constituent interfaces. Int. J. Eng. Sci 1986, 24, 1537–1552. [Google Scholar] [CrossRef]
  2. Benveniste, Y. Effective thermal-conductivity of composites with a thermal contact resistance between the constituents-nondilute case. J. Appl. Phys. 1987, 61, 2840–2843. [Google Scholar] [CrossRef]
  3. Kapitza, P.L. Collected Papers of P. L. Kapitza; ter Haar, D., Ed.; reprinted 1965; Pergamon: Oxford, UK, 1941. [Google Scholar]
  4. Pham Huy, H.; Sanchez-Palencia, E. Phenomenes de transmission à travers des couches minces de conductivité elevée. J. Math. Anal. Appl. 1974, 47, 284–309. [Google Scholar]
  5. Miloh, T.; Benveniste, Y. On the effective conductivity of composites with ellipsoidal inhomogeneities and highly conducting interfaces. Proc. R. Soc. Lond. A 1999, 455, 2687–2706. [Google Scholar] [CrossRef]
  6. Hashin, Z. Thin interphase/imperfect interface in conduction. J. App. Phys. 2001, 89, 2261–2267. [Google Scholar] [CrossRef]
  7. Javili, A.; Kaessmair, S.; Steinmann, P. General imperfect interfaces. Comput. Methods Appl. Mech. Engrg. 2014, 275, 76–97. [Google Scholar] [CrossRef]
  8. Benveniste, Y. Effective mechanical behaviour of composite materials with imperfect contact between the constituents. Mech. Mat. 1985, 4, 197–208. [Google Scholar] [CrossRef]
  9. Hashin, Z. Thin interphase/imperfect interface in elasticity with application to coated fiber composites. J. Mech. Phys. Solids 2002, 50, 2509–2537. [Google Scholar] [CrossRef]
  10. Gurtin, M.E.; Murdoch, A.I. A continuum theory of elastic material surfaces. Arch. Ration. Mech. Anal. 1975, 57, 291–323. [Google Scholar] [CrossRef]
  11. Yvonnet, J.; Le Quang, H.; He, Q.C. An XFEM/level set approach to modelling surface/interface effects and to computing the size-dependent effective properties of nanocomposites. Comput. Mech. 2008, 42, 119–131. [Google Scholar] [CrossRef] [Green Version]
  12. Yvonnet, J.; Mitrushchenkov, A.; Chambaud, G.; He, Q.C. Finite element model of ionic nanowires with size-dependent mechanical properties determined by ab initio calculations. Comput. Methods Appl. Mech. Eng. 2011, 200, 614–625. [Google Scholar] [CrossRef] [Green Version]
  13. Benveniste, Y.; Miloh, T. Imperfect soft and stiff interfaces in two-dimensional elasticity. Mech. Mat. 2001, 33, 309–323. [Google Scholar] [CrossRef]
  14. Benveniste, Y. A general interface model for a three-dimensional curved thin anisotropic interphase between two anisotropic media. J. Mech. Phys. Solids 2006, 54, 708–734. [Google Scholar] [CrossRef]
  15. Ciarlet, P.G. Mathematical Elasticity, Vol. II: Theory of Plates; North-Holland: Amsterdam, The Netherlands, 1997. [Google Scholar]
  16. Serpilli, M.; Lenci, S. Asymptotic modelling of the linear dynamics of laminated beams. Int. J. Solids Struct. 2012, 49, 1147–1157. [Google Scholar] [CrossRef] [Green Version]
  17. Serpilli, M.; Lenci, S. An overview of different asymptotic models for anisotropic three-layer plates with soft adhesive. Int. J. Solids Struct. 2016, 81, 130–140. [Google Scholar] [CrossRef]
  18. Furtsev, A.; Itou, H.; Rudoy, E. Modeling of bonded elastic structures by a variational method: Theoretical analysis and numerical simulation. Int. J. Solids Struct. 2020, 182–183, 100–111. [Google Scholar] [CrossRef]
  19. Furtsev, A.; Rudoy, E. Variational approach to modeling soft and stiff interfaces in the Kirchhoff-Love theory of plates. Int. J. Solids Struct. 2020, 202, 562–574. [Google Scholar] [CrossRef]
  20. Klarbring, A. Derivation of the adhesively bonded joints by the asymptotic expansion method. Int. J. Eng. Sci. 1991, 29, 493–512. [Google Scholar] [CrossRef]
  21. Geymonat, G.; Krasucki, F.; Lenci, S. Mathematical Analysis of a bonded joint with a soft thin adhesive. Math. Mech. Solids 1999, 16, 201–225. [Google Scholar] [CrossRef]
  22. Bessoud, A.L.; Krasucki, F.; Serpilli, M. Asymptotic analysis of shell-like inclusions with high rigidity. J. Elast. 2011, 103, 153–172. [Google Scholar] [CrossRef]
  23. Ljulj, M.; Tambaca, J. 3D structure – 2D plate interaction model. Math. Mech. Solids 2019, 4, 3354–3377. [Google Scholar] [CrossRef]
  24. Lebon, F.; Rizzoni, R. Asymptotic analysis of a thin interface: The case involving similar rigidity. Int. J. Eng. Sci. 2010, 48, 473–486. [Google Scholar] [CrossRef] [Green Version]
  25. Lebon, F.; Rizzoni, R. Asymptotic behavior of a hard thin linear interphase: An energy approach. Int. J. Solids Struct. 2011, 48, 441–449. [Google Scholar] [CrossRef] [Green Version]
  26. Rizzoni, R.; Dumont, S.; Lebon, F.; Sacco, E. Higher order model for soft and hard elastic interfaces. Int. J. Solids Struct. 2014, 51, 4137–4148. [Google Scholar] [CrossRef]
  27. Dumont, S.; Rizzoni, R.; Lebon, F.; Sacco, E. Soft and hard interface models for bonded elements. Compos. Part B Eng. 2018, 153, 480–490. [Google Scholar] [CrossRef] [Green Version]
  28. Lebon, F.; Rizzoni, R. Higher order interfacial effects for elastic waves in one dimensional phononic crystals via the Lagrange-Hamilton’s principle. Eur. J. Mech. A Solids 2018, 67, 57–70. [Google Scholar] [CrossRef] [Green Version]
  29. Serpilli, M. Mathematical modeling of weak and strong piezoelectric interfaces. J. Elast. 2015, 121, 235–254. [Google Scholar] [CrossRef]
  30. Serpilli, M. Asymptotic interface models in magneto-electro-thermo-elastic composites. Meccanica 2017, 52, 1407–1424. [Google Scholar] [CrossRef]
  31. Serpilli, M. Classical and higher order interface conditions in poroelasticity. Ann. Solid Struct. Mech. 2019, 11, 1–10. [Google Scholar] [CrossRef]
  32. Serpilli, M. On modeling interfaces in linear micropolar composites. Math. Mech. Solids 2018, 23, 667–685. [Google Scholar] [CrossRef]
  33. Serpilli, M.; Rizzoni, R.; Lebon, F.; Dumont, S. An asymptotic derivation of a general imperfect interface law for linear multiphysics composites. Int. J. Solids Struct 2019, 180–181, 97–107. [Google Scholar] [CrossRef] [Green Version]
  34. Dumont, S.; Serpilli, M.; Rizzoni, R.; Lebon, F. Numerical Validation of Multiphysic Imperfect Interfaces Models. Front. Mater. 2020, 7, 158. [Google Scholar] [CrossRef]
  35. Shi, P.; Shillor, M. Existence of a solution to the N dimensional problem of thermoelastic contact. Indiana Univ. Math. J. 1992, 17, 1597–1618. [Google Scholar]
  36. Khludnev, M.; Kovtunenko, V.A. Analysis of Cracks in Solids; WIT-Press: Southampton, UK; Boston, MA, USA, 2000. [Google Scholar]
  37. Xu, X. The N-dimensional quasistatic problem of thermoelastic contact with Barbers heat exchange conditions. Adv. Math. Sci. Appl. 1996, 6, 559–587. [Google Scholar]
  38. Showalter, R.E. Diffusion in poro-elastic media. J. Math. Anal. Appl. 2000, 251, 310–340. [Google Scholar] [CrossRef] [Green Version]
  39. Miara, B.; Suarez, J.S. Asymptotic pyroelectricity and pyroelasticity in thermopiezoelectric plates. Asympt. Anal. 2013, 81, 211–250. [Google Scholar] [CrossRef]
  40. Bonaldi, F.; Geymonat, G.; Krasucki, F.; Serpilli, M. An asymptotic plate model for magneto-electro-thermo-elastic sensors and actuators. Math. Mech. Solids 2017, 22, 798–822. [Google Scholar] [CrossRef] [Green Version]
  41. Renard, Y.; Pommier, J. Getfem++. An Open Source Generic C++ Library for Finite Element Methods; Tech. Rep.; INSA Lyon: Villeurbanne, France, 2002. [Google Scholar]
  42. Geuzaine, C.; Remacle, J.-F. Gmsh: A 3-d finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Methods Eng. 2009, 79, 1309–1331. [Google Scholar] [CrossRef]
  43. Jafarinezhad, M.R.; Eslami, M.R. Coupled thermoelasticity of FGM annular plate under lateral thermal shock. Compos. Struct. 2017, 168, 758–771. [Google Scholar] [CrossRef]
Figure 1. Initial (a); rescaled (b); and limit (c) configurations of the composite.
Figure 1. Initial (a); rescaled (b); and limit (c) configurations of the composite.
Technologies 09 00017 g001
Figure 2. The 3D geometry of the thermoelastic laminated plate represented in the plane ( x 1 , x 3 ) .
Figure 2. The 3D geometry of the thermoelastic laminated plate represented in the plane ( x 1 , x 3 ) .
Technologies 09 00017 g002
Figure 3. Applied heat flow shock.
Figure 3. Applied heat flow shock.
Technologies 09 00017 g003
Figure 4. Convergence diagrams with respect to ε / h for: t = 1 s (left); and t = 10 s (right).
Figure 4. Convergence diagrams with respect to ε / h for: t = 1 s (left); and t = 10 s (right).
Technologies 09 00017 g004
Figure 5. Evolution of the relative error with respect to the time ( ε / h = 0.01 ).
Figure 5. Evolution of the relative error with respect to the time ( ε / h = 0.01 ).
Technologies 09 00017 g005
Figure 6. Displacement u 3 and temperature Θ along the X 3 -axis, on a fiber ( X ¯ 1 = 6 , X ¯ 2 = 6 , X 3 ) , for given times.
Figure 6. Displacement u 3 and temperature Θ along the X 3 -axis, on a fiber ( X ¯ 1 = 6 , X ¯ 2 = 6 , X 3 ) , for given times.
Technologies 09 00017 g006
Figure 7. Displacement field U = ( U i ) and temperature Θ versus time t, at a given point X ¯ = ( 6.5 , 3.5 , 0 ) .
Figure 7. Displacement field U = ( U i ) and temperature Θ versus time t, at a given point X ¯ = ( 6.5 , 3.5 , 0 ) .
Technologies 09 00017 g007
Figure 8. Stresses Σ 33 and Σ 13 along the X 3 -axis, on a fiber ( X ¯ 1 = 6 , X ¯ 2 = 6 , X 3 ) , for given times.
Figure 8. Stresses Σ 33 and Σ 13 along the X 3 -axis, on a fiber ( X ¯ 1 = 6 , X ¯ 2 = 6 , X 3 ) , for given times.
Technologies 09 00017 g008
Figure 9. Evolution in time of U 3 ( X ¯ , t ) at a given point X ¯ = ( 6.5 , 3.5 , 0 . ) for a homogeneous Al-plate and an Al/PVC composite plate.
Figure 9. Evolution in time of U 3 ( X ¯ , t ) at a given point X ¯ = ( 6.5 , 3.5 , 0 . ) for a homogeneous Al-plate and an Al/PVC composite plate.
Technologies 09 00017 g009
Table 1. Thermoelastic material properties for Al and PVC Foam.
Table 1. Thermoelastic material properties for Al and PVC Foam.
Material 1: AlMaterial 2: PVC Foam
ρ 1 2700[kg/m 3 ] ρ 2 250[kg/m 3 ]
e 1 72.4[GPa] e 2 0.28[GPa]
ν 1 0.32 ν 2 0.40
α 1 40.0[m/m K] α 2 22.4[m/m K]
k 1 122.2[W/mK] k 2 0.05[W/mK]
c v 1 900[J/kg K] c v 2 1900[J/kg K]
Table 2. Relative error u ε u L 2 u L 2 .
Table 2. Relative error u ε u L 2 u L 2 .
ε h / t 5101520
0.1 4.20 × 10 3 1.89 × 10 3 1.58 × 10 3 1.60 × 10 3
0.05 7.55 × 10 4 3.44 × 10 4 2.76 × 10 4 3.45 × 10 4
0.01 2.25 × 10 5 3.48 × 10 6 4.52 × 10 6 4.56 × 10 6
Table 3. Relative error θ ε θ L 2 θ L 2 .
Table 3. Relative error θ ε θ L 2 θ L 2 .
ε h / t 5101520
0.1 1.56 × 10 3 5.16 × 10 4 7.57 × 10 4 4.97 × 10 4
0.05 1.07 × 10 4 2.83 × 10 4 9.94 × 10 5 2.27 × 10 5
0.01 9.88 × 10 8 5.58 × 10 7 1.78 × 10 8 7.43 × 10 10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Serpilli, M.; Dumont, S.; Rizzoni, R.; Lebon, F. Interface Models in Coupled Thermoelasticity. Technologies 2021, 9, 17. https://0-doi-org.brum.beds.ac.uk/10.3390/technologies9010017

AMA Style

Serpilli M, Dumont S, Rizzoni R, Lebon F. Interface Models in Coupled Thermoelasticity. Technologies. 2021; 9(1):17. https://0-doi-org.brum.beds.ac.uk/10.3390/technologies9010017

Chicago/Turabian Style

Serpilli, Michele, Serge Dumont, Raffaella Rizzoni, and Frédéric Lebon. 2021. "Interface Models in Coupled Thermoelasticity" Technologies 9, no. 1: 17. https://0-doi-org.brum.beds.ac.uk/10.3390/technologies9010017

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop