A multiscale-multiphysics strategy for numerical modeling of thin piezoelectric sheets

Flexible piezoelectric devices made of polymeric materials are widely used for micro- and nano-electro-mechanical systems. In particular, numerous recent applications concern energy harvesting. Due to the importance of computational modeling to understand the influence that microscale geometry and constitutive variables exert on the macroscopic behavior, a numerical approach is developed here for multiscale and multiphysics modeling of piezoelectric materials made of aligned arrays of polymeric nanofibers. At the microscale, the representative volume element consists in piezoelectric polymeric nanofibers, assumed to feature a linear piezoelastic constitutive behavior and subjected to electromechanical contact constraints using the penalty method. To avoid the drawbacks associated with the non-smooth discretization of the master surface, a contact smoothing approach based on B\'ezier patches is extended to the multiphysics framework providing an improved continuity of the parameterization. The contact element contributions to the virtual work equations are included through suitable electric, mechanical and coupling potentials. From the solution of the micro-scale boundary value problem, a suitable scale transition procedure leads to the formulation of a macroscopic thin piezoelectric shell element.Comment: 11 pages, 6 pages, 21 reference

Multiscale computational first order homogenization of thick shells for the analysis of out-of-plane loaded masonry walls

This work presents a multiscale method based on computational homogenization for the analysis of general heterogeneous thick shell structures, with special focus on periodic brick-masonry walls. The proposed method is designed for the analysis of shells whose micro-structure is heterogeneous in the in-plane directions, but initially homogeneous in the shell-thickness direction, a structural topology that can be found in single-leaf brick masonry walls. Under this assumption, this work proposes an efficient homogenization scheme where both the macro-scale and the micro-scale are described by the same shell theory. The proposed method is then applied to the analysis of out-of-plane loaded brick-masonry walls, and compared to experimental and micro-modeling results.Peer ReviewedPostprint (author's final draft

Multiscale Finite Element Modeling of Nonlinear Magnetoquasistatic Problems Using Magnetic Induction Conforming Formulations

In this paper we develop magnetic induction conforming multiscale formulations for magnetoquasistatic problems involving periodic materials. The formulations are derived using the periodic homogenization theory and applied within a heterogeneous multiscale approach. Therefore the fine-scale problem is replaced by a macroscale problem defined on a coarse mesh that covers the entire domain and many mesoscale problems defined on finely-meshed small areas around some points of interest of the macroscale mesh (e.g. numerical quadrature points). The exchange of information between these macro and meso problems is thoroughly explained in this paper. For the sake of validation, we consider a two-dimensional geometry of an idealized periodic soft magnetic composite.Comment: Paper accepted for publication in the SIAM MMS journa

The asymptotic homogenization elasticity tensor properties for composites with material discontinuities

The classical asymptotic homogenization approach for linear elastic composites with discontinuous material properties is considered as a starting point. The sharp length scale separation between the fine periodic structure and the whole material formally leads to anisotropic elastic-type balance equations on the coarse scale, where the arising fourth rank operator is to be computed solving single periodic cell problems on the fine scale. After revisiting the derivation of the problem, which here explicitly points out how the discontinuity in the individual constituents’ elastic coefficients translates into stress jump interface conditions for the cell problems, we prove that the gradient of the cell problem solution is minor symmetric and that its cell average is zero. This property holds for perfect interfaces only (i.e., when the elastic displacement is continuous across the composite’s interface) and can be used to assess the accuracy of the computed numerical solutions. These facts are further exploited, together with the individual constituents’ elastic coefficients and the specific form of the cell problems, to prove a theorem that characterizes the fourth rank operator appearing in the coarse-scale elastic-type balance equations as a composite material effective elasticity tensor. We both recover known facts, such as minor and major symmetries and positive definiteness, and establish new facts concerning the Voigt and Reuss bounds. The latter are shown for the first time without assuming any equivalence between coarse and fine-scale energies (Hill’s condition), which, in contrast to the case of representative volume elements, does not identically hold in the context of asymptotic homogenization. We conclude with instructive three-dimensional numerical simulations of a soft elastic matrix with an embedded cubic stiffer inclusion to show the profile of the physically relevant elastic moduli (Young’s and shear moduli) and Poisson’s ratio at increasing (up to 100 %) inclusion’s volume fraction, thus providing a proxy for the design of artificial elastic composites