Optical properties of two-dimensional polymer photonic crystals after deformation-induced pattern transformations

The photonic band structure and optical transmittance of two-dimensional periodic elastomeric photonic crystals are studied computationally to understand the effects of large strains on optical properties of the structures. The large compressive deformation patterns of the two-dimensional periodic structure studied by Mullin and coworkers [Mullin, T., Deschanel, S., Bertoldi, K., Boyce, M.C., 2007. Pattern transformation triggered by deformation. Physical Review Letters 99(8), 084301] are first reproduced using hyperelastic material models for the elastomer SU-8. Finite element analysis is then used to solve Maxwell's equations to obtain light transmittance through both the undeformed and deformed structures; simultaneously the wave equation resulting from the appropriate two-dimensional form of Maxwell's equations is solved as an eigenvalue problem to obtain the band structure. The deformation-induced shift in transmission spectrum valleys for different bands is calculated, and the changes in the width of these reflectance peaks are also obtained. The band structure calculation shows that there are no complete photonic band gaps as expected for the low dielectric contrast system. However, the effect of the observed reversible, symmetry-breaking deformation pattern is to uncouple many of the photonic bands in all three high symmetry directions, i.e. Γ–X, X–M, and Γ–M. New non-degenerate deformation-induced optical modes appear in both the real space transmittance spectra and the band structure with lower reflectance values. Analyses of the deformation pattern, the optical mode shapes, and the photonic band structure reveal that localized regions of large rotation are responsible for the significant changes in optical transmittance. The results have practical importance for the design of strain-tunable optomechanical materials for sensing and actuation.     

An endochronic plasticity formulation for filled rubber

The nature of elastomeric material demands the consideration of finite deformations, nonlinear elasticity including damage as well as rate-dependent and rate-independent dissipative properties. While many models accounting for these effects have been refined over time to do better justice to the real behavior of rubber-like materials, the realistic simulation of the elastoplastic characteristics for filled rubber remains challenging.The classical elastic-ideal-plastic formulation exhibits a distinct yield-surface, whereas the elastoplastic material behavior of filled rubber components shows a yield-surface free plasticity. In order to describe this elastoplastic deformation of a material point adequately, a physically based endochronic plasticity model was developed and implemented into a Finite Element code. The formulation of the ground state elastic characteristics is based on Arruda and Boyce (1993) eight-chain model. The evolution of the constitutive equations for the nonlinear endochronic elastoplastic response are derived in analogy to the Bergström–Boyce finite viscoelasticity model discussed by Dal and Kaliske (2009).     

Modelling large deformation behaviour under loading–unloading of semicrystalline polymers: Application to a high density polyethylene

In this work, the large deformation behaviour under monotonic loading and unloading of a high density polyethylene (HDPE) is studied. To analyze the nonlinear time-dependent response of the material, mechanical tests were conducted at room temperature under constant true strain rates and stress relaxation conditions. A physically-based inelastic model written under finite strain formulation is proposed to describe the mechanical behaviour of HDPE. In the model, the inelastic mechanisms involve two parallel elements: a visco-hyperelastic network resistance acting in parallel with a viscoelastic–viscoplastic intermolecular resistance where the amorphous and crystalline phases are explicitly taken into consideration. The semicrystalline polymer is considered as a two-phase composite. The influence of the crystallinity on the loading and unloading behaviour is investigated. Numerical results are compared to experimental data. It is shown that the model is able to accurately reproduce the experimental observations corresponding to monotonic loading, unloading and stress relaxation behaviours at different strain levels. Finally, the model capabilities to capture cyclic loading–unloading behaviour up to large strains are discussed. To demonstrate the improved modelling capabilities, simulations are also performed using the original model of Boyce et al. [Boyce, M.C., Socrate, S., Llana, P.G., 2000. Constitutive model for the finite deformation stress–strain behavior of poly(ethylene terephthalate) above the glass transition. Polymer 41, 2183–2201] modified by Ahzi et al. [Ahzi, S., Makradi, A., Gregory, R.V., Edie, D.D., 2003. Modeling of deformation behavior and strain-induced crystallization in poly(ethylene terephthalate) above the glass transition temperature. Mechanics of Materials 35, 1139–1148].     

Implicit finite element formulations for multi-phase transformation in high carbon steel

An anomalous plastic deformation observed during the phase transformation of steels was implemented into the finite element modeling. The constitutive equations for the transformation plasticity originally proposed by Greenwood and Johnson [Greenwood, G.W., Johnson, R.H., 1965. The deformation of metals under small stresses during phase transformation. Proc. Roy. Soc. A 283, 403] and further extended by Leblond et al. [Leblond, J.B., Mottet, G., Devaux, J.C., 1986a. A theoretical and numerical approach to the plastic behavior of steels during phase transformations, I. Derivation of general relations. J. Mech. Phys. Solids 34, 395–409; Leblond, J.B., Mottet, G., Devaux, J.C., 1986b. A theoretical and numerical approach to the plastic behavior of steels during phase transformations, II. Study of classical plasticity for ideal-plastic phases. J. Mech. Phys. Solids 34, 411–432; Leblond, J.B., Devaux, J., Devaux, J.C., 1989a. Mathematical modeling of transformation plasticity in steels, I: case of ideal-plastic phases. Int. J. Plasticity 5, 511–572; Leblond, J.B., 1989b. Mathematical modeling of transformation plasticity in steels, II: coupling with strain hardening phenomena. Int. J. Plasticity 5, 573–591] were modified to consider the thermo-mechanical response of generalized multi-phase steel during phase transformations from austenite at high temperature. An implicit numerical solution procedure to calculate the plastic deformation of each constituent phase was newly proposed and implemented into the general purpose implicit finite element program via user material subroutine. The new algorithms include efficient calculation of consistent tangent modulus for the transformation plasticity and application of general anisotropic yield functions without limitation to the isotropic yield function. Besides the thermo-elastic–plastic constitutive equations, non-isothermal transformation kinetics was characterized by the Johnson–Mehl–Avrami–Kolmogorov (JMAK) equation and additivity relationship for the diffusional transformation, while the model proposed by Koistinen and Marburger was used for the diffusionless transformation. Numerical verifications for the continuous cooling experiments under various loading conditions were conducted to demonstrate the applicability of the developed numerical algorithms to the high carbon steel SK5.     

Heuristic search for a predictive strain-energy function in nonlinear elasticity

In this work, a new, quasi-structural model – bootstrapped eight-chain model – is proposed as a modification to the strain energy of eight-chain model [Arruda, E.M., Boyce, M.C., 1993. A three-dimensional constitutive model for the large stretch behaviour of rubber elastic materials. J. Mech. Phys. Solids 41, 389—412] that invokes the Langevin chain statistics. This development has been led to by our heuristic search into how the strain energy of eight-chain model may be adapted in order to account better for the mechanical behaviour of elastomeric materials in both linear and nonlinear elastic regimes [Treloar, L.R.G., 1944. Stress–strain data for vulcanised rubber under various types of deformation. Trans. Faraday Soc. 40, 59–70]. The eight-chain model appears to produce very similar results in predicting biaxial stress to those of a first stretch-invariant model that gives a good fit in uniaxial extension and, thus, it is shown that the former can not be significantly enhanced within the limitation of the latter. Evaluation of predictive capability for an additive invariant-separated form of strain energy shows that an explicit inclusion of a second stretch-invariant function would not work and that any thus added term ought to be dependent on both the first and second stretch-invariants of deformation tensor, and hints that an improvement is possibly needed at low strain. The composite and filament models [Miroshnychenko, D., Green, W.A., Turner, D.M., 2005. Composite and filament models for the mechanical behaviour of elastomeric materials. J. Mech. Phys. Solids 53 (4), 748–770] have their strain-energy functions in that suggested form and cope very well with predicting the experimental data of Treloar (1944). We use the form of strain energy for the filament model, that proved to be successful, to bootstrap the strain energy of eight-chain model in order to improve the performance of the latter at low strain. Thus, we derive a new model – bootstrapped eight-chain model – that requires only two material parameters – a rubber modulus and a limiting chain extensibility. The proposed model is quasi-structural due to bootstrapping and it retains the best traits and corrects the faults of the eight-chain model, conforming more closely to the classical experimental data of Treloar (1944).     

Numerical modelling of the plasticity induced during diffusive transformation. An ensemble averaging approach for the case of random arrays of nuclei

The grounds of a numerical modelling of the mechanical consequences of diffusive phase transformation in solids have been established by Leblond [Leblond, J.B., Mottet, G., Devaux, J.C., 1986. A theoretical and numerical approach to the plastic behavior of steels during phase transformations I: derivation of general relations, J. Mech. Phys. Solids 34 (4) 395–409] and Ganghoffer [Ganghoffer, J.F., Denis, S., Gautier, E., Simon, A., Sjöström, S., 1993. Finite element calculation of the micromechanics of a diffusional transformation, Eur. J. Mech. A Solids 12 (1) 21–32]: this modelling resorts to the FE method to evaluate the stress and strain fields which ensure the mechanical equilibrium between a diffusionaly growing phase and its parent phase. It has been the subject of a thorough analysis in [Barbe, F., Quey, R., Taleb, L., 2007. Numerical modelling of the plasticity induced during diffusive transformation. Case of a cubic array of nuclei, Eur. J. Mech. A Solids 26, 611–625] which has evidenced the main limit underlying this modelling with regard to physics, relative to the fact that nuclei are implicitly positioned according to a periodic array. The present work proposes, in details, an extension to the case of nuclei instantaneously appearing at random positions in a quasi infinite homogeneous medium. It constitutes a fundamental step towards a numerical modelling explicitly taking into account the crystalline plasticity and the morphology of the transforming medium.     

Mechanically tunable phononic band gaps in three-dimensional periodic elastomeric structures

Three-dimensional periodic structures have many applications in acoustics and their properties are strongly related to structural details. Here we demonstrate through simulations the ability to tune the phononic band gaps of 3D periodic elastomeric structures using deformation. The elastomeric nature of the material makes the transformation of the band gaps a reversible and repeatable process, providing avenues for the design of tunable 3D phononic crystals such as sonic switches.     

Electro-elastomers: Large deformation analysis of silicone membranes

Electro-elastomers are large strain smart materials capable of both sensing and actuation. Typical electro-elastomer setups consist of either a silicone or acrylic membrane sandwiched between two compliant grease electrodes. Silicone electro-elastomers have maximum elastic strains between 200% and 350%. Acrylic electro-elastomers are more widely employed due to larger actuation strains but are softer than silicone and have a lower force output [Goulbourne, N.C., Frecker, M., Mockensturm, E.M., Snyder, A.J., 2003. Modeling of a dielectric elastomer diaphragm for a prosthetic blood pump. In: Proceedings of SPIE, Smart Structures and Materials: EAPAD, San Diego; Goulbourne, N.C., Mockensturm, E.M., Frecker, M., 2005b. Quasi-static and dynamic inflation of a dielectric elastomer membrane. In: Proceedings of SPIE, Smart Structures and Materials: EAPAD, San Diego]. A numerical formulation for the large deformation response of electro-elastomer membranes subject to electromechanical loading is derived in this paper. The approach is based on modifying the elastic membrane theory of Green, Adkins, and Rivlin [Adkins, J.E., Rivlin, R.S., 1952. Large elastic deformations of isotropic materials IX. The deformation of thin shells. Philosophical Transactions of the Royal Society of London. Series A Mathematical and Physical Sciences 244, 505–531; Green, A.E., Adkins, J.E., 1970. Large Elastic Deformations. Oxford University Press, London]. The electro-elastic stress state is defined as the combination of the electrical Maxwell stress and the mechanical stress for hyperelastic materials [Goulbourne, N.C., Mockensturm, E.M., Frecker, M., 2005a. A nonlinear model for dielectric elastomer membranes. ASME Journal of Applied Mechanics 72, (6) 899–906]. This paper augments our previous work by presenting a mathematical solution procedure for simulating the field responsive behavior of silicone electro-elastomers configured for both in-plane and out-of-plane deformation. Thin axisymmetric membranes subject to electromechanical loads are the focus of this investigation. The numerical analysis shows that there is a delicate balance between the electrical and the mechanical portions of the stress, which must be maintained for the overall stress to remain tensile and by extension the electro-elastomer to remain stable. It is shown that at very high voltages the stress can become negative ultimately leading to transducer failure. For sensing applications, the varying capacitive behavior of electro-elastomers is used to extract information about the membrane’s deformed state.     

An inclusion theory for the propagation of martensite band in NiTi shape memory alloy wires under tension

In this paper, we present analytical modelling for the pure mechanical response of uniaxial tensioned NiTi wire that experiences stress-induced martensitic transformation via a propagating martensite band at the superelastic temperature regime. The model aims to predict the overall behavior of the SMA wire as a structural response containing propagating instabilities. Based on the systematic experimental investigation of Shaw and Kyriakides (Shaw, J.A., Kyriakides, S., 1995. Thermomechemical aspects of NiTi. J. Mech. Phys. Solids 43, T243–1281 and Shaw, J.A., Kyriakides, S., 1997. On the nucleation and propogation of phase transformation fronts in a NiTi alloy. Acta Mater. 45(2), 638–700), the wire is modeled as an elastic rod containing a single cylindrical transformation inclusion with a uniform axisymmetric eigenstrain. The analytical expression of the free energy of this special matrix-inclusion system is formulated and the length of the martensite band is identified as the key variable describing the transformation process of the system. Theoretical predictions on the peak stress and the subsequent steady-state propagation stress of the wire during forward and reverse transformations are provided and compared with the available experimental data. Specimen size effect on the nominal stress-strain curves and general deformation features of the wire are discussed.     

On modeling the micro-indentation response of an amorphous polymer

Interest in instrumented indentation experiments as a means to estimate mechanical properties has grown rapidly in recent years. Although numerous nano/micro-indentation experimental studies on polymeric materials have been reported in the literature, a corresponding methodology for extracting material property information from the experimental data does not exist. This situation for polymeric materials exists primarily because baseline numerical analyses of sharp indentation using appropriate large deformation constitutive models for the nonlinear viscoelastic–plastic response of these materials appear not to have been previously reported in the literature. An existing, widely used theory for amorphous polymers (e.g. [Boyce, M., Parks, D., Argon, A.S., 1988. Large inelastic deformation of glassy polymers. Part 1: Rate dependent constitutive model. Mechanics of Materials 7, 15–33; Arruda, E.M., Boyce, M.C., 1993. Evolution of plastic anisotropy in amorphous polymers during finite straining. International Journal of Plasticity 9, 697–720]) has been recently found to lack sufficient richness to enable one to quantitatively reproduce the major features of the indentation load-versus-depth curves for some common amorphous polymers [Gearing, B.P., 2002. Constitutive equations and failure criteria for amorphous polymeric solids. Ph.D. thesis, Massachusetts Institute of Technology].This study develops a new continuum model for the viscoelastic–plastic deformation of amorphous polymeric solids. We have applied the constitutive model to capture salient features of the mechanical response of the amorphous polymeric solid poly(methyl methacrylate) (PMMA) at ambient temperature and stress states under which this material does not exhibit crazing. We have conducted compression-tension strain-controlled experiments, as well as stress-controlled compression-creep experiments, and these experiments are used to calibrate the material parameters in the constitutive model for PMMA.We have implemented our constitutive model in a finite-element computer program, and using this finite-element program we have simulated micro-indentation experiments on PMMA. We show that our constitutive model and finite element simulations reproduce the experimentally-measured indentation load-versus-depth response with reasonable accuracy.     

Analysis of Lattices with Non-linear Interphases

Anti-plane deformation of square lattices containing interphases is analyzed. It is assumed that lattices are linear elastic but not necessarily isotropic, whereas interphases exhibit non-linear elastic behavior. It is demonstrated that such problems can be treated effectively using Green’s functions, which allow to eliminate the degrees of freedom outside of the interphase. Illustrative numerical examples focus on the determination of applied stresses leading to lattice instability.     

Micromechanics, macromechanics and constitutive modeling of the elasto-viscoplastic deformation of rubber-toughened glassy polymers

Under certain conditions, such as sufficiently low temperatures, high loading rates and/or highly triaxial stress states, glassy polymers display an unfavorable characteristic—brittleness. A technique used for reducing the brittleness (increasing the fracture toughness) of these materials is rubber toughening. While there is significant qualitative understanding of the mechanical behavior of rubber-toughened polymers, quantitative modeling tools for the large-strain deformation of rubber-toughened glassy polymers are largely lacking.In this paper, we develop a suite of numerical tools to investigate the mechanical behavior of rubber-toughened glassy polymers, with emphasis on rubber-toughened polycarbonate. The rubber particles are modeled as voids in view of their deformation-induced cavitation early during deformation. A three-dimensional micromechanical model of the heterogeneous microstructure is developed to study the effects of initial rubber particle (void) volume fraction on the underlying elasto-viscoplastic deformation mechanisms in the material, and how these mechanisms influence the macroscopic response of the material. A continuum-level constitutive model is developed for the large-strain elasto-viscoplastic deformation of porous glassy polymers, and it is calibrated against micromechanical modeling results for porous polycarbonate. The constitutive model can be used to study various boundary value problems involving rubber-toughened (porous) glassy polymers. As an example, the case of an axisymmetric notched bar is simulated for the case of polycarbonate with varying levels of initial porosity. The quality of the constitutive model calibration is assessed using a multi-scale modeling approach.     

A Variational Model for Reconstructive Phase Transformations in Crystals,and their Relation to Dislocations and Plasticity

1We study the reconstructive martensitic transformations in crystalline solids (i.e., martensitic transformations in which the parent and product lattices have arithmetic symmetry groups admitting no finite supergroup), the best known example of which is the bcc–fcc transformation in iron. We first describe the maximal Ericksen-Pitteri neighborhoods in the space of lattice metrics, thereby obtaining a quantitative characterization of the weak transformations, which occur within these domains. Then, focusing for simplicity on a two-dimensional setting, we construct a class of strain-energy functions admitting large strains in their domain, and which are invariant under the full symmetry group of the lattice. In particular, we exhibit an explicit energy suitable for the square-to-hexagonal reconstructive transformation in planar lattices. We present a numerical scheme based on atomic-scale finite elements and, by means of our constitutive function, we use it to analyze the effects of transformation cycling on a planar crystal. This example illustrates the main phenomena related to the reconstructive martensitic phase changes in crystals: in particular, the formation of dislocations, vacancies and interstitials in the lattice.     

A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

多孔极化PZT95/5铁电陶瓷单轴压缩力学响应与放电特性

采用添加造孔剂的方法制备了四种不同孔隙率PZT95/5铁电陶瓷,对其进行电场极化,随后开展了准静态单轴压缩实验,讨论了畴变、相变以及孔隙率对极化PZT95/5铁电陶瓷的力学响应与放电特性的影响. 研究结果表明：(1)多孔极化PZT95/5铁电陶瓷非线性力学响应行为主要归因于畴变和相变的共同作用,与内部孔洞变形和坍塌基本无关;(2) 在准静态单轴压缩下极化PZT95/5铁电陶瓷的去极化机制是畴变和相变的共同作用;(3) 孔隙率对极化PZT95/5铁电陶瓷的弹性模量、压缩强度有明显的影响,而对断裂应变的影响较小;(4)极化PZT95/5铁电陶瓷畴变和相变开始的临界应力都随着孔隙率的增大而线性衰减,但相变开始的临界体积应变却不依赖孔隙率;(5)极化PZT95/5铁电陶瓷电荷饱和释放量随着孔隙率呈线性减小,但孔隙率对电荷释放速率基本没有影响。     

Numerical analysis and elastic–plastic deformation behavior of anisotropically damaged solids

The present paper is concerned with the numerical modelling of the large elastic–plastic deformation behavior and localization prediction of ductile metals which are sensitive to hydrostatic stress and anisotropically damaged. The model is based on a generalized macroscopic theory within the framework of nonlinear continuum damage mechanics. The formulation relies on a multiplicative decomposition of the metric transformation tensor into elastic and damaged-plastic parts. Furthermore, undamaged configurations are introduced which are related to the damaged configurations via associated metric transformations which allow for the interpretation as damage tensors. Strain rates are shown to be additively decomposed into elastic, plastic and damage strain rate tensors. Moreover, based on the standard dissipative material approach the constitutive framework is completed by different stress tensors, a yield criterion and a separate damage condition as well as corresponding potential functions. The evolution laws for plastic and damage strain rates are discussed in some detail. Estimates of the stress and strain histories are obtained via an explicit integration procedure which employs an inelastic (damage-plastic) predictor followed by an elastic corrector step. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. A variety of large strain elastic–plastic-damage problems including severe localization is presented, and the influence of different model parameters on the deformation and localization prediction of ductile metals is discussed.     

Effective elasticity tensor of a periodic composite

The effective elasticity tensor of a composite is defined to be the four-tensor C which relates the average stress to the average strain. We determine it for an array of rigid spheres centered on the points of a periodic lattice in a homogeneous isotropic elastic medium. We first express C in terms of the traction exerted on a single sphere by the medium, and then derive an integral equation for this traction. We solve this equation numerically for simple, body-centered and face-centered cubic lattices with inclusion concentrations up to 90% of the close-packing concentration. For lattices with cubic symmetry the effective elasticity tensor involves just three parameters, which we compute from the solution for the traction. We obtain approximate asymptotic formulas for low concentrations which agree well with the numerical results. We also derive asymptotic results for C at high inclusion concentrations for arbitrary lattice geometries. We find them to be in good agreement with the numerical results for cubic lattices. For low and moderate concentrations the approximate results of Nemat-Nasseret al., also agree well with the numerical results for cubic lattices.     

Transformation elastodynamics and cloaking for flexural waves

The paper addresses an important issue of cloaking transformations for fourth-order partial differential equations representing flexural waves in thin elastic plates. It is shown that, in contrast with the Helmholtz equation, the general form of the partial differential equation is not invariant with respect to the cloaking transformation. The significant result of this paper is the analysis of the transformed equation and its interpretation in the framework of the linear theory of pre-stressed plates. The paper provides a formal framework for transformation elastodynamics as applied to elastic plates. Furthermore, an algorithm is proposed for designing a broadband square cloak for flexural waves, which employs a regularised push-out transformation. Illustrative numerical examples show high accuracy and efficiency of the proposed cloaking algorithm. In particular, a physical configuration involving a perturbation of an interference pattern generated by two coherent sources is presented. It is demonstrated that the perturbation produced by a cloaked defect is negligibly small even for such a delicate interference pattern.     

A micro-macro approach to rubber-like materials. Part III: The micro-sphere model of anisotropic Mullins-type damage

A micromechanically based non-affine network model for finite rubber elasticity and viscoelasticity was discussed in Parts I and II [Miehe, C., Göktepe, S., Lulei, F., 2004. A micro-macro approach to rubber-like materials. Part I: The non-affine micro-sphere model of rubber elasticity. J. Mech. Phys. Solids 52, 2617-2660; Miehe, C., Göktepe, S., 2005. A micro-macro approach to rubber-like materials. Part II: Viscoelasticity model for polymer networks. J. Mech. Phys. Solids, published on-line, doi:10.1016/j.jmps.2005.04.006.] of this work. In this follow-up contribution, we further extend the micro-sphere network model such that it incorporates a deformation-induced softening commonly referred to as the Mullins effect. To this end, a continuum formulation is constructed by a superimposed modeling of a crosslink-to-crosslink (CC) and a particle-to-particle (PP) network. The former is described by the non-affine elastic network model proposed in Part I. The Mullins-type damage phenomenon is embedded into the PP network and micromechanically motivated by a breakdown of bonds between chains and filler particles. Key idea of the constitutive approach is a two-step procedure that includes (i) the set up of micromechanically based constitutive models for a single chain orientation and (ii) the definition of the macroscopic stress response by a directly evaluated homogenization of state variables defined on a micro-sphere of space orientations. In contrast to previous works on the Mullins effect, our formulation inherently describes a deformation-induced anisotropy of the damage as observed in experiments. We show that the experimentally observed permanent set in stress-strain diagrams is achieved by our model in a natural way as an anisotropy effect. The performance of the model is demonstrated by means of several numerical experiments including the solution of boundary-value problems.