An electro-mechanically coupled computational multiscale formulation for electrical conductors

Motivated by experimental findings on deformation induced microcracks in thin metal films and by their influence on the effective macroscopic electrical conductivity, a computational multiscale formulation for electrical conductors is proposed in this contribution. In particular, averaging theorems for kinematic quantities and for their energetic duals are discussed, an extended version of the Hill–Mandel energy equivalence condition is proposed and suitable boundary conditions for the microscale problem are elaborated. The implementation of the proposed framework in a two-scale finite element environment is shown and representative boundary value problems are studied in two- and three-dimensional settings.


Introduction
Advances in material science and fabrication technologies enable the development of flexible electronic devices like wearable sensors [1,22] and foldable displays [3,15], which are in the focus of many engineering applications. To ensure the functionality of these devices, detailed knowledge is required on the electro-mechanical material properties for various loading conditions. However, the experimental characterization of material thin films, especially the study of failure mechanisms, is difficult and requires the development of advanced measuring technologies as discussed in [4]. Computational multiscale homogenization schemes are promising numerical approaches to support these developments as they allow material microstructures and their evolutions to be resolved in numerical simulations. On the other hand, complex deformation processes which are difficult to study experimentally are accessible using suitable numerical multiscale formulations that were calibrated on the basis of experiments.
More specifically speaking, processes at different material length-and timescales are considered in finite element-based computational multiscale simulations with the evaluation of classic constitutive material models being replaced by finite element calculations of the underlying material microstructures. The latter ones are represented by representative volume elements (RVEs) that take distinct microstructural features like grains and microcracks into account [9,27]. The boundary conditions on the microscale are prescribed based on the macroscopic material state. Vice versa, effective energetic duals to the kinematic quantities on the macroscale are calculated based on their microscopic counterparts by using computational homogenization schemes. First-order computational multiscale formulations for mechanical problems have in the meanwhile been wellestablished and are for example discussed in [10,11,16,[19][20][21]. Homogenization methods for generalized continua, in particular, second-order homogenization methods that account for the macroscopic deformation tensor and for its gradient, are elaborated in [12,17,18,23]. Moreover, extensions to thermomechanical coupling are presented in [2,6,24,28,29] and electro-mechanically coupled multiscale formulations for piezo-and ferroelectric materials are discussed in [13,14,26].
The present contribution establishes an electro-mechanically coupled computational multiscale framework for electrical conductors. In contrast to the theoretical developments on (dielectric) electro-active solids discussed in [13,26], the electro-mechanical coupling in the electrical conductors considered is rather implicit. In particular, changes in electrical conductivity due to mechanically induced microcracks motivate the present contribution. Against this background, assume for now a homogeneous, quasi-one-dimensional electrical problem for which the effective macroscopic electrical resistance is given by Whereas the geometrical contribution R geo to the electrical resistance is defined by the macroscopic length L and cross-sectional area A of the conductor, the specific resistance R spe which is defined as the inverse conductivity κ is an effective macroscopic material parameter that accounts for the (possibly inhomogeneous) material microstructure. The different contributions to the effective macroscopic resistance are exemplified in Fig. 1 for a material featuring microstructural imperfections. Furthermore, the influence of plastic deformation and of grain boundaries on the electrical conductivity may be taken into account. With regard to the influence of the electrical field quantities on the mechanical field quantities, it is observed that the dissipation associated with the electric current contributes to the balance equation of energy and that the induced changes in temperature, in turn, influence the mechanical properties. Finally, mechanically induced deformations directly influence the geometric part of the electrical resistance in finite deformation processes. The article is organized as follows: After a brief summary of the governing field equations for electromechanical problems in Sect. 2, Sect. 3 focuses on the electro-mechanically coupled multiscale formulation. In particular, averaging theorems for kinematic quantities and for their energetic duals are elaborated, generalized Hill-Mandel conditions are derived and suitable boundary conditions are discussed. Based on these developments, a finite element implementation of the proposed theory is presented in Sect. 4. A study of representative boundary value problems in Sect. 5 eventually shows the applicability of the proposed formulation.

Notation
Let α, β, γ and δ denote arbitrary first-order tensors, let the standard dyadic product be indicated by ⊗ and let single tensor contractions be given by [α ⊗ β] · γ ⊗ δ = β · γ [α ⊗ δ]. Based on the latter definitions, the generalized dyadic products are introduced to shorten notation. Moreover, these serve as the basis for the introduction of higher-order identity tensors based on suitable combinations of second-order identity tensors I. Double tensor contractions are defined as [α ⊗ β] : γ ⊗ δ = α · γ [β · δ] and gradient, divergence and curl operators, Fig. 1 Differences in the geometric dimensions and in the underlying microstructures of the idealized wire-like structures (left, right) manifest themselves in different geometric and specific resistances compared to the reference specimen (middle). Microstructure images courtesy of Megan J. Cordill, Erich Schmid Institute of Materials Science, Austrian Academy of Sciences indicated by ∇•, by ∇ · • and by ∇ × •, respectively, are used in sense of right-gradient, right-divergence and right-curl operations.

Continuum thermodynamics
This section provides a summary of the thermodynamic fundamentals of the (thermo-)electro-mechanically coupled multi-field formulation. To this end, let the region that is occupied by the body under consideration at some reference time t 0 ∈ R be denoted by B 0 ⊂ R 3 . Moreover, define the deformation field as the point mapping ϕ (X, t) : B 0 ×R → B t ⊂ R 3 that relates the referential position of a particle X ∈ B 0 to its counterpart x ∈ B t in the current configuration B t ⊂ R 3 at time t ≥ t 0 . As is customary for infinitesimal deformation processes which are assumed in the following, the displacement field u = ϕ (X, t) − X is additionally introduced and the reference and current configuration will not be distinguished, e.g., when defining integrals or derivatives, such that the respective superscripts are omitted in the following.

Mechanical field equations
Considering mass to be conserved, the mechanical response is governed by the balance equation of linear momentum with the mass density per unit volume ρ, the velocity vectoru and the body force density (per unit mass) f . By invoking Cauchy's theorem to relate the stress vector t to the stress tensor σ and to the outward unit surface normal n according to and by applying Gauss's theorem, (2) can be localized Using (3) and (4), the balance equation of angular momentum with x = x − x ref denoting the difference vector to a fixed but otherwise arbitrary reference point x ref , reduces to the symmetry condition of the stress tensor

Electrical field equations
General electro-magnetic problems are described in terms of Maxwell's equations. Being more specific, these are Gauss's magnetic law Faraday's law of induction Ampère's circuital law with the electric field vector e, the dielectric displacement vector d, the electric current density vector j , the magnetic field vector h, the magnetic flux density vector b and with the free charge density ρ f . Unlike ∂B, area A needs not to be, and is generally not, closed. However, evaluating (7d) for a closed surface and inserting (7a) gives rise to the continuity equation for the electric current which is of particular importance for the developments to be presented in this contribution. Localizing the set of equations (7) by using Gauss's and Stokes' theorem yields Gauss's electric law Ampère's circuital law ∇ × h = j +ḋ (9d) and the localization of (8) results in Focusing on (quasi-)stationary electrical processes of conducting materials, Faraday's law (9c) reduces to ∇ × e = 0 (11) and the continuity equation for the electric current density is given by while the evaluation of Gauss's electric law (9a) yields the free charge density at every point of the continuum.

Conservation of energy
The balance equation of energy which governs the temperature evolution is given by where, in addition to the mechanical working by body forces and surface tractions, thermal and electrical contributions are assumed to contribute to the (mass-specific) internal energy density e. The thermal contributions are given in terms of the heat source (per unit mass) r , which may, e.g., be related to chemical reactions, as well as in terms of the heat flux vector q. The electrical contribution is given in terms of the electric current density vector j and its energetic dual, the electric field vector e, see [5,8,25]. Localizing (13) by invoking Gauss's theorem and by making use of (3), (4) and (6) yields where the definition of the small strain deformation tensor was used.

Dissipation inequality
By introducing the (mass specific) entropy density s, the dissipation inequality characterizing the (thermo-)electro-mechanically coupled continuum reads The localization of (16) by using Gauss's theorem yields Moreover, introducing the (mass-specific) Helmholtz free energy density function ψ, making use of the convexconcave Legendre(-Fenchel) transformation with the dual variables θ and −s defined as and inserting the local form of the balance equation of energy (14) yields In addition to the well-established Clausius-Planck and Fourier inequalities, (20) poses restrictions on the specific form of the constitutive equation relating the electric field vector to the electric current density vector.
In particular, it is observed that the (volume-specific) dissipation contribution j · e is associated with the electric current. The physical processes of electrodynamics, analyzed in the present contribution, are thus fundamentally different from the (possibly) reversible processes that occur in electro-active solids and that are studied in, e.g., [26].

Multiscale modelling
In computational multiscale formulations effective quantities at the macroscale are related to their counterparts at the microscale. This procedure is intrinsically based on averaging theorems for stresses, strains, electric current densities, and electric field strengths as discussed in Sect. 3.1. In addition, requiring micro-macro energy equivalence gives rise to the extended Hill-Mandel conditions presented in Sect. 3.2. Different sets of boundary conditions which a priori fulfil the averaging theorems and the extended Hill-Mandel conditions are finally established in Sect. 3.3. Due to the structural similarity of the continuity equation for the electric current density and the balance equation of linear momentum, and since the electric field is derived from a potential analogous to the derivation of the strain tensor from the displacement field, the ensuing derivations for the electrodynamic problem are similar to those of mechanical problems as documented in, e.g., [10,16,[19][20][21]. Moreover, the structure of the underlying set of equations shares similarities with thermo-mechanically coupled problems as discussed in, e.g., [24] and with electro-mechanically coupled problems (for dielectric solids) as discussed in, e.g., [26]. Despite these structural similarities, there are significant differences between the present theory and the developments on electro-active solids presented in, e.g., [26]. The derivations presented in [26] are for example based on the assumptions of electrostatics while the present formulation relies on the fundamentals of electrodynamics such that different balance equations are solved, different field quantities occur, and different constitutive equations are postulated. In particular, the electric current is intrinsically dissipative as opposed to the (possibly) reversible processes in dielectric solids analyzed in [26].
In order to distinguish micro-and macroscale quantities, subscripts m(icro) and M(acro) are introduced. Moreover, quasi-static and quasi-stationary processes are assumed at the microscale. By additionally assuming negligible body forces, the set of balance equations to be solved at the microscale takes the form with Faraday's law (11) being accounted for by the introduction of an electric potential φ m such that The assumptions that give rise to (21) are well-established for the mechanical problem and are for example discussed in [28].

Averaging theorems
In computational multiscale formulations, effective macroscopic quantities are defined as volume averages of their microscopic counterparts. For the evaluation of the generalized Hill-Mandel conditions, cf. Sect. 3.2, and from an implementation point of view it is moreover beneficial to rewrite the occurring volume integrals as surface integrals using Gauss's theorem. Invoking the identity and introducing the RVE-volume the effective macroscopic strain tensor may be specified in terms of boundary displacements, namely, Analogously, the effective macroscopic electric field vector is given by Moreover, the governing equation for the effective macroscopic stress tensor reads and the one for the electric current density vector results in

Hill-Mandel conditions
The micro-macro energy equivalence condition, known in the literature as the Hill-Mandel condition, states that the volume average of the microscopic power in the RVE equals the local power on the macroscale. Focusing on the electro-mechanical problem introduced in Sect. 2 and denoting variations of field quantities • by δ•, the classic Hill-Mandel condition for the mechanical problem reads and the one for the electrical problem is given by Comprehensive discussions of the Hill-Mandel conditions for thermo-mechanically coupled problems, which are not in the focus of the present contribution, are presented in, e.g., [6,28,29].

Boundary conditions
Different types of boundary conditions for the microstructural representative volume element can be defined which ensure that the Hill-Mandel conditions and the averaging theorems are fulfilled a priori. In the following, the classic three types of boundary conditions, i.e., affine and periodic as well as uniform flux boundary conditions, are briefly recapitulated for mechanical problems and the extension to electrical problems is shown.

Affine boundary conditions
Based on the macroscopic deformation state, affine displacement boundary conditions for the representative volume element are given by Analogously, taking into account the macroscopic electric potential and the macroscopic electric field strength, affine boundary conditions for the electric potential are defined as The particular form (32) has been chosen in accordance with thermo-mechanical multiscale formulations as discussed in [2,6]. However, in contrast to the temperature, the electric potential is assumed to enter the constitutive equations only in terms of its gradient, i.e., the electric field vector, such that alternative definitions that do not account for φ M are possible. The boundary conditions (31) and (32) a priori fulfil the averaging theorem for the strain tensor (25) and for the electric field vector (26) as well as the mechanical and the electrical Hill-Mandel condition, (29) respectively (30). More specifically speaking, by evaluating (25) for the boundary displacements given by (31) one arrives at and by inserting (32) into (26) at By additionally inserting (32) into the Hill-Mandel condition for the mechanical problem (29) and by making use of (27) one finds Likewise, the evaluation of (30) for (32), using (28), yields

Periodic boundary conditions
Indicating quantities at two opposing parts of the boundary, ∂B + m and ∂B − m , with superscripts • + and • − , the periodicity conditions for the mechanical problem are given by the ones for the electrical problem read and the classic geometric constraint for the unit outward surface normal vectors is stated as The set of equations (37) and (38) implies periodicity of the displacements and of the electric potential as well as anti-periodicity of the tractions and of the electric current densities. Inserting (37a) into the averaging theorem for the strain tensor (25), using (39) and applying Gauss's theorem in the same way as in (33), yields Analogously, by inserting (38a) into (26) one finds Evaluating the mechanical and electrical Hill-Mandel conditions, (29) and (30), for the periodic boundary conditions (37) and (38) eventually yields

Uniform flux boundary conditions
As opposed to affine and periodic boundary conditions, Neumann-type boundary conditions in terms of tractions and electric current densities are prescribed at the boundary of the representative volume element when uniform flux boundary conditions are applied. For the mechanical boundary value problem these are given by and those for the electrical problem read The equality of the macroscopic stress tensor σ M and of the averaged microscopic stress tensor for uniform traction boundary conditions follows immediately by inserting (44) into (27), to be specific Analogously, the insertion of (45) into (28) yields The evaluation of the Hill-Mandel conditions (29) and (30) furthermore results in and in where (25) and (26) were used.

Finite element implementation
This section focuses on the finite element implementation of the electro-mechanically coupled multiscale formulation proposed in Sect. 3. In particular, the homogenization of effective macroscopic stresses and electric current densities is discussed in Sect. 4.1, and the corresponding algorithmic tangent stiffness tensors are derived in Sect. 4.2. The ensuing derivations rely on and are direct extensions of well-established computational homogenization procedures. For a detailed elaboration of these fundamentals, the reader is referred to, e.g., the pioneering work [19] which focuses on the purely mechanical case. The extension to thermo-mechanically coupled multi-field problems is moreover discussed in [24,28] and electro-mechanically coupled problems (of dielectric materials) are for instance studied in [13,14]. The subsequent derivations apply to affine and periodic boundary conditions. In the case of periodic boundary conditions, however, the generalized stiffness matrix K m and the generalized reaction force vector f m must be substituted by the reduced generalized stiffness matrix K * m and by the reduced generalized reaction force vector f * m that result from the application of linear constraints when enforcing (anti-)periodicity. Specifically speaking, the system is rewritten in each iteration step in the form using the transformation matrix T that relates dependent X md and independent X mi degrees of freedom, see [30].

Homogenization
Following standard procedure, the discrete version of (27) is given by with n pn denoting the number of nodes where the corresponding degrees of freedom are prescribed, and with (i) x m and (i) f u m denoting the position vector and the reaction force vector of node i. Analogously, the discrete version of (28) is given by with (i) f φ m denoting the generalized reaction force of node i, corresponding to the electrical problem.

Generalized algorithmic tangent stiffness tensors
In order to derive generalized algorithmic tangent stiffness operators, changes in the generalized reaction forces f mp due to perturbations of the prescribed degrees of freedom X mp are studied. To this end, the partitioned system of equations with subscripts • p and • f referring to prescribed and free degrees of freedom is considered. Solving the second set of equations for X mf and inserting the result into the first set of equations results in For affine as well as for periodic boundary conditions, the kinematic relations and hold for nodes where the respective degrees of freedom are prescribed. In the case of affine boundary conditions, these are the boundary nodes while the set consists of three corner nodes in a two-dimensional setting, respectively, of four corner nodes in a three-dimensional setting, when periodic boundary conditions are applied. By partitioning (56) into mechanical and electrical contributions and after inserting the kinematic relations (57) and (58), it is observed that changes in the generalized reaction force vectors at the microscale due to changes in ε M and e M are given by and In accordance with (53) and (54), superscripts • u and • φ refer to mechanical and electrical field equations and fields, respectively. Moreover, it is assumed that the electric potential enters the field equations and constitutive equations only via its gradient, i.e., via the electric field vector, such that changes in φ M are neglected in (59), (60) and in the following. A dependence of the constitutive equations on the electric potential would lead to problems analogous to those of the thermal case, as discussed in [29]. Eventually, the insertion of (59) into (53) gives rise to the definition of the macroscale algorithmic tangent stiffness tensors Likewise, inserting (60) into (54) yields

Representative simulation results
This section focuses on the application of the proposed electro-mechanical multiscale formulation to two-and three dimensional quasi-static and quasi-stationary boundary value problems. To this end, the material model that is used at the microscale is briefly discussed in Sect. 5.1. In a second step, the calculation of effective electrical conductivity tensors for two-dimensional microstructures is discussed in Sect. 5.2, before tension test samples with different microstructures are subjected to detailed three-dimensional multiscale analyses in Sect. 5.3.

Microscale material models
On the microscale, constitutive models that characterise the material behavior need to be developed subject to the restrictions posed by the dissipation inequality (20). Moreover, the specific form of the constitutive equations is decisive for the coupling between the electrical and the mechanical field equations. Different physical effects that result in an explicit and implicit coupling of the field equations were outlined in Sect. 1 and motivate the present contribution. However, since the focus of the present work lies on the fundamentals of the electro-mechanically coupled computational multiscale formulation and not on the development of a particular microscale material model, we restrict ourselves to isotropic, linear electrical and mechanical microscale material models and postpone the analysis of more elaborated material models to future works. Thus, for the particular material model considered, the electrical and the mechanical field equations remain uncoupled. Moreover, superscripts • m referring to the microscale are neglected in the following for the sake of brevity. By introducing the volume-specific Helmholtz free energy density function Ψ (ε, θ, •) = ρ ψ (ε, θ, •), the evaluation of the dissipation inequality (20) yields the standard definition of the stress tensor In conjunction with the specific quadratic form of the volume specific Helmholtz free energy density function with E denoting the fourth-order elasticity tensor and with E and ν denoting the Young's modulus and Poisson's ratio, the evaluation of (63) results in In addition, by introducing the positive definite electrical conductivity tensor with the idealized scalar-valued conductivity κ, a suitable choice for the electric current density vector that is in accordance with (20) is given by

Two-dimensional representative simulations
This section focuses on the application of the proposed electro-mechanical multiscale formulation in a twodimensional plane strain setting. To this end, idealized material parameters of copper are assumed at the microscale, and effective macroscopic conductivity tensors are exemplarily calculated for different microstructures. The material parameters are summarized in Table 1 and three different microstructures to be analyzed are depicted in Fig. 2.
The microstructures are assumed to be periodic and the representative volume elements which form the basis for the finite element simulations are indicated by dashed lines. For the sake of clarity, only the dimensions in e 1 -direction are provided. However, since symmetry with respect to the e 1 -and e 2 -direction or to the space diagonals is assumed, the remaining dimensions follow implicitly.
The application of the homogenization scheme with periodic boundary conditions as discussed in Sect. 4 to the microstructures depicted in Fig. 2  with superscripts referring to the respective microstructure. Comparing the effective macroscopic conductivity tensor of a material with circular voids (69a) with the conductivity tensor of an idealized material (67), it is observed that the effective macroscopic constitutive response remains isotropic whereas a reduction in the conductivity of about 45% is observed. In contrast, the effective macroscopic conductivity tensor (69b) that results from a microstructure with a groove-shaped void as depicted in Fig. 2b is significantly anisotropic. However, the principal material axes are aligned with the e 1 -and e 2 -coordinate axes such that the conductivity tensor is in diagonal form. In the case of the microstructure of diagonal type depicted in Fig. 2c, the e 1 -and e 2coordinate axes are not aligned with the principal material axes. Thus, the effective macroscopic conductivity tensor (69c) is not in diagonal form but features off-diagonal components.
The previous examples demonstrate the applicability of the electro-mechanical multiscale formulation to extract effective macroscopic conductivity tensors for given microstructures. This approach may be interesting Tension test specimen of type E according to DIN 50125, see [7], used in the electro-mechanical multiscale finite element simulations. Dark grey color indicates regions where Dirichlet boundary conditions are applied. The light grey colored 60 mm × 10 mm × 6 mm region in the center of the specimen resembles material featuring microscale imperfections. All dimensions are given in mm for electrical engineering applications to developed tailored microstructures with specific electrical properties. Moreover, regarding electro-mechanical coupling, changes in effective macroscopic conductivities that are accessible via non-destructive testing methods may be related to microcracks that are induced by mechanical loadings.

Three-dimensional representative simulations
As an extension of the two-dimensional sample boundary value problems discussed in Sect. 5.2 to a threedimensional setting, this section focuses on the effective macroscopic responses of tensile test specimens with different material microstructures. The tensile test specimen of type E that serves as the basis for the simulations is sketched in Fig. 3, with supports at the left and right boundaries being indicated by dark grey color. Specifically speaking, homogeneous Dirichlet boundary conditions for the electrical and mechanical fields are assumed at the top and bottom surface nodes of the left boundary. Likewise, homogeneous Dirichlet boundary conditions for the mechanical field are applied in e 2 -and e 3 -direction at the top and bottom nodes of the right boundary, while both the displacement in e 1 -direction as well as the electric potential take values different from zero, see Table 2.
In order to study the influence of different microstructures on the effective electrical and mechanical properties at the macroscale, simulation results of a tensile test specimen with a perfect microstructure as depicted in Fig. 4b are compared to those of a tensile test specimen with an imperfect microstructure. Specifically speaking, imperfections on the microscale are assumed in the light grey colored 60 mm × 10 mm × 6 mm region in the center of the specimen and are given in terms of a spherical void with a radius of 0.23 RVE-edge length as depicted in Fig. 4d. In accordance with the two-dimensional simulation results presented in Sect. 5.2, isotropic, linear constitutive relations are assumed for the mechanical and electric problem at each material point at the microscale. Furthermore, the same set of material parameters as in the two-dimensional simulations is used, see Table 1, and periodic boundary conditions are applied at the microscale.
The distributions of the electric current density at the micro-and macroscale are provided in Fig. 4 for tensile test samples with perfect and imperfect microstructures. By comparing the distributions of the electric current density at the macroscale, Fig. 4a, c, it is observed that the material imperfection manifests itself in a reduction of the electric current density in the center region of the specimen. This effect can be explained by  , d), near the center of the specimen, are depicted. The simulation results for a tensile test specimen with a perfect microstructure (a, b) differ significantly from the ones with an imperfect microstructure that contains voids (c, d).
Gray-colored arrows indicate the electric current density vector. Only one half of the representative volume element is shown the microscale simulation results depicted in Fig. 4b, d which show that the material imperfection acts as an obstacle for the electric current. In addition, the force in e 1 -direction and the electric current that both occur as generalized reaction forces dual to the applied displacement and to the applied electric potential difference, respectively, are summarized in Table 2. In accordance with the simulation results of Fig. 4, the effective macroscopic electric current takes smaller values in the case of a material with an imperfect microstructure. The same observation is made for the reaction force in e 2 -direction.

Closure
Summary Motivated by advances in flexible electronic technologies and by the desire to develop non-destructive testing methods, an electro-mechanical multiscale formulation for conductors is proposed in this contribution. After a brief recapitulation of the governing set of electro-mechanical field equations, scale-bridging relations for the electro-mechanical problem are proposed. In addition, a finite element implementation of the electromechanical multiscale formulation is discussed, and different types of boundary condition are shown to a priori fulfil the extended Hill-Mandel energy equivalence conditions. Several numerical examples are studied in two-and three-dimensional settings to show the applicability of the proposed formulation. In particular, the calculation of effective macroscopic conductivity tensors for given two-dimensional microstructures is discussed, which is of interest for electrical engineering applications when developing tailored microstructures with specific electrical properties. Moreover, the study of tensile test samples with different microstructures, which are subjected to electro-mechanical boundary conditions, reveals the influence of the microstructure on the effective macroscopic electrical and mechanical material properties. The interrelation of the latter is of special interest in the development of nondestructive testing methods, and is in the focus of future works.
Outlook Multiscale approaches are based on an accurate simulation of processes that occur at a lower scale. For this reason, more elaborated electro-mechanically coupled material models are to be developed and to be used at the microscale, e.g., in order to account for the influence of mechanically induced microcracks on the electrical conductivity. Against this background and with tailored microscale material models being beyond the scope of the present contribution, assume for now a deformation-dependent conductivity tensor of the form and with material parameter κ * controlling the decrease in electrical conductivity with increasing deformation. The particular form of (70) is based on the idea that damage evolution is intrinsically energy-driven and that cracks represent obstacles for the electric current. Of course, the simplified format considered here does not include any history variables.
To study the latter material model in a multiscale environment, a two-dimensional plane strain one-element test as depicted in Fig. 5 is analyzed, with a periodic microstructure according to Fig. 2a being assumed. For the mechanical sub-problem the material model (63)-(66) is adopted. Moreover, the material parameters are chosen according to Table 1 and κ * = 0.01 E is assumed.
The decrease in the 11-and 22-coefficient of the effective macroscopic conductivity tensor as a function of deformation in terms of the prescribed macroscale displacementū is provided in Fig. 6a, with a significantly different decrease of the electrical conductivity in the two spatial directions being observed. The latter finding can be explained by the inhomogeneous microscale deformation as exemplarily shown in Fig. 6b. For the analyzed load state, the severe deformation zones that occur at the top and bottom of the circular void result in a significant decrease of the electrical conductivity in e 1 -direction.
The derivation of S is motivated by geometrically induced changes in resistance of an idealized wire according to (1). Additionally, linearized versions of the transport theorems for infinitesimal line elements