No abstract

No abstract

SUMMARYThis paper presents a general framework for the macroscopic, continuum-based formulation and numerical implementation of dissipative functional materials with electro-magneto-mechanical couplings based on incremental variational principles. We focus on quasi-static problems, where mechanical inertia effects and time-dependent electro-magnetic couplings are a priori neglected and a time-dependence enters the formulation only through a possible rate-dependent dissipative material response. The underlying variational structure of non-reversible coupled processes is related to a canonical constitutive modeling approach, often addressed to so-called standard dissipative materials. It is shown to have enormous consequences with respect to all aspects of the continuum-based modeling in macroscopic electro-magnetomechanics. At first, the local constitutive modeling of the coupled dissipative response, i.e. stress, electric and magnetic fields versus strain, electric displacement and magnetic induction, is shown to be variational based, governed by incremental minimization and saddle-point principles. Next, the implications on the formulation of boundary-value problems are addressed, which appear in energy-based formulations as minimization principles and in enthalpy-based formulations in the form of saddle-point principles. Furthermore, the material stability of dissipative electro-magneto-mechanics on the macroscopic level is defined based on the convexity/concavity of incremental potentials. We provide a comprehensive outline of alternative variational structures and discuss details of their computational implementation, such as formulation of constitutive update algorithms and finite element solvers. From the viewpoint of constitutive modeling, including the understanding of the stability in coupled electro-magneto-mechanics, an energy-based formulation is shown to be the canonical setting. From the viewpoint of the computational convenience, an enthalpy-based formulation is the most convenient setting. A numerical investigation of a multiferroic composite demonstrates perspectives of the proposed framework with regard to the future design of new functional materials.

This article is concerned with the modeling of the magnetic shape memory alloy (MSMA) constitutive response caused by the reorientation of martensitic variants under mechanical and magnetic fields. The presented model is able to better capture the complexity of the magnetomechanical MSMA behavior by accounting not only for the mechanism of field-induced variant reorientation, but also the magnetization rotation away from magnetic easy axes and the magnetic domain wall motion at low stress and magnetic field levels. Following the general formulation of the model, reduced versions of the constitutive equations are derived for three specific loading cases: (1) magnetic-field-induced variant reorientation at constant stress; (2) stress-induced variant reorientation at constant magnetic field; (3) variant reorientation under collinear magnetic field and stress with perpendicular bias field. For each of these cases the nonlinear and hysteretic strain and magnetization response of MSMAs are predicted and compared to experimental data where available. The relation of critical stresses and magnetic fields for the activation of the reorientation process are visualized in a variant reorientation diagram. The captured loading-history-dependent macroscopic material response is explained in detail by connecting it to the evolution of the crystallographic and magnetic microstructure as represented by a set of internal state variables.

A non-local gradient-enhanced damage-plasticity formulation is proposed, which prevents the loss of well-posedness of the governing field equations in the post-critical damage regime. The non-locality of the formulation then manifests itself in terms of a non-local free energy contribution that penalizes the occurrence of damage gradients. A second penalty term is introduced to force the global damage field to coincide with the internal damage state variable at the Gauss point level. An enforcement of Karush–Kuhn–Tucker conditions on the global level can thus be avoided and classical local damage models may directly be incorporated and equipped with a non-local gradient enhancement. An important part of the present work is to investigate the efficiency and robustness of different algorithmic schemes to locally enforce the Karush–Kuhn–Tucker conditions in the multi-surface damage-plasticity setting. Response simulations for representative inhomogeneous boundary value problems are studied to assess the effectiveness of the gradient enhancement regarding stability and mesh objectivity.

No abstract

From a mathematical point of view, phase field theory can be understood as a smooth approximation of an underlying sharp interface problem. However, the smooth phase field approximation is not uniquely defined. Different phase field approximations are known to converge to the same sharp interface problem in the limiting case-if the thickness of the diffuse interface converges to zero. In this respect and focusing on numerics, a question that naturally arises is as follows: What are the convergence rates of the different phase field models? The paper deals precisely with this question for a certain family of phase field models. The focus is on an Allen-Cahn-type phase field model coupled to continuum mechanics. This model is rewritten into a unified variational phase field framework that covers different homogenization assumptions in the diffuse interfaces: Voigt/Taylor, Reuss/Sachs and more sound homogenization approaches falling into the range of rank-one convexification. It is shown by means of numerical experiments that the underlying phase field model-that is, the homogenization assumption in the diffuse interface-indeed influences the convergence rate.see [3], is comprised of a double-well contribution in the first term, defining distinct energetic minima for two phases-at order-parameter field values of p = 0 and p = 1, respectively-and the second, gradient-based term that attributes an energetic cost to the formation of interfaces. Assuming (2), the Euler -Lagrange equation of (1) can be written in the form, see also [19],NUMERICAL CONVERGENCE STUDY IN PHASE FIELD MODELING ) denotes the functional (variational) derivative. Physically speaking, this describes an order-parameter field that evolves in such a manner that the system is driven towards an energy minimizing equilibrium

No abstract