Mahmoud Ali
email: mahmoud.ali@univ-orleans.fr

Thomas Sayet

Alain Gasser
email: alain.gasser@univ-orleans.fr

Eric Blond

Computational Homogenization of Elastic-Viscoplastic Refractory Masonry with Dry Joints

Refractory masonry with dry joints is widely used in the steel-making industry for the linings of several high-temperature components (> 1500 °C) including steel ladles and furnaces. To properly optimize the design and performance of these linings, thorough numerical models that consider the presence of joints, joints closure and reopening and the nonlinear elastic-viscoplastic behaviour (creep and stress relaxation) of refractories at high temperature are required. The present study reports on the formulation, numerical implementation, and application of a homogenized multi-scale elastic-viscoplastic model for the simulation of refractory masonry linings with dry joints. Refractory bricks are

considered to exhibit linear elasticity as well as rate-dependent plasticity. Four joint patterns are predefined based on the state of bed and head joints. The homogenized elastic-viscoplastic behaviour of each joint pattern is determined using finite element based nonlinear homogenization approach. The transition criterion between the four patterns are defined in terms of macroscopic stresses and strains. Verification of the developed homogenized constitutive laws is carried out by comparing the numerical results of the detailed micro models (brick and joints are considered) with the homogeneous equivalent material models. Furthermore, comparisons with experimental results of refractory masonry walls subjected to biaxial compression load at room and high temperature are carried out. Good agreements between the experimental and numerical results are obtained.

Then, the validated models are employed to predict the mechanical behavior of refractory masonry structures subjected to different loading conditions. The present numerical model is able to simulate the orthotropic, compressible, rate-dependent homogenized behaviour of mortarless refractory masonry structures, and accounts for joints closure and reopening due to loading and unloading.

Keywords

Introduction

Refractory masonry with dry joints is extensively used for the linings of several high temperature components such as steel ladles, furnaces and rotary kilns due to their high thermal, mechanical and chemical stability, good thermo-mechanical and thermo-chemical properties [START_REF] Blond | Thermo-chemomechanical modeling of refractory behavior in service: Key points and new developments[END_REF][START_REF] Coz Díaz | Design and finite element analysis of a wet cycle cement rotary kiln[END_REF][START_REF] Qi | Preparation and characterization of microporous mullite-corundum refractory aggregates with high strength and closed porosity[END_REF][START_REF] Smets | Behaviour of magnesia-carbon refractories in vacuum-oxygen decarburisation ladle linings[END_REF][START_REF] Wu | A comparative study on the microstructures and mechanical properties of a dense and a lightweight magnesia refractories[END_REF]. When used in steel ladles, they are subjected to high thermal gradients, thermal shock, cyclic thermal heating and cooling, high thermo-mechanical stresses, slag attack, and harsh chemical environment [START_REF] Blond | Effect of slag impregnation on thermal degradations in refractories[END_REF][START_REF] Fruhstorfer | Erosion and corrosion of alumina refractory by ingot casting steels[END_REF][START_REF] Andreev | Role of fatigue in damage development of refractories under thermal shock loads of different intensity[END_REF][START_REF] Chen | Corrosion and penetration behaviors of slag/steel on the corroded interfaces of Al2O3-C refractories: Role of Ti3AlC2[END_REF][START_REF] Dai | Corrosion mechanism and protection of BOF refractory for high silicon hot metal steelmaking process[END_REF][START_REF] Zou | Slag corrosion-resistance mechanism of lightweight magnesia-based refractories under a static magnetic field[END_REF]. Therefore, the lifespan of steel ladles internal linings is very short (several castings). Moreover, the design and optimization of these linings is still an engineering challenge due to the complex interaction between thermal fields, chemistry and the nonlinear mechanical behaviour of refractories at high temperature [START_REF] Blond | Thermo-chemomechanical modeling of refractory behavior in service: Key points and new developments[END_REF][START_REF] Blond | Effect of slag impregnation on thermal degradations in refractories[END_REF][START_REF] Samadi | Statistical study of compressive creep parameters of an alumina spinel refractory[END_REF].

In a typical industrial scale steel ladle, the internal linings are built with dry joints masonry (also called mortarless masonry). Thousands of refractory bricks are arranged periodically. Between the contact surfaces of the bricks, there are small gaps called dry joints. The macroscopic mechanical behaviour of these structures is orthotropic, nonlinear and different from the behaviour of the material of the bricks due to the following reasons:

first, material discontinuity introduced by the presence of joints, and second, joints closure and reopening caused by cyclic loading and unloading [START_REF] Gasser | Modelling of joint effects on refractory lining behaviour[END_REF][START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF]. To properly optimize the design and performance of these linings, rigorous numerical models are required especially when experimental work remains a challenging topic due to technical complexities and severe working conditions. The numerical model should consider the orthotropic and nonlinear behaviour of the linings at high temperature and should be efficient in terms of computational time and cost to enable full scale structural modelling of the linings.

However, most of previous studies on mortarless masonry typically focused on masonry used in civil structures (usually operates at room temperature) [START_REF] Zahra | Characterisation and strategies for mitigation of the contact surface unevenness in dry-stack masonry[END_REF][START_REF] Bui | Discrete element modelling of the inplane and out-of-plane behaviour of dry-joint masonry wall constructions[END_REF][START_REF] Chewe Ngapeya | Overcome of bed-joint imperfections and improvement of actual contact in dry-stacked masonry[END_REF][START_REF] Nodargi | A variational-based fixed-point algorithm for the limit analysis of dry-masonry block structures with non-associative Coulomb friction[END_REF][START_REF] Lourenço | Dry joint stone masonry walls subjected to in-plane combined loading[END_REF][START_REF] Thamboo | Development of design methodology for mortarless masonry system: Case study -a resettlement housing colony[END_REF][START_REF] Martínez | Experimental and numerical evaluation of reinforced dry-stacked concrete masonry walls[END_REF][START_REF] Smoljanović | A finite-discrete element model for dry stone masonry structures strengthened with steel clamps and bolts[END_REF][START_REF] Tempesta | Safety evaluation of masonry arches. A numerical procedure based on the thrust line closest to the geometrical axis[END_REF]. Therefore, the bricks are considered to exhibit isotropic linear elasticity. For example, several experimental and numerical studies are carried out to investigate the impact of contact surfaces roughness on the contact area between the bricks [START_REF] Zahra | Characterisation and strategies for mitigation of the contact surface unevenness in dry-stack masonry[END_REF][START_REF] Chewe Ngapeya | Overcome of bed-joint imperfections and improvement of actual contact in dry-stacked masonry[END_REF] or load bearing capacity of the wall [START_REF] Chewe Ngapeya | Experimental and analytical analysis of the load-bearing capacity Pu of improved dry-stacked masonry[END_REF]. In these studies, the micro finite element modelling (FEM) approach is adopted (i.e. bricks and joints are considered and only few bricks and joints are used as a physical model of the problem). Also, limit analysis and discrete element modelling of the collapse mechanism of masonry walls subjected to in-plane and out-of-plane loading conditions can be found in the literature [START_REF] Bui | Discrete element modelling of the inplane and out-of-plane behaviour of dry-joint masonry wall constructions[END_REF][START_REF] Nodargi | A variational-based fixed-point algorithm for the limit analysis of dry-masonry block structures with non-associative Coulomb friction[END_REF][START_REF] Foti | DEM modeling and experimental analysis of the static behavior of a dry-joints masonry cross vaults[END_REF][START_REF] Smoljanović | Numerical analysis of 3D drystone masonry structures by combined finite-discrete element method[END_REF][START_REF] Casapulla | In-plane frictional resistances in dry block masonry walls and rocking-sliding failure modes revisited and experimentally validated[END_REF][START_REF] Portioli | Large displacement analysis of dry-jointed masonry structures subjected to settlements using rigid block modelling[END_REF]. It should be noted that the micro FEM approach is only suitable for modelling few bricks and joints because considering all the bricks and joints of a wall (hundreds of bricks) leads rapidly to too expensive computational costs and solution convergence problems.

In the case of refractory masonry, very limited number of studies could be found in the literature. Room and high temperature compression tests on two refractory bricks were carried out to investigate the impact of applied stress level on the gradual closure of joints and the contact area between the bricks [START_REF] Allaoui | Digital image correlation measurements of mortarless joint closure in refractory masonries[END_REF][START_REF] Andreev | Compressive behaviour of dry joints in refractory ceramic masonry[END_REF]. It has been found that the contact area between the bricks increases with the increase of the applied compression stress.

Furthermore, dry joints close gradually (i.e. a decrease in joint thickness) with the gradual increase of the applied stress. With regard to numerical studies, Hou et al. [START_REF] Hou | A Method for Steel Ladle Lining Optimization Applying Thermomechanical Modeling and Taguchi Approaches[END_REF] developed a 2D micro FE model to investigate the effect of linings thickness and material properties on the thermo-mechanical behaviour of steel ladle refractory linings. Due to the selected micro FEM approach, only 2D cut of few bricks and one joint is considered. Nguyen et al. [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF] developed and validated a homogeneous equivalent material model, by replacing the bricks and joints with an equivalent material, for the simulation of mortarless refractory masonry structures. The bricks are considered to obey isotropic linear elasticity. The model considers the influence of joints closure on the homogeneous mechanical response of the structure. Their results indicate that the macroscopic behaviour of the masonry structure is orthotropic and nonlinear due to gradual closure of joints. It should be noted that the macroscopic model developed by Nguyen et al. [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF] can be used to simulate large scale structures, such as steel ladles and furnaces, with reasonable computational costs.

From the above discussion, it could be inferred that most of previous studies on masonry with dry joints focus on civil structure masonry using micro FEM approach and consider elastic behaviour of the bricks. However, these models are not suitable to model large sized structures such as steel ladles and furnaces. In addition, only one multi scale numerical model has been developed in the literature and it is suitable to simulate only the elastic or thermoelastic behaviour of the masonry. At high temperature, elastic models usually leads to overestimation of stress levels due to the elastic viscoplastic behaviour of refractories at high temperature [START_REF] Blond | Thermo-chemomechanical modeling of refractory behavior in service: Key points and new developments[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF]. More studies to understand the effects of joints, joints closure and reopening combined with viscoplasticity on the mechanical behaviour of refractory masonry linings at high temperature are required.

The present paper aims at performing computational nonlinear homogenization of elastic viscoplastic mortarless refractory masonry structures and developing a multi scale numerical model for their simulation. The developed model is efficient in terms of computational cost as it replaces all the bricks and joints with an elastic viscoplastic homogeneous equivalent material model. The present work is an extension of the previous work developed by Nguyen et al. and Ali et al. [14,[START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF]. As compared to the work in the literature [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF], the addition of viscoplasticty to the developed multi scale numerical model promotes better understanding of the nonlinear mechanical behaviour of refractory masonry structures at high temperature and, therefore, assists in the optimization of their design and performance.

The present paper is organized as following: in section 2, description of the mortarless refractory masonry structure is presented. In section 3, periodic nonlinear computational homogenization of the refractory mortarless masonry structure and identification of the effective elastic and viscoplastic parameters are presented and discussed. Verification of the developed models and results of masonry structures subjected to different loading cases are presented and discussed in section 5. Key findings of the present study and conclusions are given in section 6.

Mortarless masonry structure

In the present study, the mortarless refractory masonry structure illustrated in figure 1 is considered. The height and diameter of a typical industrial steel ladle are around 5 and 4.5 m, respectively [START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF]. The internal linings (also called working lining) of the ladle are built up from thousands of tapered shape refractory bricks. However, for the sake of Two types of joints are defined on the basis of their orientation: bed joints with initial thickness g 0,bed (horizontal joints) and head joints with initial thickness g 0,head (vertical joints). Under cyclic thermal heating and cooling or mechanical loading and unloading, these joints can close and reopen. It should be noted that the mechanical response of mortarless masonry structures varies with the closure and reopening of joints. Previous experimental studies show that the stiffness of mortarless masonry structure increases with the closure of joints [START_REF] Allaoui | Digital image correlation measurements of mortarless joint closure in refractory masonries[END_REF][START_REF] Andreev | Compressive behaviour of dry joints in refractory ceramic masonry[END_REF]. This phenomenon should be considered when developing thorough numerical models for the simulation of mortarless masonry linings. For the purpose of considering the impact of joints closure and reopening on the homogenized elastic-viscoplastic behaviour of mortarless masonry structures, four possible joint patterns are predefined (see figure 2). Each pattern is defined based on the state of both bed and head joints (i.e. open or closed) and represents different periodic masonry structure with different equivalent behaviour. The four joint patterns are as follows [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF]:

1. Pattern C: head and bed joints are closed (all joints are closed). In the present study, bauxite-based refractory ceramic is considered as the base material of the bricks. At high temperature, the constitutive material is assumed to undergo small deformations (i.e. the displacement gradients are small, each component is small as compared to unity, |u i,j | << 1) and to exhibit an isotropic linear elasticity as well as an isotropic steady state rate-dependent plasticity (secondary creep). Under these assumptions, the total strain second order tensor (ε t ) can be described by the symmetric part of the displacement vector (u) with respect to the position vector (X) of the material point according to:

ε t = 1 2   ∂u ∂X + ∂u ∂X T   (1)
The total strain tensor can be additively decomposed into elastic (ε e ) and viscoplastic (ε vp ) second order strain tensors according to:

ε t = ε e + ε vp (2)
The second order linear elastic strain tensor can be determined using Hooke's law for isotropic linear elasticity according to:

ε e = 1 + ν Y σ - ν Y tr(σ)I (3)
where, ν is the Poisson's ratio of the bulk material of the bricks and Y is the Young's modulus of the bulk material of the bricks. σ and I refer respectively to the second order local stress tensor and the second order identity tensor.

The steady state isotropic viscoplastic strain rate second order tensor ( εvp ) is expressed as [START_REF] Lemaitre | Mechanics of Solid Materials[END_REF]:

εvp = 3 2 A (σ eq ) n-1 σ d (4)
where A, n are the creep parameters of the bulk material of the bricks. σ eq and σ d denote the Von Mises equivalent stress and the second order deviatoric stress tensor, respectively.

The deviatoric stress tensor and Von Mises equivalent stress are expressed as:

σ d = σ - 1 3 tr(σ)I ( 5
)
σ eq = 3 2 σ d : σ d (6)
The required material parameters are taken from [START_REF] Blond | Modelling of high temperature asymmetric creep behavior of ceramics[END_REF] and reported in table 1. For the sake of simplicity, no strain hardening of the constitutive material of the bricks is considered. In addition, symmetric creep behaviour is assumed (i.e. creep parameters in tension and compression are equal). Further details on the nonlinear homogenization analysis and the identification of the homogenized elastic viscoplastic behaviour of the four patterns and the transition criteria from one joint pattern to another are described in sections 3 and 4. 3 Periodic nonlinear homogenization of mortarless masonry structures

Joint pattern C

As indicated previously in the case of joint pattern C, head and bed joints are closed (see figure 2). Previous experimental studies show that when joints are closed, the mechanical behaviour of the masonry is similar to the mechanical behaviour of the base material of the bricks. Thus, the macroscopic elastic-viscoplastic behaviour of pattern C is isotropic, similar to that of the base material of the bricks and can be described by equations 1 -6.

Joint pattern B

Micro modelling

In the case of pattern B, head joints are open, and bed joints are closed (see figure 2). The existence of open joints in the masonry results in a reduction in their effective stiffness [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF][START_REF] Ali | Thermomechanical modelling of refractory mortarless masonry wall subjected to biaxial compression[END_REF][START_REF] Rekik | Numerical homogenization model for effective creep properties of microcracked masonry[END_REF]. Thus, the macroscopic elastic-viscoplastic behaviour of pattern B is different from that of the base material of the bricks. For the purpose of studying the homogenized mechanical behaviour of periodic structures, like pattern B, and computing the effective elastic viscoplastic parameters, the finite element (FE) -based homogenization technique can be adopted [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Tsuda | Duplex model for homogenized elastic-viscoplastic behavior of plate-fin structures at high temperatures[END_REF][START_REF] Goto | Elasto-viscoplastic analysis for negative through-the-thickness Poisson's ratio of woven laminate composites based on homogenization theory[END_REF][START_REF] Tikarrouchine | Fully coupled thermo-viscoplastic analysis of composite structures by means of multi-scale three-dimensional finite element computations[END_REF].

To carry out homogenization of joint pattern B using FE technique, a periodic unit cell (UC) with volume (V UC ) has been selected from the periodic structure as illustrated in figure 3. Then, 3D-FE simulations have been performed on the UC to characterise its homogenized elastic viscoplastic response. To account for the periodicity of the structure, periodic boundary conditions have been applied to the UC's boundaries [START_REF] Xia | A unified periodical boundary conditions for representative volume elements of composites and applications[END_REF][START_REF] Omairey | Development of an ABAQUS plugin tool for periodic RVE homogenisation[END_REF][START_REF] Henyš | Comparison of current methods for implementing periodic boundary conditions in multi-scale homogenisation[END_REF]. In contrast to homogeneous boundary conditions, periodic boundary conditions satisfy both boundary displacement periodicity and boundary traction periodicity. Moreover, homogeneous boundary conditions are over constrained conditions, violate stress periodicity of the boundaries and, therefore, are not recommended to be used for periodic UC subjected to shear loading conditions [START_REF] Xia | A unified periodical boundary conditions for representative volume elements of composites and applications[END_REF][START_REF] Omairey | Development of an ABAQUS plugin tool for periodic RVE homogenisation[END_REF]. Several FE numerical experiments of uniaxial tensile tests along the x, y and z directions and simple shear tests in the xy, xz, and yz planes (see figure 4) have been performed.

The numerical models have been developed using ABAQUS software and the UC has been meshed with 3D quadrilateral elements with 1 mm size. For all cases, the UC has been subjected to an average or macroscopic stress (shear or normal) of 1.5 MPa.

The base material of the bricks is considered to exhibit an isotropic linear elasticity, an isotropic steady state rate-dependent plasticity (secondary creep) and to obey the elastic and viscoplastic constitutive equations given in section 2.

To characterise the homogenized elastic behaviour of the UC and to avoid the contribution of the viscoplastic strain when calculating the effective elastic stiffness of the structure, first, purely elastic numerical simulations have been carried out. Then, viscoplastic numerical analysis have been performed on the UC to characterise its homogenized viscoplastic response. From figure 4, it can be seen that the periodicity of boundary displacements and stresses are satisfied due to the applied periodic boundary conditions. From the uniaxial tensile and simple shear numerical tests of the UC subjected to macroscopic normal or shear stress of 1.5 MPa, the macroscopic strain tensor and the effective material parameters can be found.

Macroscopic constitutive equations

Macroscopic elastic behaviour

The first goal is to define a relation between macroscopic stress tensor and macroscopic elastic strain tensor (i.e. determine the effective elastic stiffness tensor of pattern B).

Due to the existence of open joints, the macroscopic elastic behaviour of pattern B is orthotropic and can be described using equation 7. In the case of orthotropic linear elastic materials, the macroscopic stress tensor is related to the macroscopic elastic strain tensor through the macroscopic form of Hooke's law as following [START_REF] Reddy | Mechanics of laminated composite plates and shells: theory and analysis[END_REF]:

Σ = C e : E e (7)
where Σ and E e are the second order macroscopic stress and elastic strain tensors while, C e is the fourth order macroscopic elastic stiffness tensor.

In equation 7, the second and fourth order tensors are symmetric and, therefore, they can be respectively reduced to 6 × 1 and 6 × 6 matrices using the Voigt notations [START_REF] Tikarrouchine | Fully coupled thermo-viscoplastic analysis of composite structures by means of multi-scale three-dimensional finite element computations[END_REF].

The effective elastic fourth order tensor is characterized by 9 nonzero components and can be written in a matrix form, with respect to the principal material coordinate system, using Voigt notations as follows [START_REF] Reddy | Mechanics of laminated composite plates and shells: theory and analysis[END_REF]:

C e =                    
0 0 0 0 0 0 G 12 0 0 0 0 0 0 G 13 0 0 0 0 0 0 G 23                     [M P a], (8) with,
= 1 -ν 12 ν 21 -ν 23 ν 32 -ν 31 ν 13 -2 ν 12 ν 23 ν 31 E 11 E 22 E 33 (9)
In the above equations, E and ν denote respectively the macroscopic or effective Young's modulus and the macroscopic or effective Poisson's ratio. G is the macroscopic or effective shear modulus.

According to the average theorems, when a uniform stress is applied to the boundary surfaces of the UC, the macroscopic stress is equal to the applied stress (i.e. σ = Σ) [START_REF] Tikarrouchine | Fully coupled thermo-viscoplastic analysis of composite structures by means of multi-scale three-dimensional finite element computations[END_REF].

Another possible approach for calculating the macroscopic stress is to integrate the local stress fields over the volume of the UC (V UC ) according to Hill's definition as follows [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Iltchev | Computational homogenisation of periodic cellular materials: Application to structural modelling[END_REF][START_REF] Hill | The essential structure of constitutive laws for metal composites and polycrystals[END_REF]:

Σ = σ = 1 V UC V UC σdV ( 10
)
From the simulated combination of uniaxial tensile and simple shear tests, the local strains can be obtained from solving the boundary value problem of the UC subjected to average normal or shear stress of 1.5 MPa. From the local elastic strains, the macroscopic elastic strains can be calculated. Then, using equations 7 -9, the 9 nonzero components of C e can be found and the relation between the macroscopic stress tensor and macroscopic elastic strain tensor can be established.

In previous studies, two main approaches have been used to calculate the macroscopic strain fields [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Iltchev | Computational homogenisation of periodic cellular materials: Application to structural modelling[END_REF][START_REF] Hori | On two micromechanics theories for determining micro-macro relations in heterogeneous solids[END_REF][START_REF] Fritzen | Computational homogenization of elasto-plastic porous metals[END_REF][START_REF] Matsuda | Effects of fiber distribution on elastic-viscoplastic behavior of long fiber-reinforced laminates[END_REF][START_REF] Molla | Efficient modeling of metallic interconnects for thermo-mechanical simulation of SOFC stacks: Homogenized behaviors and effect of contact[END_REF][START_REF] Ohno | Homogenized elastic-viscoplastic behavior of plate-fin structures with two pore pressures[END_REF]. The first approach is based on calculating the macroscopic total strain by integrating the local total strain (ε t ) over the volume of the UC according to [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Hori | On two micromechanics theories for determining micro-macro relations in heterogeneous solids[END_REF][START_REF] Fritzen | Computational homogenization of elasto-plastic porous metals[END_REF][START_REF] Matsuda | Effects of fiber distribution on elastic-viscoplastic behavior of long fiber-reinforced laminates[END_REF]:

E t = ε t = 1 V UC V UC εdV (11)
In the second approach, the macroscopic elastic strains are calculated from the average change in the displacements of the UC's corners (i.e. by dividing the change in displacements of the corners by the initial dimensions of the UC) [START_REF] Iltchev | Computational homogenisation of periodic cellular materials: Application to structural modelling[END_REF][START_REF] Molla | Efficient modeling of metallic interconnects for thermo-mechanical simulation of SOFC stacks: Homogenized behaviors and effect of contact[END_REF][START_REF] Ohno | Homogenized elastic-viscoplastic behavior of plate-fin structures with two pore pressures[END_REF]. If there are neither gaps nor cracks in the UC, both techniques converge to the same value. However, in the case of joint pattern B, the macroscopic elastic strains cannot be computed by integrating the micro (or local) strains over the volume of the UC (first approach). As, the presence of open joints within the UC leads to discontinuous displacement fields (see figure 5), because there is neither viscosity nor rigidity in the volume of the gaps (i.e. local strains within the gaps are unknown) [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Iltchev | Computational homogenisation of periodic cellular materials: Application to structural modelling[END_REF]. Therefore, in the present study, the second approach has been used to calculate the macroscopic strains. The calculated 9 nonzero components of C e are reported in table 2. Now, the relationship between macroscopic stresses and macroscopic elastic strains is established. Further details on defining the macroscopic viscoplastic behaviour law are given below.

Macroscopic viscoplastic behavior

The second goal is to define a macroscopic constitutive law to describe the viscoplastic The orthotropic homogenized steady-state viscoplastic behaviour of heterogeneous solids can be described using the macroscopic form of equation 4 developed by Tsuda et al. [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Tsuda | Duplex model for homogenized elastic-viscoplastic behavior of plate-fin structures at high temperatures[END_REF][START_REF] Ohno | Homogenized elastic-viscoplastic behavior of plate-fin structures with two pore pressures[END_REF][START_REF] Ohno | Homogenized elastic-viscoplastic behavior of anisotropic open-porous bodies with pore pressure[END_REF]. The main advantage of this constitutive law is that using the same parameters of the base material, the orthotropic homogenized viscoplastic behaviour of heterogeneous solids can be found. The constitutive law is obtained through averaging the local fields (σ, ε vp , σ eq ) over the volume of the UC and using a localization tensor according to:

Ė vp = 1 2 A (Σ eq ) n-1 N : Σ (12)
with A and n are materials parameters of the bricks, reported in table 1, Σ eq and N are the macroscopic equivalent stress and a fourth order tensor with the same meaning of the localization tensor. The N tensor accounts for the orthotropy and compressibility (in the x direction, see figure 6) of the structure and enables using the creep parameters of the constitutive material (i.e. bridge between the micro and macro scales). The N tensor is symmetric and can be represented by 6 × 6 matrix using the Voigt notations.

N is characterized by 9 nonzero components and can be represented, with respect to the principal material coordinate system, as:

N =                     N 11 N 12 N 13 0 0 0 N 12 N 22 N 23 0 0 0 N 13 N 23 N 33 0 0 0 0 0 0 N 44 0 0 0 0 0 0 N 55 0 0 0 0 0 0 N 66                     (13)
For heterogeneous solids, the macroscopic equivalent stress is written in terms of the macroscopic stress tensor and N as [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF][START_REF] Ohno | Homogenized elastic-viscoplastic behavior of plate-fin structures with two pore pressures[END_REF][START_REF] Xue | Constitutive model for quasi-static deformation of metallic sandwich cores[END_REF]:

Σ eq = 1 2 Σ T : N : Σ (14)
Here Σ T is the transpose of the macroscopic stress second order tensor. The macroscopic equivalent viscoplastic strain rate ( Ėvp eq ) is defined in terms of N and Ė vp as follows [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF]:

Ėvp eq = 2 Ė vp T : N -1 : Ė vp ( 15
)
By combining equations 12 -14 and using the combinations of uniaxial and simple shear finite element numerical tests presented before, the 9 non-zero components of the tensor N can be calculated and the macroscopic creep law can be defined.

Under uniaxial loading and periodic boundary conditions the shear components of the macroscopic stress tensor vanish. Therefore, the first 6 components of N can be defined as follows [START_REF] Tsuda | Homogenized elasticviscoplastic behavior of plate-fin structures at high temperatures: Numerical analysis and macroscopic constitutive modeling[END_REF]:

-Uniaxial tension along the x-direction

N 11 = 2 1 Σ xx 2 Σ xx Ėvp xx A 2 n+1
(16)

N 12 = Ėvp yy Ėvp xx N 11 ( 17
)
N 13 = Ėvp zz Ėvp xx N 11 (18)
-Uniaxial tension along the y-direction:

N 22 = 2 1 Σ yy 2 Σ yy Ėvp yy A 2 n+1 (19) N 12 = Ėvp xx Ėvp yy N 22 ( 20
)
N 23 = Ėvp zz Ėvp yy N 22 (21)
-Uniaxial tension along the z-direction:

N 33 = 2 1 Σ zz 2 Σ zz Ėvp zz A 2 n+1 (22) N 13 = Ėvp xx Ėvp zz N 33 (23)
N 23 = Ėvp yy Ėvp zz N 33 (24)
Under shear loading and periodic boundary conditions, the normal components of Σ vanish and the shear components of N are written as follows:

-Simple shear in the xy-plane:

N 44 = 2 1 Σ xy 2 Σ xy Ėvp xy A 2 n+1 (25)
-Simple shear in the xz-plane:

N 55 = 1 Σ xz 2 Σ xz Ėvp xz A 2 n+1 (26)
-Simple shear in the yz-plane:

N 66 = 2 1 Σ yz 2 Σ yz Ėvp yz A 2 n+1 (27)
In the previous equations, the values of macroscopic normal and shear stresses are equal to the applied average stresses. The macroscopic viscoplastic strain rates are determined from the time derivative of the macroscopic viscoplastic strain (E vp calculated from the displacement of the corners of the UC as described before). The obtained values of the 9 none-zero components of N are reported in table 2. It can be observed that the y and z components of the tensor are almost the same (due to the absence of joints in the z direction, only one brick in the depth of the wall, and the absence of open joints in the y direction). To compare the components of the N tensor of joint pattern B and those of the constitutive material, the components of N tensor have been identified for joint pattern C, using the same approach, and reported in table 2. It can be seen that, in the case of pattern B, the y and z components of N tensor are almost equal to those of the constitutive material of the bricks.

Joint pattern O

As presented earlier, in the case of pattern O, both head and bed joints are open. It can be seen from figure 2 that the masonry structure is composed of an array of separated bricks and the masonry structure is disconnected. For the purpose of computing the effective mechanical parameters of the structure (elastic and viscoplastic), one cannot use finite element-based homogenization approach (micro modelling) presented earlier because the boundary value problem of the UC is not clearly defined. However, the mechanical parameters of the constitutive material can be used to define the effective mechanical properties of joint pattern O [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF][START_REF] Ali | Thermomechanical modelling of refractory mortarless masonry wall subjected to biaxial compression[END_REF]. Due to the existence of open joints in the x and y directions, the effective elastic and viscoplastic mechanical parameters of the masonry structure in the x and y directions are very small. Therefore, the effective Young's modulus and the components of N in directions x and y are zero ( E x = E y = 0, see the coordinate system in figure 3). In addition, the effective Young's modulus and N components in z direction are equal to those of the constitutive material due to the absence of joints in the z direction (only one brick in the z direction). The C e and N tensors of joint pattern O are reported in table 2. To facilitate the numerical computations and to avoid numerical singularities and solution divergence, a very small value has been assigned to E x , E y , ν xy , ν xz , ν yz and G xy instead of zero.

Joint pattern H

In the case of joint pattern H, head joints are closed whereas bed joints are open (see figure 2). The masonry structure is composed of an array of separated courses of bricks (in the y-direction). Thus, the structure has zero macroscopic stiffness in the y-direction (due to the presence of open joints) while it has stiffness in the x (as all joints are closed in this direction) and z-directions (only one brick in the z direction) [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF][START_REF] Ali | Thermomechanical modelling of refractory mortarless masonry wall subjected to biaxial compression[END_REF]. Similar to pattern O, one cannot use finite element-based homogenization approach presented earlier because the boundary value problem of the UC is not clearly defined. However, the mechanical parameters of the constitutive material can be used to define the effective mechanical properties of joint pattern H. Since bed joints are open, the macroscopic stiffness of the structure in direction y is very small and E y is zero. Thus, the components of N in the y direction are also equal to zero. However, the effective Young's modulus and N components in x and z directions are equal to those of the constitutive material due to the absence of joints in the z direction and the absence of open joints in x direction. The C e and N tensors of joint pattern H are reported in table 2.

Joints closure and reopening criteria

As discussed above, each joint pattern has different homogeneous elastic viscoplastic response. In normal operating conditions, the masonry is subjected to cyclic thermal or mechanical loading and unloading. Thus, head or bed or both joints may close and reopen.

Consequently, the masonry structure changes from one joint pattern to another leading to a change in the macroscopic elastic viscoplastic behaviour of the structure. In the present work, this change has been taken into account by using suitable joints closure, reopening and pattern transition criteria.

Before loading, head and bed joints are open and the masonry structure can be fully described by pattern O. Under compression loads, the thickness of head or bed or both joints decreases gradually from the initial value (g 0 , around 0. As the equivalent material properties of the masonry are piecewise constant, the displacement increment at every point in the structure linearly depends on the increment in the macroscopic quantities [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Ali | Transient thermo-mechanical analysis of steel ladle refractory linings using mechanical homogenization approach[END_REF]. Thus, the instantaneous thickness of the head and bed joints can be written in terms of the macroscopic components of total strains in the x (E t xx ) and y (E t yy ) directions as:

g head = g 0, head + M x E t xx ( 29
)
g bed = g 0, bed + M y E t yy (30)
Where g bed and g head are the instantaneous thickness of bed and head joints, respectively. g 0,bed and g 0,head denote the initial thickness of bed and head joints, respectively. M y and M x are parameters with the same meaning of the localization tensor, and they depend on the dimensions of the brick (M y = h b and M x = l b ). E t yy and E t xx are the macroscopic total strains (elastic and viscoplastic strains) in the y and x directions.

With regard to joints reopening criterion, head and bed joints can reopen if the normal stress to the surface of the joint (head or bed) is higher than zero (i.e. tensile stress). In

Results and discussion

Comparisons between the micro and macro modelling

The main aim of this section is to assess the accuracy of the developed macroscopic elastic viscoplastic numerical model. For this purpose, the equivalent material model (hereafter will be called macro model) has been implemented in the finite element code, Abaqus [START_REF]ABAQUS/standard: User's Manual[END_REF],

with a user material subroutine UMAT. Then, it has been used to simulate uniaxial, simple shear and biaxial numerical experiments. Two different loading conditions have been considered namely, constant stress (creep) and constant strain loading (stress relaxation).

For both loading conditions, comparisons between results from the macro and micro models are performed. In the case of macro modelling, the bricks and the joints are replaced by an equivalent homogeneous medium and the macro constitutive laws (equations 7 -27) are used. Whereas, in the case of micro modelling, the bricks and joints are considered and the micro constitutive laws (equations 1 -6) are used.

For joint pattern B, comparisons between average total strains (elastic and viscoplastic strains) obtained from the micro and macro models for both constant stress uniaxial tensile load in the x, y and z directions and simple shear load with respect to the xy, xz and yz planes are shown in figure 8. The comparisons are presented for three constant stress levels in the respective direction. Good agreements between the results obtained from the macro and micro models can be observed. In addition, it can be observed that the macroscopic elastic viscoplastic behaviour of pattern B is orthotropic due to the presence of open joints.

In It should be noted that the previous constant stress loading conditions (uniaxial and simple shear) are used to identify the elastic viscoplastic parameters of the structure and, therefore, the agreements between the micro and macro models results are predicted. To In the second loading case, comparisons between average stresses obtained from the micro and macro models for constant strain rate uniaxial tensile loading in the x, y and z directions and simple shear with respect to the xy, xz and yz planes are shown in figure 10. The comparisons are presented for three constant strain rate levels. Good agreements between the results obtained from the macro and micro models can be observed. Overall, the stresses increase sharply due to load application and then, during holding, there is a decay due to stress relaxation. Again, similar values of macroscopic stresses have been observed in the case of uniaxial loading in the y and z directions. However, Σ yy and Σ zz are higher than Σ xx , because the effective stiffness of the structure in the x direction is lower than that in the y and z directions.

Modelling of masonry structures subjected to different loading conditions

The previous results have shown the capability of the developed model to describe the orthotropic and rate-dependent homogenized behaviour of the four periodic joint patterns of mortarless masonry structure described before. The focus now is to model masonry structures subjected to different loading conditions and to take into account joints closure and reopening due to loading and unloading and their impact on the elastic viscoplastic behaviour of masonry structures.

In order to validate the present numerical models, comparative studies with experimental results reported in [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Triessnig | Determination of Thermo-Mechanical Properties at the RHI Refractories Technology Center Leoben[END_REF][START_REF] Prietl | The evaluation of refractory linings thermo-mechanical properties[END_REF][START_REF] Prietl | Determination of material-specific parameters of refractory materials and linings under uniaxial and biaxial load conditions for the non-ferrous metal industry[END_REF] For the purpose of measuring the displacement in the x (perpendicular to the contact surfaces of head joints) and y (perpendicular to the contact surfaces of bed joints) directions, four linear variable differential transformers (LVDT) are used (see Fig. 11).The testing procedures were as following: first a pre-compression load was applied to the two moving ceramic plates at the same time and then, it was stopped in the direction for which the LVDTs detected a displacement while continuing in the second direction until the LVDTs, in the corresponding direction, detected a displacement. Finally, monotonically increasing loads were applied to both directions at the same time. For the high temperature test, the test field was covered by a heating hood. The temperatures of the top face (hot face) and bottom face (cold face) were measured. During heating, the masonry wall was free to expand. Once the wall reached the thermal steady state, the compression loads were applied similar to the room temperature test. The temperatures of the hot and cold faces during load application were around 1200 °C and 800°C, respectively. Further details about the experimental setup are reported in references [START_REF] Nguyen | Mechanical homogenisation of masonry wall without mortar[END_REF][START_REF] Triessnig | Determination of Thermo-Mechanical Properties at the RHI Refractories Technology Center Leoben[END_REF][START_REF] Prietl | The evaluation of refractory linings thermo-mechanical properties[END_REF][START_REF] Prietl | Determination of material-specific parameters of refractory materials and linings under uniaxial and biaxial load conditions for the non-ferrous metal industry[END_REF]. The finite element models of the room and high temperature tests were developed using ABAQUS software and the wall was meshed with 3D quadrilateral elements. The four ceramic loading plates (fixed and moving) are replaced by rigid plates. The masonry (bricks and dry joints) is replaced by an equivalent material model whose mechanical properties depend on the state of bed and head joints (open or closed). The effective elastic and viscoplastic parameters of the four joint patterns are computed using FE based nonlinear homogenization technique, described in section 3.2, and are reported in tables A1 and A2 in the appendix. For room temperature test, the mechanical behavior of the masonry is assumed to be elastic while, for high temperature test it was assumed to be elastic-viscoplastic. In both cases, the macroscopic behavior changes with the change of joint pattern (i.e., closure of bed, head or both joints).

Comparisons between experimental and numerical stress-strain diagrams of masonry walls subjected to biaxial compression at room temperature and 1200 °C are reported in figure 12. The strains, in the x and y directions, are calculated from the average displacements measured by the LVDTs and the stresses were determined from the resulting reaction forces of the two moving ceramic loading beams. It can be noticed that the present multi scale numerical model can reproduce, with reasonable accuracy, the mechanical behavior of the masonry structure under biaxial compression at either room temperature or 1200 °C. It can be noticed from figure 12 that the macroscopic mechanical behavior in both

x and y directions is orthotropic. This can be attributed to that the number of joints in the x and y directions, the length and height of the bricks are different. Furthermore, the macroscopic behavior of the wall is nonlinear due to change of the macroscopic mechanical parameters of the masonry with the gradual closure of the joints. As compared to room temperature test, the strains (and displacements, see figures 13 and 14) in the x and y directions are higher. This can be attributed to that the decrease of the material stiffness During the holding step, bed joints remain closed (figure 17-c). Finally, during unloading, few bed joints near the fixed rigid plate reopen and, thus, the masonry structure change from pattern B back to pattern O. Similarly, in the case of x-LHU, with the gradual increase of the applied displacement in the x direction, head joints close gradually while bed joints remain open and, therefore, the masonry structure changes from pattern O to pattern H. Also, during holding step, there are no joints opening. However, during unloading, some head joints near the fixed rigid plate reopen and, as a result, there is a change to pattern O. In both cases (x and y-LHU), as expected, joints near the moving rigid plates close before those near the fixed rigid plate. Regarding time variations of stresses in the masonry structures, it has been noticed that, under the same applied displacement load (for both x and y-LHU), macroscopic stresses in the x direction (Σ xx ) are higher than stresses in the y direction (Σ yy ). This behaviour clearly illustrates the macroscopic orthotropic behaviour of the structure. In the It should be noted that, the peak in Σ xx is much higher as compared to the peak in Σ yy . Also, the time period during which Σ yy is almost zero is longer as compared to this of Σ xx . This can attributed to that the number of head joints is less as compared to the number of bed joints (4 head joints and 14 bed joints, see figure 16). Therefore, the closure of all head joints (due to gradual increase of applied x displacement) and the increase of macroscopic stiffness in the x direction of the material occurs faster as compared to the closure of all bed joints (due to gradual increase of applied y displacement) and the increase of the material stiffness in the y direction.

Schematics of masonry structure geometry subjected to constant stress uniaxial compression in the x and y directions as well as the boundary conditions are presented in figure 16. Time variations of homogenized total strains in masonry structures subjected to constant uniaxial stress LHU conditions in the x, first loading case, and y, second loading case, directions are presented in figure 21. In both cases, the total strains in the normal directions to the surfaces of the dry joints are reported (i.e. E t xx and E t yy ). The loading steps are: first, applying constant stress of 1.5 MPa, second, holding for 50, 000 seconds and finally, unloading. Overall, it has been observed that total strains, in the loading direction, increases instantaneously due to load application, then increases at constant rate due to creep and finally it decreases slightly due to unloading. In addition, after load removal, the recovered strain (in the loading direction) is very small as compared to the instantaneous strain (due to the applied load at t = 0 sec). This can be attributed to that, first, during unloading, only few joints reopen and their final thickness is very small as compared to the initial joint thickness, second, the permanent deformation resulting from the viscoplastic behaviour of the structure.

Moreover, in the first loading case (x-LHU), the resulting instantaneous strain, in the loading direction (E t xx at t = 0 seconds), is less as compared to the resulting instantaneous strain, in the loading direction (E t yy at t = 0 seconds), in the second loading case (y-LHU).

The reason for this is that the number of bed joints is higher than the number of head joints. Consequently, under the same applied stress level, all head joints close faster than bed joints and the increase of the macroscopic stiffness of the masonry in the x direction (during load application step) occurs faster as compared to the increase of the structure stiffness in the y direction. This leads to lower values of resulting instantaneous strain under the same applied stress level. joints are closed and they remain closed until the end of the holding step. After unloading, some joints reopen and the recovered strain in the y direction is very small as compared to the instantaneous strain. In addition, there is no strain recovery in the x direction, because it is mainly caused by the viscoplasticity of the structure.

Conclusion

In the present study, computational homogenization of elastic viscoplastic refractory masonry with dry joints is performed. Several finite element simulations of mortarless masonry walls subjected to constant stress and constant strain uniaxial loading, holding and unloading are performed. It has been shown that the macroscopic elastic viscoplastic behaviour of mortarless masonry walls is orthotropic and nonlinear due to joints closure, reopening, creep and stress relaxation.

In addition, after unloading only few reopen and the recovered strain is very small as compared to the instantaneous applied strain due to the viscoplastic behaviour of the structure and the final small joints thickness after unloading.

The presented numerical model takes into account the orthotropic, nonlinear, elastic viscoplastic behaviour of refractory masonry structures at high temperature and the presence of dry joints. It can be used to model and predict the homogeneous mechanical response of many high-temperature large sized industrial applications such as steel ladles and furnaces.

  simplicity, cuboid bricks are considered in the present study. Bricks with height h b , width l b and depth d b are periodically arranged in running bond texture. Dry joints with initial thickness (g 0 << h b , l b , d b ) are separating the bricks from each other. The dimensions of the brick are: l b = 124 mm, h b = 76 mm, d b = 185 mm, and the initial thickness of the joints is 0.1 mm. Often, these joints are resulting from the surface unevenness, shape and dimensions errors of the bricks during production. However sometimes, for instance rotary kiln and blast furnace, the joints are designed by attaching cardboard spacers to the bricks during the installation to compensate for thermal expansion effects.

Figure 1 :

 1 Figure 1: (a) Steel ladle lined with mortarless refractory masonry [33]. (b) Schematic of mortarless masonry wall showing the periodically arranged bricks and the gaps between them. Vertical gaps are called head joints and horizontal gaps are called bed joints.

2 .

 2 Pattern B: bed joints are closed, while head joints are open. 3. Pattern H: head joints are closed, while bed joints are open. 4. Pattern O: head and bed joints are open (all joints are open).

Figure 2 :

 2 Figure 2: Schematics of possible joint patterns of mortarless refractory masonry structure (a) pattern O, (b) pattern H, (c) pattern B, and (d) pattern C. The four joint patterns are defined based on the state of bed and head joints (open or closed).

Figure 3 :

 3 Figure 3: (a) Schematic of periodic mortarless masonry structure in pattern B. (b) Periodic unit cell selected from the structure and used for the nonlinear finite element based homogenization analysis.

Figure 4 :

 4 Figure 4: Von Mises stress (in MPa) distribution in a periodic unit cell subjected to periodic boundary conditions and (a) uniaxial tension along the x-direction (σ xx =1.5 MPa), (b) uniaxial tension along the y-direction (σ yy =1.5 MPa) and (c) simple shear in the xy plane (τ xy =1.5 MPa).

Figure 5 :

 5 Figure 5: Displacement discontinuity in a periodic unit cell subjected to uniaxial tension along the x direction. The displacement discontinuity is caused by the presence of open head joints (vertical gaps).

  behaviour of joint pattern B. Time variations of the macroscopic total strains (instantaneous elastic and viscoplastic) in a periodic UC subjected to uniaxial and simple shear constant stress loading conditions (σ ii = τ ij = 1.5 MPa) are shown in figure 6. It can be seen that the macroscopic elastic viscoplastic behaviour of joint pattern B is orthotropic (due to the presence of joints) and the homogenized mechanical responses of the material in the y and z directions are similar (as there are no joints in the z direction, and joints in the y direction are closed). In order to describe the macroscopic viscoplastic behaviour of joint pattern B, an orthotropic creep equation is required. It should be noted that equation 4 is suitable only for describing the local creep behaviour of isotropic materials such as joint pattern C.

Figure 6 :

 6 Figure 6: Macroscopic total strains in a periodic unit cell subjected to periodic boundary conditions as well as constant stress uniaxial tension (in the x, y and z directions) and simple shear loading (in the xy, xz and yz planes) (Σ ii = Σ ij =1.5 MPa). The macroscopic elastic viscoplastic behaviour of the structure is orthotropic.

  other words, a head joint can reopen if σ xx > 0 and a bed joint can reopen if σ yy > 0. As the macroscopic stresses are linearly dependent on local stresses (see equation 10), joints reopening criterion can be rewritten in terms of macroscopic stresses using localization tensor Σ yy -→ σ yy and Σ xx -→ σ xx . The joints closure, reopening and transition criteria from joint pattern O to pattern H, joint pattern H to pattern B, joint pattern H to pattern C, joint pattern B to pattern C, and the opposite are presented in figure 7.

Figure 7 :

 7 Figure 7: Schematics of all possible joint patterns of mortarless refractory masonry structure and joints closure and reopening criteria due to loading and unloading. Joints closure criteria are written in terms of macroscopic total strains. Joints reopening criteria are defined in terms of macroscopic stresses.

Figure 8 :

 8 Figure 8: Comparisons between average total strains obtained from the detailed micro model (bricks and joints are considered) and macro model (bricks and joints are replaced by equivalent homogeneous medium) for constant stress uniaxial tension (in the x, y and z directions) and simple shear loadings (in the xy, xz and yz planes).

  further examine the macro model and the identification technique, comparisons between the results of the micro and macro model under two different loading conditions (from those used in the identification) are carried out. The two loading conditions are biaxial constant stress loading and constant strain rate uniaxial loading. In the biaxial constant stress loading case, the following combination of compressive and tensile stresses are considered: compressive stress of 3 MPa in the x direction and tensile stress of 1 MPa in the y direction (Σ xx = -3 M P a and Σ yy = 1 M P a). Comparisons between average total strains obtained from the micro and macro models are presented in figure9. Good agreements between the micro and macro models are observed.

Figure 9 :

 9 Figure 9: Comparisons between average total strains obtained from detailed micro model (bricks and joints are considered) and macro model (bricks and joints are replaced by equivalent homogeneous medium) for constant stress biaxial loading (compression in the x direction and tension in the y direction).

Figure 10 :

 10 Figure 10: Comparisons between average macro stresses obtained from detailed micro model (bricks and joints are considered) and macro model (bricks and joints are replaced by equivalent homogeneous medium) for constant strain rate uniaxial tension (in the x, y and z directions) and simple shear loadings (in the xy, xz and yz planes).

Figure 11 :

 11 Figure 11: Setup of the biaxial compression tests a) at room temperature and b) 1200 °C of Magnesia Chromite refractory masonry wall [14, 56-58].

  at high temperature and viscoplasticity. Joint patterns by the end of load application step of the two tests are presented in figure 15. Higher percentage of closed bed and head joints are noticed in the case of the test at 1200 °C as compared to room temperature test.

Figure 12 :

 12 Figure 12: Stress-strain diagrams of refractory masonry wall subjected to biaxial compression load at a) room temperature and b) 1200 °C: experimental and numerical results.

Figure 13 :

 13 Figure 13: Horizontal displacement fields, in the x direction, in refractory masonry wall subjected to biaxial compression load at a) room temperature and b) 1200 °C. The word fixed indicates the location of the fixed rigid plates, while the arrows indicate the location and load direction of the moving rigid plates. The interactions between the contact surfaces of the wall and the fixed and moving rigid plates are considered using frictional contact.

Figure 14 :

 14 Figure 14: Vertical displacement fields, in the y direction, in refractory masonry wall subjected to biaxial compression load at a) room temperature and b) 1200 °C. The word fixed indicates the location of the fixed rigid plates, while the arrows indicate the location and load direction of the moving rigid plates. The interactions between the contact surfaces of the wall and the fixed and moving rigid plates are considered using frictional contact.

Figure 15 :

 15 Figure 15: Joints pattern in a refractory masonry wall subjected to biaxial compression load at a) room temperature and b) 1200 °C. The word fixed indicates the location of the fixed rigid plates, while the arrows indicate the location and load direction of the moving rigid plates. The interactions between the contact surfaces of the wall and the fixed and moving rigid plates are considered using frictional contact.

Figure 16 :

 16 Figure 16: Schematics of masonry structures subjected to uniaxial compression in a) y-direction and b) x-direction: geometry and boundary conditions.

Figure 17 :

 17 Figure 17: Gradual closure and reopening of joints in masonry wall subjected to uniaxial compression loading (displacement), holding and unloading in the y and x directions. (a) All joints are open (pattern O) at time = 0 seconds. (b) Gradual closure of joints due to the applied displacement. (c) All bed or head joints remain closed during the holding step. (d) Some joints reopen due to unloading.

Figure 18 :

 18 Figure 18: Time variations of a) global displacement and b) stress of masonry wall subjected to uniaxial compression displacement loading, holding and unloading in the x and the y directions.In both cases (x and y), the applied displacement load is equal. However, resulting stresses are different due to the orthotropic elastic-viscoplastic behaviour of the masonry structure.

Figure 19 :

 19 Figure19: Contours of viscoplastic strain and final head joints thickness (in mm) after unloading in mortarless masonry structure subjected to constant uniaxial compression displacement loading, holding and unloading in the x direction. After unloading, few head joints reopen (g head > 0) and their thickness is very small as compared to the initial thickness (0.1 mm).

Figure 20 :

 20 Figure20: Contours of viscoplastic strain and final bed joints thickness (in mm) after unloading, in mortarless masonry structure subjected to constant uniaxial compression displacement loading, holding and unloading in the y direction. After unloading, few bed joints reopen (g bed > 0) and their thickness is very small as compared to the initial thickness (0.1 mm).

Figure 21 :

 21 Figure21: Time variations of macroscopic total strains in masonry structure subjected to constant stress loading, holding and unloading in the x and y directions. The total strains in the normal directions to the surfaces of the dry joints are reported. In both cases (x and y), the applied stress load is equal (1.5 MPa). However, resulting total strains are different due to the orthotropic elastic-viscoplastic behaviour of the masonry structure.

Figure 22 :

 22 Figure 22: Gradual closure and reopening of joints in masonry wall subjected to 1.5 MPa constant stress LHU in the y direction. Points a to e are indicated in figure 21. (a) Gradual closure of bed joints due to the applied load. (b) All bed joints are closed by the end of the loading step. (c) Head joints close gradually due to creep. (d) All bed and head joints are closed. (e) Joints reopening due to unloading.

  Four possible joint patterns are predefined based on the state of bed and head joints. Each joint pattern has different homogeneous elastic viscoplastic response which has been identified using nonlinear homogenization technique. The homogenized elastic viscoplastic response of the first three patterns (patterns pattern O, pattern H and pattern B) is orthotropic due to the presence of open joints whereas, the homogenized response of the last joint pattern (pattern C) is isotropic and similar to the mechanical response of the base material of the bricks because all joints are closed. It has been shown that the presence of open joints leads to a significant reduction in the effective elastic stiffness of the structure and orthotropic viscoplastic behaviour. Suitable macroscopic orthotropic constitutive laws are used to describe the homogenized elastic viscoplastic behaviour of the joint patterns. To assess the accuracy of the developed macro model, several finite element simulations of uniaxial, simple shear and biaxial loadings are carried out. Then, comparisons between average mechanical responses from the detailed micro model and the macro model are performed. Good agreements between the results from the micro and macro models are observed. During operation, the masonry structure changes from one joint pattern to another leading to a change in the macroscopic elastic viscoplastic behaviour of the structure. This change is taken into account by employing suitable joints closure, reopening and pattern transition criteria. The joints closure and reopening criteria are written in terms of macroscopic stresses and strains.

  Effective viscoplastic 4 th order tensor of the four joint patterns of Magnesia Chromite refractory masonry wall at different temperatures.

  : nonlinear homogenization, computational homogenization, refractories,

	Ė vp 2 nd order macroscopic steady state	g head Instantaneous thickness of head
		viscoplastic strain rate		joints
	E t	2 nd order macroscopic total strain	h b	Brick height
	masonry, creep, modelling tensor	l b	Brick width
	Nomenclature I 2 nd order identity tensor	n	Creep stress exponent
	E	Macroscopic or effective Young's	T	Temperature
	Subscripts modulus	ε t	2 nd order local total strain tensor
	b G	Brick Macroscopic or effective shear mod-	ε e V UC Periodic unit cell volume 2 nd order local elastic strain tensor
	ulus bed Bed joint d A Creep power law multiplier Deviatoric eq Equivalent d b Brick depth	ε vp 2 nd order local viscoplastic strain Y Young's modulus of the bulk mate-rial tensor Σ eq Macroscopic equivalent stress Ėvp Macroscopic equivalent viscoplastic eq strain rate
	E ij head Head joint Components of macroscopic strain	σ eq Superscripts Local equivalent stress
	tensor i, j, k, l Indices of the dimension g 0 Initial thickness of dry joints	Σ ij e	Components of macroscopic stress tensor Elastic
	U C Periodic unit cell g bed Instantaneous thickness of bed	σ ij t	Components of local stress tensor Total
	x, y, z Coordinates joints	ν vp	Macroscopic or effective Poisson's Viscoplastic
	Greek Letters		ratio
	ν	Poisson's ratio of the bulk material	Roman Letters
	ε vp 2 nd order local viscoplastic strain	u	Displacement vector
		rate tensor	X	Position vector
	Σ	2 nd order macroscopic stress tensor	C e 4 th order macroscopic or effective
	σ	2 nd order local stress tensor		elastic stiffness tensor
	σ d	2 nd order local deviatoric stress ten-	N	4 th order viscoplastic localization
		sor		tensor

Table 1 :

 1 Material parameters of bauxite refractory at T = 1000 °C[START_REF] Blond | Modelling of high temperature asymmetric creep behavior of ceramics[END_REF].

	Property

  1ν 23 ν 32 E 22 E 33 ν 21 + ν 31 ν 23 E 22 E 33 ν 31 + ν 21 ν 32 E 22 E 33 12 + ν 13 ν 32 E 11 E 33 1ν 31 ν 13 E 11 E 33 ν 32 + ν 31 ν 12 E 11 E 33

	0	0	0
	0	0	0

ν ν 13 + ν 12 ν 23 E 11 E 22 ν 23 + ν 13 ν 21 E 11 E 22 1ν 12 ν 21 E 11 E 22

Table 2 :

 2 Effective elastic and viscoplastic 4 th order tensors of the four joint patterns.

	Pattern		C (MPa)					N (-)
	1	0	0	0	0	0	0	0	0	0	0	0
		1	0	0	0	0		0	0	0	0	0
	O		11000	0 0	0 0	0 0			2.01	0 0	0 0	0 0
		sym.			4583	0		sym.			6.01	0
						4583						6.01
	11485	0	2291	0	0	0	2.01	0	-0.99	0	0	0
		1	0	0	0	0		0	0	0	0	0
	H		11485	0 0	0 0	0 0			2.01	0 0	0 0	0 0
		sym.			4583	0		sym.			6.01	0
						4583						6.01
	4430 1113 1103	0	0	0	6.25 -1.10 -0.49	0	0	0
		11738 2560	0	0	0		2.01 -0.99	0	0	0
	B		11733	0 3602	0 0	0 0			1.99	0 7.18	0 0	0 0
		sym.			3315	0		sym.			8.25	0
						4583						6.01
	12222 3055 3055	0	0	0	2.01 -0.99 -0.99	0	0	0
		12222 3055	0	0	0		2.01 -0.99	0	0	0
	C		12222	0 4583	0 0	0 0			2.01	0 6.01	0 0	0 0
		sym.			4583	0		sym.			6.01	0
						4583						6.01

  In addition, the structure may change from pattern O to pattern C if both head and bed joints close. Based on the instantaneous thickness of the joints g, head and bed joints are considered to be open as:

	Joint open	if	g bed or head > 0	(28)

1 to 0.2 mm) to zero and the structure changes from pattern O to either pattern H (when head joints close) or pattern B (when bed joints close).

  the case of uniaxial loading in the y and z directions, values of the macroscopic total strains (E t yy and E t zz ) are almost the same because there are neither open joints in the y direction nor in the z direction. However, E t yy and E t zz are less than E t xx , because the presence of open joints in the x direction leads to a decrease in the effective stiffness of the structure and compressibility in the respective direction.

  of Magnesia Chromite refractory masonry walls subjected to biaxial compression load at room temperature and 1200 °C are carried out. The biaxial compression test setup is shown in figure 11. Magnesia Chromite bricks with length of 124 mm, height of 76 mm and depth of 185 mm are periodically arranged in a running bond

texture. The mechanical parameters of the constitutive material of the bricks are reported in table

3

. The size of the wall is 1100 × 1100 × 185 mm 3 and it is made up of 15 courses, all of them are active. The wall is surrounded by four ceramic plates, two of them are fixed while, the other two can move thanks to two hydraulic jacks.

Table 3 :

 3 Mechanical properties of Magnesia Chromite bricks at different temperatures[START_REF] Andreev | Compressive behaviour of dry joints in refractory ceramic masonry[END_REF][START_REF] Triessnig | Determination of Thermo-Mechanical Properties at the RHI Refractories Technology Center Leoben[END_REF][START_REF] Prietl | Determination of material-specific parameters of refractory materials and linings under uniaxial and biaxial load conditions for the non-ferrous metal industry[END_REF][START_REF] Jin | Compressive creep testing of refractories at elevated loads-Device, material law and evaluation techniques[END_REF]].

	Temperature (°C)	Y (GPa)	ν (-)	n (-)	A (S -1 M P a -n )
	25	28.5	0.2		---
	800	15.4	0.2		---
	1000	10.6	0.2		---
	1100	9.8	0.2	2.86	1.18 × 10 -16
	1200	9.1	0.2	2.86	2.77 × 10 -16

  yy is very small as during this period bed and head joints are open (material stiffness is very small). Then, after bed joints closure, Σ yy increases at higher rate to reach a peak of 2.5 MPa. Finally, a stress decay has been noticed due to stress relaxation.

case of x-LHU, during the first 75 seconds Σ xx values are very small, because during this period gradual closure of head joints is occurring and the masonry structure is changing from pattern O to pattern B (i.e. material stiffness is very low). After closure of joints Σ xx increases at higher rate to reach a peak of around 7.5 MPa by the end of loading step. Then during holding, a decrease in Σ xx can be noticed due to stress relaxation and, finally, after unloading Σ xx returns to zero. Similarly, in the case of y-LHU, at the first 170 seconds Σ

Acknowledgment

The authors are grateful to the Centre de Calcul Scientifique en région Centre Val de Loire (France) for the computation time and technical support.

Funding

This work was supported by the funding scheme of the European Commission, Marie Sk lodowska Curie Actions Innovative Training Networks in the frame of the project ATHOR -Advanced THermomechanical mOdelling of Refractory linings -764987 Grant.

Conflicts of interest

None to declare.

Appendix

Table A1:

Effective elastic 4 th order tensor of the four joint patterns of Magnesia Chromite refractory masonry wall at different temperatures.

Pattern