Influence of Boundary Conditions on Numerical Homogenization of High Performance Concrete

Concrete is the most widely used construction material nowadays. We are concerned with the computational modelling and laboratory testing of high-performance concrete (HPC). The idea of HPC is to enhance the functionality and sustainability of normal concrete, especially by its greater ductility as well as higher compressive, tensile, and flexural strengths. In this paper, the influence of three types (linear displacement, uniform traction, and periodic) of boundary conditions used in numerical homogenization on the calculated values of HPC properties is determined and compared with experimental data. We take into account the softening behavior of HPC due to the development of damage (micro-cracks), which finally leads to failure. The results of numerical simulations of the HPC samples were obtained by using the Abaqus package that we supplemented with our in-house finite element method (FEM) computer programs written in Python and the homogenization toolbox Homtools. This has allowed us to better account for the nonlinear response of concrete. In studying the microstructure of HPC, we considered a two-dimensional representative volume element using the finite element method. Because of the random character of the arrangement of concrete’s components, we utilized a stochastic method to generate the representative volume element (RVE) structure. Different constitutive models were used for the components of HPC: quartz sand—linear elastic, steel fibers—ideal elastic-plastic, and cement matrix—concrete damage plasticity. The numerical results obtained are compared with our own experimental data and those from the literature, and a good agreement can be observed.


Introduction
Significant progress in developing high-performance concrete (HPC) has been observed since the invention of a concrete from reactive powders (RPC) in the 1990s [1][2][3]. Concretes of this type are characterized by a much higher compressive strength, flexural tensile strength, ductility, and resilience to an aggressive environment than traditional concrete. The use of such concretes in the construction industry allows one to build challenging engineering structures (bridges and towers) and high-rising buildings. Compressive and tensile strengths and modulus of elasticity are the main mechanical parameters that affect the design properties of HPC. An accurate determination of these parameters is very important from a practical viewpoint of engineering and design for efficient and safe HPC constructions. Due to the complex composition of high-strength concrete and therefore complexity of its microstructure, which is finally formed during the hydration process, determining the abovementioned parameters is not easy [4]. Publications about the properties of cement-based composites, a name that can be also used for HPC, are concerned in most cases with experimental testing. They contain mainly parametric analysis on the influence of the volumetric fraction of individual components on the HPC response, analysis of different orientations and shape factor of fibres [5][6][7][8][9][10][11], or analysis of different types and shapes of fibres [12][13][14][15][16]. The material parameters we observe at the macro-scale are a result of phenomena taking place in the microstructure of the material, i.e. on meso-, micro-, and nano-scales. A thorough analysis of these phenomena also gives an opportunity to trace different aspects of the deformation process of HPC. One of the basic phenomena that significantly affect endurance of HPC is the formation and propagation of cracks in the microstructure of material [17][18][19].
There are many methods used to describe effective properties of heterogeneous materials, of which the most intuitive is direct homogenization [20]. It requires an estimate of the value of the averages within the area occupied by the representative volume element (RVE). The concept of RVE was defined and used by Hashin [21], where both elastic and viscoelastic properties as well as failure criteria for fibre composites where considered; by Hill [22] for elastic media; and more recently by Ostoja-Starzewski [23] in the context of plasticity of random media. The averaging operation is preceded by a solution of the boundary value problem defined for RVE with homogeneous boundary conditions. Unlike direct methods, indirect homogenization does not require estimation of the average values in the RVE cell. Indirect methods should be understood as a whole group of techniques based on Eshelby's theory [24]. They were used in designing composite materials. Many examples of practical applications of these methods can be found in [21], including the self-consistent scheme [25], which is developed for prediction of elastic moduli of two-phase composites assuming that any (ellipsoidal) inclusion is embedded in a homogeneous medium and subject to boundary conditions at infinity. Another way to describe heterogeneity caused by the introduction of a discrete phase into the homogeneous base, i.e., matrix, is characterized by the Mori-Tanaka method [26]. Another technique to be mentioned is differential indirect method [20]. It consists of a volumetric balance matrix and discrete phase participation. A large group of homogenization methods is based on the minimum principle of potential energy, and they belong to the so-called variational homogenization methods. With their help, the lower or upper limits of parameter values are often effectively estimated [20].
Relatively recent, introduced at the turn of the 1970s and the 1980s by Bensoussan et al. [27] and Sanchez-Palencia [28], a method of asymptotic expansion for multiple scales was used [29]. It is based on a mathematical theory of homogenization, in which disturbances of unknown functions are assumed to be so small that they can be presented as an asymptotic expansion with respect to a small parameter, which is a characteristic dimension of the microstructure. This method was used to analyse materials that are physically nonlinear [30,31].
For analyzing inhomogeneous materials where the geometry and distribution of microstructure components cannot be precisely known, stochastic methods [32][33][34] can be used that also provide an estimation of properties of a surrogate material model on a macro-scale.
The last group of homogenization methods to be mentioned here are numerical methods, presented, e.g., in [35][36][37][38] and in particular for cement pastes in [39,40]. They are free from the restrictions that are imposed in strictly analytical techniques, for example, limiting the shape of a particle to the sphere or ellipse or being unable to analyze problems that are geometrically and/or physically nonlinear. In recent years, together with the increase in computing power of computers, the methods of numerical homogenization have flourished. These methods make use of computational algorithms aimed at estimating numerically the relationship between strain and stress at a macroscopic material point based on separate calculations carried out on a representative RVE that is assigned to the point and reflects the material's microstructure there. Calculations can be carried out using various numerical techniques, of which the finite element method is the most popular one.
Modelling and computation of (high performance) concrete, owing to its complex multi-phase microstructure [5,[41][42][43] and exhibiting irreversible and brittle behavior (deterioration and fracture), is an extremely difficult problem. From a computational point of view, the difficulties are manifested when solving a boundary value problem for concrete specimens or concrete structural members in the strain-softening regime, where well-posedness of the problem is lost and the pathological mesh-dependence of a solution by the finite element method occurs. In order to overcome the computational difficulties and to avoid spurious numerical results, some forms of regularization of the boundary value problem corresponding to nonlinear constitutive laws for concrete have been developed. The regularization can be achieved by incorporating into the constitutive law some nonlocal information about the deformation process from a vicinity of the material point, either by accounting for gradients of selected quantities (damage and plastic strain) or by using a phase-field approach; see Bažant et al. [44], de Borst et al. [45], Ramm et al. [46], smf Kaliske et al. [47], for example.
In this paper, we follow the procedure of a numerical two-scale homogenization approach to model the behavior of high-performance concrete, paying special attention to the boundary conditions imposed on the RVE that embodies the concrete's microstructure. In fact, we treat the HPC and its microstructure as a composite consisting of four (or six) phases: cement matrix, fine (and thick) quartz sand, steel micro fibres (and steel fibres), and air voids. The following constitutive models were used for the components of HPC: quartz sand-linear elastic, steel fibers-ideal elastic-plastic, and cement matrix-concrete damage plasticity (CDP). Thus, the interphase and interfacial phenomena are not treated explicitly, being per assumption accounted for implicitly by the CDP model of cement matrix. In two-scale numerical homogenization, the macro-scale and micro-scale are distinguished, and at each of the scales, the corresponding boundary value problem (BVP) is solved using the finite element method. The boundary conditions imposed on the solution of the macro-scale BVP are a direct result of the support constraints of a body (structural element) and therefore are uniquely defined, whereas the boundary conditions imposed on the solution of the micro-scale BVP defined on RVE may be assumed in a number of ways. We solve the micro-scale BVP using three types of boundary conditions imposed on the boundary of RVE: uniform displacement boundary conditions (DBC), traction boundary conditions (TBC), and periodic boundary conditions (PBC).
The influence of boundary conditions used in the numerical homogenization of nonhomogeneous media on the obtained values of material parameters and the stability of computational schemes have been studied by many researchers, mainly for elastic composites, e.g., in [48][49][50]. A mixed formulation of the BVP on the micro-scale was developed in [48] using weakly periodic boundary conditions, which allowed a part of microstructure to be completely cut loose by cracks and in which the idea of an elastic cohesive zone was then used as a regularization method. In order to remove rigid body motions of the elastic RVE with TBC, the concept of semi-Dirichlet boundary conditions was introduced in [49], the latter enforcing non-homogeneous displacements at selected points simultaneously satisfying the Neumann-type conditions. Numerical tests for elastic woven composites have revealed high sensitivities of the in-plane extensional modulus and Poisson's ratio to the type of boundary conditions, and a mix of TBC and PBC proved to be best in representing an experimental strain field [50].
This contribution is a continuation of our papers [51][52][53][54], where a linear-elastic twoscale model of HPC and its experimental validation were presented. Herein, we extend our previous analyses to the case of nonlinear elasto-plastic behavior with damage. Our aim is to examine the impact of the particular type of boundary conditions on the response of two-dimensional RVE of high-performance concrete.

Research Methodology
The results of the numerical simulations described herein relate to the experimental findings obtained in our own laboratory tests [53] and those contained in [55]. We restrict our numerical computations to the local BVP on the RVE of high-performance concrete, checking the influence of different boundary conditions imposed on the boundary of RVE on the simulated response of the tested HPC. Accounting for randomness of the RVE microstructure, we generated positions of different phases, each in the amount as in the mix recipes in Table 1, by a pseudo-random number generator (for details, see [51]). The solution of the elasto-plastic damage BVP is solved by the finite element software package, Simulia Abaqus [56], supplemented with our in-house written script in Python and the Homtools package [57].

Microstructure of Modeled HPC
The microstructure of high-performance concrete has been studied by many authors, and the topic is well-recognized in the world literature, e.g., [5,[41][42][43]58]. In this contribution, the microstructure of HPC is defined by means of a two-dimensional representative volume element (RVE) and the finite element method is used for modelling. The idea adopted for the numerical finite element analysis of the RVE microstructure is shown in Figure 1. The method with which we generated the considered RVE was described in [51,54]. The actual geometry of the HPC ingredients was approximated by square finite elements of dimensions 0.2 × 0.2 mm of the amount given in the concrete mix proportions table (Table 1). The size of the entire RVE was 10 × 10 mm. The RVE consists of several main micro-ingredients: OS 36 and OS 38-fine quartz sand, OS 30-thick quartz sand, steel (micro) fibres, and cement matrix. The constitutive models used for the components of HPC were quartz sand-linear elastic, steel fibersideal elastic-plastic, and cement matrix-concrete damage plasticity (Section 2.4). The RVEs of mixtures I and II were generated for numerical homogenization (Figure 2). air voids (yellow): empty space (no finite elements), 4 wt%.

Parameters of Concrete Damage Plasticity (CDP) Model
For modelling the cement-based matrix, i.e., a hardened mixture of cement, silica fume, water, and superplasticizer, we adapted the concrete damage plasticity (CDP) model as described in the influential papers [59][60][61] and Abaqus documentation [56]. Let us recall that the CDP model captures the nonlinear behavior of concrete by accounting for simultaneous development of permanent (plastic) deformation and elastic stiffness degradation using two damage variables, one for tensile damage (d t ) and the other for compressive damage (d c ), to account for different damage responses of concrete in tension and in compression.
The plastic-damage model [56] assumes nonassociated potential flow with the Drucker-Prager hyperbolic function that is continuous and smooth and asymptotically approaches the linear Drucker-Prager yield condition [62]. The evolution of the elastic-plastic-damage deformation process is complex from the numerical viewpoint because it is described by inequality relations and Kuhn-Tucker complementarity conditions (loading/unloading conditions), e.g., [61,63].
The parameters of the CDP model and their values used in the calculations are gathered in Table 2, wherein β-the internal friction angle of concrete. In the CDP model, β is defined as the inclination angle of the Drucker-Prager surface asymptote to hydrostatic axis of the meridional plane; m-eccentricity of the surface of the plastic potential. This is the distance measured along the hydrostatic axis between the apex of the Drucker-Prager hyperbola and the intersection of the asymptote of this hyperbola, calculated in practice as a ratio of tensile strength to strength for compression; f b0 / f c0 -number specifying the compressive strength ratio in a two-axis state for the strength in a single-axis state; K c -parameter defining the shape of the surface of the plastic potential on a deviatoric plane; η-viscoplasticity parameter, used to regularize the concrete constitutive equations.
The used stress-strain relationships for the cement matrix in compression and in tension are shown in Figures 3 and 4, respectively, with the compressive strength equal to 200 MPa and the tensile strength equal to 20 MPa. These relations are described with the function proposed by Saenz [64] as where and The lower index γ ∈ {c, t} stands for a type of stress, with c standing for compression and t standing for tension, respectively. Furthermore, ε γu is the ultimate strain in the matrix and f γu is its corresponding stress, whereas f γm is the extremal stress sustained by the matrix with its corresponding strain ε γ1 . In the calculations, we assumed for compression E = 55 GPa, f cm = 200 MPa, f cu = 196 MPa, ε cu = 0.0028, and ε c1 = 0.0025 and for tension f tm = 20 MPa, f tu = 19 MPa, ε tu = 0.00028, and ε t1 = 0.00025. The values of the parameters of Equation (1) calculated by Equation (2) for compression and tension are collected in Table 3.   The evolution of the compression damage scalar d c for the cement matrix as a function of strain and that of the tension damage scalar d t are illustrated in Figures 5 and 6 and Equation (4), respectively.

Numerical Homogenization and Boundary Conditions
In the method of two-scale numerical homogenization, the response of a material at the macro-scale is determined by an analysis of the behavior of the material's internal structure at the micro-scale. In the case of nonlinear material behavior, achieving the equilibrium at each scale together with the compatibility of information between the two-scales requires multiple exchanges of needed information between the scales after solving a corresponding nonlinear boundary value problems formulated at each of the scales. On the micro-scale level, the distributions of micro-stresses and micro-strains were calculated, which via homogenization provide the needed information of averaged macroscopic quantities to the macro-scale. In this contribution, our analysis was restricted to the solution of the local nonlinear BVP defined on the two-dimensional representative volume element (RVE). When the characteristic microscopic length was one order smaller than the characteristic macroscopic length, we could take into consideration only effects of the first order.
The analysis was carried out within the realm of linear kinematics. The notation that a bar above a symbol denotes a macroscopic variable or quantity was used. Let x denote the vector of coordinates of a material point at the macroscopic level, and let x stand for coordinates of points within the RVE of volume V defined in the material point. Further, let us assume that the microscopic displacement field u can be additively split into a linear part εx and a fluctuating part r: where ε(x, t) is the macro-strain tensor and t is a time-like parameter. At down-scaling (macro-micro transition) and solving the local BVP at RVE, the strain tensor ε is treated as a known quantity (data). At up-scaling (micro-macro transition), the elements of the macro-strain tensor ε can be defined as mean values of the corresponding micro-strains ε averaged over the RVE: where n is the unit outward normal of the boundary of RVE, Γ = ∂V, and ⊗ is a tensor product of vectors. The Formula (6) shows that the macroscopic strain tensor can be expressed by micro-displacements u at the boundary of the RVE. However, Equation (6) is valid provided that the zero gradient condition of the microscopic displacement fluctuation field r is satisfied: Fulfilling the above condition ensures that deformation of the RVE boundary in the medium sense is in accordance with the pre-set macro-strain ε.
Macro-stresses can be defined, similar to macro-strains, as mean values of the microstresses σ: This relationship can be derived from Hill's theorem [25], which says that the work done by macro-stresses on the corresponding macro-strains is equal to the mean value of the work performed by micro-stresses on the corresponding micro-strains: where the symbol • stands for averaging over the RVE volume.
The finite element method is applied for numerical calculations, with the finite element designated as the CPS4R Abaqus system library [56]. The RVE area was discretized with 2500 finite elements, each with dimensions of 0.2 × 0.2 mm; see Figure 2.
Numerical tests were carried out enforcing three types of boundary conditions on the boundary of RVE: linear displacement boundary conditions (DBC), uniform traction boundary conditions (TBC), and periodic boundary conditions (PBC). Before starting the numerical simulations of the inhomogeneous nonlinear BVP for high-performance concrete, the work and proper interaction between the Homtools and the Abaqus packages were verified with the example of RVE for a homogeneous material with a hole (Section 3.1).

Linear Displacement Boundary Conditions (DBC)
This method consists of applying on the boundary of RVE the displacement field that would occur if the strain were uniform inside the RVE. For the considered linear kinematics, the boundary conditions can be defined by the formulas in Equations (5) and (7): in which ε(u) is the micro-strain field and x is the position vector of point x on boundary Γ. There is no restriction concerning the use of this method, except that no rigid part must intersect the boundary, and holes are permitted.

Uniform Traction Boundary Conditions (TBC)
This method consists of applying on the boundary of RVE the stress vector field that would occur if the stress were uniform inside the RVE. For the considered linear kinematics, the boundary conditions can be defined as follows: where σ is the Cauchy stress tensor and n denotes the outward unit normal. There is no restriction concerning the use of this method, except that no holes must intersect the boundary.

Periodic Boundary Conditions (PBC)
Enforcing the PBC is theoretically relevant for periodic media, which can be defined by a periodicity cell and the associated periodicity vector of translation. The periodic homogenization process consists of assuming that the strains and stresses are periodic at the level of the periodicity cell (which is defined as the RVE). The periodicity of stresses and strains leads to specific periodic boundary conditions for the localization problem on the RVE. In order to introduce periodic boundary conditions, the boundary of RVE is decomposed into two opposing parts, Γ + and Γ − , such that Γ = Γ + ∪ Γ − . Each point x + on Γ + is associated with a unique point x − on Γ − , and the unit normal vectors at these boundaries satisfy n − = −n + . Then, the PBCs are defined as follows: where t ± = σn ± is the stress vector. For this type of boundary condition, u is forced to be periodic and t is forced to be antiperiodic. Note that LDBCs satisfy the periodicity of u only whereas TBCs satisfy the antiperiodicity of t only. There is no restriction concerning the use of this method; periodic holes and rigid parts intersecting the boundary are permitted.

Test Example: Homogeneous and Linear Elastic RVE with Hole
A shear test of a homogeneous linearly elastic RVE with a hole in the centre was carried out in order to check if the functioning of the Homtools package within the Abaqus calculation environment is correct. This common testing procedure preceded relevant simulations of the complex heterogeneous nonlinear BVP. It allowed us to detect the possible pitfalls in the implementation of the method of numerical homogenization. In the test, the following values of some parameters were assumed: RVE dimensions-10 × 10 mm; hole diameter-5 mm; Young's modulus of linearly elastic material, E = 200 GPa; and Poisson ratio, ν = 0.3. Figure 7 shows the results of the shear test with the imposed macro-strain ε = {0, 0, 1} for different boundary conditions and the corresponding deformation modes of the RVE with the distribution of reduced micro-stresses of Huber-von Mises-Hencky (HMH). The obtained results are in good agreement with similar calculations done by other researchers, e.g., [38]. The solution of homogeneous material does not depend on the boundary condition, so the RVE without a hole was also tested. The homogeneous field of micro-stresses was obtained with the same values for the three different types of the boundary condition. It additionally confirmed the correct implementation of the used numerical homogenization scheme into the Abaqus calculation environment.

Compression Test for Mixtures I and II
A compression test ε = {−1, 0, 0} (Voigt notation) was performed for non-homogeneous RVE with nonlinear materials as a first test. The curves presented in Figure 8 show dependence between macro-strain and macro-stress for the micro-structure model of HPC (mixture I). In all presented results, the following notation was used, Mx_XXX_y, where x indicates the number of mixtures 1 or 2; XXX is the type of boundary condition: DBC, PBC, or TBC; and y is the test class (c-compression, t-tensile, and s-shear).
Analyzing the graphs in Figure 8, one can notice that use of the periodic boundary condition in the deformation range of RVE ε = (1.75 ÷ 3.20) × 10 −3 leads to an upper estimate of the compressive macro-stress. In the case of smaller deformations, the upper estimate was obtained using the displacement boundary condition. The difference between PBC and DBC is not significant; both types of conditions on the RVE boundary give similar results in the analyzed case. However, the traction boundary condition definitely leads to the lowest estimate of the macroscopic value, and under this type of boundary condition, the analyzed RVE of HPC exhibits responses that is characteristic for brittle materials-the plastic part clearly visible for the PBC and DBC does not exist for the TBC. After reaching the deformation of ε = 1.2 × 10 −3 with an imposed traction boundary condition, the rapid brittle failure of micro-structures appears. Our own experimental results were compared with the determined numerical ones to validate our computational model. In the Figure 8, the average compressive strength of six concrete specimens 100 × 100 × 100 mm made from mixture I [53] was marked (experiment_M1_c = 106 MPa). The stress-strain path could not be recorded; therefore, in Figures 8-11, only the value of compressive or tensile strength is depicted by a dotted line. Additionally the stress-strain characteristic of concrete made by the authors of the paper [55] is presented. The recipe for that concrete is twin-like in comparison to ours for mixture I. The correlation between the experimental and the numerical results is satisfying. A compression test ε = {−1, 0, 0} was also performed for the model of micro-structure HPC made from mixture II. The results are shown in Figure 9. In this example, the results obtained with the periodic and the displacement boundary conditions are practically identical to the deformation level of ε = 2.22 × 10 −3 . The upper estimate of macroscopic values is obtained when using displacement boundary conditions in this case. Use of the traction boundary conditions also leads to an lower estimate of the macro-value. Six concrete specimens 100 × 100 × 100 mm made from mixture II [53] were also examined ex-perimentally. The average compressive strength (marked in Figure 9 as experiment_M2_c) is equal 141 MPa, which is in very good agreement with our own numerical result, similar to the experimental results taken from the paper [55] (marked in Figure 9).

Tensile Test for Mixture I and II
The next test performed for the non-homogeneous RVE with material nonlinearity was a tensile test ε = {1, 0, 0}. The curves presented in Figure 10 depict the obtained relation between macro-stress and macro-strain for HPC modelled as RVE-mixture I and, in Figure 11, the similar response for RVE-mixture II. Analogical to that in the compression test, the upper estimate of macro-stresses is provided for mixture I when using the periodic boundary condition and, for mixture II, when using displacement boundary condition.
The RVE deformation range in the tensile test is an order of magnitude lower than that in the compression test, which is an obvious effect of the assumed values of tensile properties for the CDP model of cement matrix. A detailed description of this phenomenon is outside the present research, but it is a very important aspect of HPC, which deserves a separate study. The authors intend to investigate the tensile response of HPC in the future. The average of the experimentally determined tensile strengths [53] is also presented in Figures 10 and 11. The tensile strength was determined in the flexural test for both mixture I-experiment_M1_t = 12.5 MPa-and mixture II-experiment_M2_t = 18.0 MPa. The agreement of the experimental findings with the numerical results is very good. Figure 12 illustrates some distributions of the compressive and the tensile damage variables. The green areas indicate the places determined numerically where cracks appeared in the HCP micro-structures. A continuity of the micro-structure is breached, and a distribution of the micro-cracks is random. The latter is due to the randomly generated structure of RVE, which reflects the real conditions of HCP manufacturing. The maps of damage variable distribution shown in Figure 12

Shear Test for Mixture I and II
The third performed numerical simulation represents a shear test ε = {0, 0, 1}. The results of the test are shown in Figures 13 and 14. The upper estimate of macroscopic shear stress in the shear test is provided for both cases of the RVE structure (mixture I and II) by displacement boundary conditions. The lower estimate is reached similar to that in previous cases with traction boundary conditions. The shear strength is not experimentally researched by the authors. Figures 13 and 14 contain this parameter taken from [55], where the shear strength is determined for a similar concrete and on a level of 19.1 MPa, which is depicted in Figures 13 and 14 by a dotted line. In this case, the difference between the experimental and the numerical results is higher by about 25% but could be accepted.
Additionally, the extreme values of the macro-stresses and the corresponding macrostrains determined in the numerical tests are collected in Table 4.

Concluding Remarks
In this paper, the technique of numerical homogenization was used for determining the macroscopic response of high-performance concrete modelled as a nonhomogeneous multi-phase nonlinear composite. The local plastic-damage boundary value problem defined on a representative volume element of the composite with complex micro-structure was solved by making use of the finite element method. Three types of boundary condition were imposed on the 2-dimensional RVE with nonlinear components: linear displacement boundary conditions (DBC), uniform traction boundary condition (TBC), and periodic boundary condition (PBC). The impact of the boundary conditions on the determined macroscopic response was studied. The obtained results show that the applied boundary conditions exert a significant influence on the calculated values of macro-quantities and that premature numerical stability problems may occur in the case of TBC.
Detailed conclusions concerning the investigated influence of a type of boundary condition on the macro-results could be summarized as follows: Validation of the numerical results by our own laboratory experiments and the experimental data from the literature shows that the presented model is capable of providing good estimations of macroscopic features for the investigated deformation process. In the first step, we restricted ourselves to an analysis of the micro-scale BVP on the RVE. In the future, it would be useful to enrich the presented analysis by using the concept of modelling error and adaptivity [65,66]. This is especially useful in the circumstance, where the behavior of a structural member as a whole is the objective and then the transition of information between the macro-and micro-scales must be performed iteratively at all integration points of the macro-finite element mesh.