On Love-type waves in a finitely deformed magnetoelastic layered half-space

In this paper, the propagation of Love-type waves in a homogeneously and finitely deformed layered half-space of an incompressible non-conducting magnetoelastic material in the presence of an initial uniform magnetic field is analyzed. The equations and boundary conditions governing linearized incremental motions superimposed on an underlying deformation and magnetic field for a magnetoelastic material are summarized and then specialized to a form appropriate for the study of Love-type waves in a layered half-space. The wave propagation problem is then analyzed for different directions of the initial magnetic field for two different magnetoelastic energy functions, which are generalizations of the standard neo-Hookean and Mooney?CRivlin elasticity models. The resulting wave speed characteristics in general depend significantly on the initial magnetic field as well as on the initial finite deformation, and the results are illustrated graphically for different combinations of these parameters. In the absence of a layer, shear horizontal surface waves do not exist in a purely elastic material, but the presence of a magnetic field normal to the sagittal plane makes such waves possible, these being analogous to Bleustein?CGulyaev waves in piezoelectric materials. Such waves are discussed briefly at the end of the paper.     

A Finite Deformation Microsphere Model for Magneto-Visco-Elastic Response in Magnetorheological Elastomers

In this work we present a novel approach to the modeling of magnetorheological elastomers (MREs) for finite deformations. Keeping in mind the composite nature at the microscale, we employ the microsphere model as an effective tool to capture the constitutive response of the material. The microsphere model has been successfully applied to the modelling of rubber-like materials. Here, we extend this approach by taking into account the effect of the magnetic dipole-dipole interactions on the orientation of the polymer chains. Thus, the presented microsphere model is directly motivated by considering the underlying phenomena at the microscale level. Finally the material model is embedded in a finite element framework and the results of a boundary value problems is presented. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Purity and Quasiequality of Abelian Groups

We study the Cohn purity in an abelian group regarded as a left module over its endomorphism ring. We prove that if a finite rank torsion-free abelian group G is quasiequal to a direct sum in which all summands are purely simple modules over their endomorphism rings then the module _E(G) G is purely semisimple. This theorem makes it possible to construct abelian groups of any finite rank which are purely semisimple over their endomorphism rings and it reduces the problem of endopure semisimplicity of abelian groups to the same problem in the class of strongly indecomposable abelian groups.     

A global damping factor for a non-invasive form finding algorithm

Inverse form finding based on the finite element method (FEM) aims in determining the optimal material (undeformed) configuration when knowing the target spatial (deformed) configuration in a discretized setting. The strategy is to iteratively update the material coordinates and recompute the spatial configuration by a FEM simulation until the computed spatial nodal positions are close enough to a priori given spatial nodal positions. A form finding algorithm is utilized, which is purely based on geometrical considerations and can be coupled with arbitrary external FEM software via subroutines in a non-invasive fashion. At large deformations degenerated elements can occur when updating the material coordinates. Evaluating the mesh quality of the updated material configuration and adjusting a global damping factor before recomputing the next spatial configuration helps to avoid mesh distortions. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Limiting absorption principle for the magnetic Dirichlet Laplacian in a half-plane

We consider the Dirichlet Laplacian in the half-plane with a constant magnetic field. Due to the translational invariance, this operator admits a fiber decomposition and a family of dispersion curves that are real analytic functions. Each of them is simple and monotonically decreasing from positive infinity to a finite value, which is the corresponding Landau level. These finite limits are thresholds in the purely absolutely continuous spectrum of the magnetic Laplacian. We prove a limiting absorption principle for this operator, both outside and at the thresholds. Finally, we establish analytic and decay properties for functions lying in the absorption spaces. We point out that the analysis carried out in this article is rather general, and can be adapted to a wide class of fibered magnetic Laplacians with thresholds in their spectrum that are finite limits of their band functions.     

Characterisation of filled rubber with a pronounced non-linear viscoelasticity

This contribution presents the characterisation of an incompressible carbon black-filled elastomer as one characteristical example for highly filled rubber. It shows a strongly pronounced non-linear viscoelastic behaviour and the most important characteristic is the extremely long relaxation time which has to be taken into account. The material model is developed with respect to uniaxial tension data. The basis in the development of a phenomenological model is given by the basic elasticity. For this evaluation the long term relaxation behaviour results in a complex experimental procedure. Therefore, special attention has to be paid according to an optimised experimental process in order to get the necessary reference data in an adequate and reproduceable way [1]. With this model basis further investigations are taken into account concerning the time-dependent viscoelasticity. Therefore, cyclic deformations from zero up to a maximum of deformation are considered for different strain rates. Furthermore, the relaxation behaviour is investigated for multiple strain levels. The phenomena which are observed in the experimental results yield in a purely viscoelastic model, based on a rheological analogous model consisting of an equilibrium spring and several Maxwell-elements which contain nonlinear relations for the relaxation times of the dashpot elements [1,2]. The material model's numerical realisation is accomplished in two ways. Because of its numerical simplicity especially according to the parameter identification the model is restricted only to the simple case of uniaxial tension. A second, alternative implementation is executed providing the benefit that more complex deformation conditions can also be taken into account. Therefore, the general, three-dimensional finite model is implemented in an open-source Finite Element library [3]. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Absolutely Continuous Spectrum for Random Schrödinger Operators on Tree-Strips of Finite Cone Type

A tree-strip of finite cone type is the product of a tree of finite cone type with a finite set. We consider random Schrödinger operators on these tree-strips, similar to the Anderson model. We prove that for small disorder, the spectrum is almost surely, purely, absolutely continuous in a certain set.     

Geometrical Aspects of the Incorporation of Free Space in Magnetomechanics at Finite Strains

This paper addresses the modeling of finite magnetoelasticity with a particular focus on the incorporation of free space energy storage and magnetic body forces. First, a coordinate-invariant formulation of an extended Neo-Hookean model with magnetostrictive coupling is presented. Then, a finite element model for finite magnetomechanics based on a variational stationarity principle is formulated. The influence of magnetic body forces is studied in an exemplary boundary value problem. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite element method for viscoelastic fluids flows: approximation of the Maxwell model

We discuss about the finite element approximation of viscoelastic fluid flow. We consider a fluid obeying the Oldroyd model and particularly we study the purely viscoelastic case, the so-called Maxwell model, important in practice for the applications. We examine two kinds of methods used for the approximation of the Maxwell model: method using a splitting technique and finite element method satisfying inf-sup conditions relating tensor and velocity. We present numerical results for these methods and we discuss about their stability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.      

Topology and material orientation optimization based on evolution equations

Many modern high-performance materials have inherent anisotropic elastic properties and its local material orientation can be considered to be an additional design variable for the topology optimization [1–3]. We extend our previous model for topology optimization with variational controlled growth [4–6] for linear elastic anisotropic materials, for which the material orientation is introduced as an additional design variable. We solve the optimization problem purely with the principles of thermodynamics by minimizing the Gibbs energy. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerically efficient method for the MAP/D/1/K queue via rational approximations

The Markovian Arrival Process (MAP), which contains the Markov Modulated Poisson Process (MMPP) and the Phase-Type (PH) renewal processes as special cases, is a convenient traffic model for use in the performance analysis of Asynchronous Transfer Mode (ATM) networks. In ATM networks, packets are of fixed length and the buffering memory in switching nodes is limited to a finite numberK of cells. These motivate us to study the MAP/D/1/K queue. We present an algorithm to compute the stationary virtual waiting time distribution for the MAP/D/1/K queue via rational approximations for the deterministic service time distribution in transform domain. These approximations include the well-known Erlang distributions and the Padé approximations that we propose. Using these approximations, the solution for the queueing system is shown to reduce to the solution of a linear differential equation with suitable boundary conditions. The proposed algorithm has a computational complexity independent of the queue storage capacityK. We show through numerical examples that, the idea of using Padé approximations for the MAP/D/1/K queue can yield very high accuracy with tractable computational load even in the case of large queue capacities.This work was done when the author was with the Bilkent University, Ankara, Turkey and the research was supported by TÜBITAK under Grant No. EEEAG-93.     

A Constitutive Model for Magnetostrictive Materials - Theory and Finite Element Implementation

A 3D macroscopic constitutive law for hysteresis effects in magnetostrictive materials is presented and a finite element implementation is provided. The novel aspect of the thermodynamically consistent model is an additive decomposition of the magnetic and the strain field in a reversible and an irreversible part. Employing the irreversible magnetic field is advantageous for a finite element implementation, where the displacements and magnetic scalar potential are the nodal degrees of freedom. To consider the correlation between the irreversible magnetic field and the irreversible strains a one-to-one relation is assumed. The irreversible magnetic field determines as internal variable the movement of the center of a switching surface. This controls the motion of the domain walls during the magnetization process. The evolution of the internal variables is derived from the magnetic enthalpy function by the postulate of maximum dissipation, where the switching surface serves as constraint. The evolution equations are integrated using the backward Euler implicit integration scheme. The constitutive model is implemented in a 3D hexahedral element which provides an algorithmic consistent tangent stiffness matrix. A numerical example demonstrates the capability of the proposed model to reproduce the ferromagnetic hysteresis loops of a Terfenol-D sample. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Semi-analytic modelling of transversely isotropic magneto-electro-elastic materials under frictional sliding contact

Functionally graded magneto-electro-elastic (FGMEE) materials has been increasingly used in engineering applications, particularly in smart material or intelligent structure systems. This paper proposes a semi-analytical approach for sliding frictional contact problem between a rigid insulating sphere and a transversely isotropic FGMEE film and half-space based on frequency response functions (FRFs). Multilayered approximation is used to model the functionally graded material (FGM), and the FRFs for each MEE layer are derived explicitly. The unknown coefficients in FRFs are formulated by two matrix equations, and their efficient solution process is proposed. Based on the obtained FRFs, a highly efficient semi-analytical model (SAM) is developed which is able to solve the three-dimensional frictional contact of FGMEE materials with arbitrary layer designs. The model is validated with finite element method and the literature. Furthermore, the pressure/stress distribution and electric/magnetic potential are studied in different FGM designs to investigate the influence of material layout.     

Interpretation of parameters in phase field models for fracture

Phase field fracture models typically feature a length parameter, which controls the width of the diffuse transition zone between broken and undamaged material. In the limit case of a vanishing length parameter, these models converge to a sharp crack formulation. From this point of view, the length scale parameter is a purely auxiliary numerical quantity. However, the study of the stability of homogeneous solutions in a one dimensional setting permits a different interpretation. Since the length parameter is directly related to the critical stress at which the homogeneous solution becomes unstable and crack nucleation occurs, it can be related to the strength of the material. In this regard, the length parameter itself may be seen as a material parameter. These analytical findings are approved by finite element simulations. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stabilized nodal‐based finite element methods for the magnetic induction problem

We consider the time‐dependent magnetic induction model as a step towards the resistive magnetohydrodynamics model in incompressible media. Conforming nodal‐based finite element approximations of the induction model with inf‐sup stable finite elements for the magnetic field and the magnetic pseudo‐pressure are investigated. Based on a residual‐based stabilization technique proposed by Badia and Codina, SIAM J. Numer. Anal. 50 (2012), pp. 398–417, we consider a stabilized nodal‐based finite element method for the numerical solution. Error estimates are given for the semi‐discrete model in space. Finally, we present some examples, for example, for the magnetic flux expulsion problem, Shercliff's test case and singular solutions of the Maxwell problem. Copyright © 2015 John Wiley & Sons, Ltd.     

Longitudinal and Transverse Waves in Finite Elastic Strain. Hadamard and Green Materials

This work is concerned with the propagation of purely longitudinaland purely transverse waves in homogeneously deformed isotropicelastic materials. Two types of compressible material are alsodiscussed. A Hadamard material, so called by John in the hyperelasticcase, is one in which longitudinal waves may propagate in everydirection when the material is homogeneously deformed. A secondmaterial, called a "Green material" is introduced. In it twotransverse waves can propagate in every direction when the materialis homogeneously deformed. It is seen that a Mooney materialis the only isotropic incompressible elastic material in whichtwo transverse waves can propagate in every direction when itis homogeneously deformed, while the pressure stays constantthroughout the material. The propagation of finite amplitudewaves in these materials is discussed. Finally, it is shownthat the only motions which can be maintained in all homogeneouscompressible elastic Hadamard materials under the action ofsurface forces alone, are necessarily homogeneous and accelerationless.     

Accelerating the convergence of the method of alternating projections

The powerful von Neumann-Halperin method of alternating projections (MAP) is an algorithm for determining the best approximation to any given point in a Hilbert space from the intersection of a finite number of subspaces. It achieves this by reducing the problem to an iterative scheme which involves only computing best approximations from the individual subspaces which make up the intersection. The main practical drawback of this algorithm, at least for some applications, is that the method is slowly convergent. In this paper, we consider a general class of iterative methods which includes the MAP as a special case. For such methods, we study an ``accelerated' version of this algorithm that was considered earlier by Gubin, Polyak, and Raik (1967) and by Gearhart and Koshy (1989). We show that the accelerated algorithm converges faster than the MAP in the case of two subspaces, but is, in general, not faster than the MAP for more than two subspaces! However, for a ``symmetric' version of the MAP, the accelerated algorithm always converges faster for any number of subspaces. Our proof seems to require the use of the Spectral Theorem for selfadjoint mappings.