Authors:Artur V. Dmitrenko Pages: 1 - 9 Abstract: Abstract The stochastic equations of continuum are used for determining the hydraulic drag coefficients. As a result, the formulas for the hydraulic drag coefficients dependent on the turbulence intensity and scale instead of only on the Reynolds number are proposed for the classic flows of an incompressible fluid along a smooth flat plate and a round smooth tube. It is shown that the new expressions for the classical drag coefficients, which depend only on the Reynolds number, should be obtained from these new general formulas if to use the well-known experimental data for the initial turbulence. It is found that the limitations of classical empirical and semiempirical formulas for the hydraulic drag coefficients and their deviation from the experimental data depend on different parameters of initial fluctuations in the flow for different experiments in a wide range of Reynolds numbers. On the basis of these new dependencies, it is possible to explain that the differences between the experimental results for the fixed Reynolds number are caused by the difference in the values of flow fluctuations for each experiment instead of only due to the systematic error in the processing of experiments. Accordingly, the obtained general dependencies for the smooth flat plate and the smooth round tube can serve as the basis for clarifying the results of experiments and the experimental formulas, which used for continuum flows in different devices. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0514-1 Issue No:Vol. 29, No. 1 (2017)

Authors:Elena Bonetti; Francesco Freddi; Antonio Segatti Pages: 31 - 50 Abstract: Abstract In this paper, we consider a model describing evolution of damage in elastic materials, in which stiffness completely degenerates once the material is fully damaged. The model is written by using a phase transition approach, with respect to the damage parameter. In particular, a source of damage is represented by a quadratic form involving deformations, which vanishes in the case of complete damage. Hence, an internal constraint is ensured by a maximal monotone operator. The evolution of damage is considered “reversible”, in the sense that the material may repair itself. We can prove an existence result for a suitable weak formulation of the problem, rewritten in terms of a new variable (an internal stress). Some numerical simulations are presented in agreement with the mathematical analysis of the system. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0520-3 Issue No:Vol. 29, No. 1 (2017)

Authors:Stefan Schindler; Julia Mergheim; Marco Zimmermann; Jan C. Aurich; Paul Steinmann Pages: 51 - 75 Abstract: Abstract A two-scale material modeling approach is adopted in order to determine macroscopic thermal and elastic constitutive laws and the respective parameters for metal matrix composite (MMC). Since the common homogenization framework violates the thermodynamical consistency for non-constant temperature fields, i.e., the dissipation is not conserved through the scale transition, the respective error is calculated numerically in order to prove the applicability of the homogenization method. The thermomechanical homogenization is applied to compute the macroscopic mass density, thermal expansion, elasticity, heat capacity and thermal conductivity for two specific MMCs, i.e., aluminum alloy Al2024 reinforced with 17 or 30 % silicon carbide particles. The temperature dependency of the material properties has been considered in the range from 0 to \(500{\,}^\circ \mathrm {C}\) , the melting temperature of the alloy. The numerically determined material properties are validated with experimental data from the literature as far as possible. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0515-0 Issue No:Vol. 29, No. 1 (2017)

Authors:Lidiia Nazarenko; Swantje Bargmann; Henryk Stolarski Pages: 77 - 96 Abstract: Abstract The objective of this work is to present an approach allowing for inclusion of the complete Gurtin–Murdoch material surface equations in methods leading to closed-form formulas defining effective properties of particle-reinforced nanocomposites. Considering that all previous developments of the closed-form formulas for effective properties employ only some parts of the Gurtin–Murdoch model, its complete inclusion constitutes the main focus of this work. To this end, the recently introduced new notion of the energy-equivalent inhomogeneity is generalized to precisely include all terms of the model. The crucial aspect of that generalization is the identification of the energy associated with the last term of the Gurtin–Murdoch equation, i.e., with the surface gradient of displacements. With the help of that definition, the real nanoparticle and its surface possessing its own distinct elastic properties and residual stresses are replaced by an energy-equivalent inhomogeneity with properties incorporating all surface effects. Such equivalent inhomogeneity can then be used in combination with any existing homogenization method. In this work, the method of conditional moments is used to analyze composites with randomly dispersed spherical nanoparticles. Closed-form expressions for effective moduli are derived for both bulk and shear moduli. As numerical examples, nanoporous aluminum is investigated. The normalized bulk and shear moduli of nanoporous aluminum as a function of residual stresses are analyzed and evaluated in the context of other theoretical predictions. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0521-2 Issue No:Vol. 29, No. 1 (2017)

Authors:Diego Grandi; Ulisse Stefanelli Pages: 97 - 116 Abstract: Abstract We address a finite-plasticity model based on the symmetric tensor \(\varvec{P}^\top \! \varvec{P}\) instead of the classical plastic strain \(\varvec{P}\) . Such a structure arises by assuming that the material behavior is invariant with respect to frame transformations of the intermediate configuration. The resulting variational model is lower dimensional, symmetric and based solely on the reference configuration. We discuss the existence of energetic solutions at the material-point level as well as the convergence of time discretizations. The linearization of the model for small deformations is ascertained via a rigorous evolution- \(\Gamma \) -convergence argument. The constitutive model is combined with the equilibrium system in Part II where we prove the existence of quasistatic evolutions and ascertain the linearization limit (Grandi and Stefanelli in 2016). PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0522-1 Issue No:Vol. 29, No. 1 (2017)

Authors:Alireza Mohammadzadeh; Henning Struchtrup Pages: 117 - 144 Abstract: Abstract Heat transfer in solids is modeled by deriving the macroscopic equations for phonon transport from the phonon-Boltzmann equation. In these equations, the Callaway model with frequency-dependent relaxation time is considered to describe the Resistive and Normal processes in the phonon interactions. Also, the Brillouin zone is considered to be a sphere, and its diameter depends on the temperature of the system. A simple model to describe phonon interaction with crystal boundary is employed to obtain macroscopic boundary conditions, where the reflection kernel is the superposition of diffusive reflection, specular reflection and isotropic scattering. Macroscopic moments are defined using a polynomial of the frequency and wave vector of phonons. As an example, a system of moment equations, consisting of three directional and seven frequency moments, i.e., 63 moments in total, is used to study one-dimensional heat transfer, as well as Poiseuille flow of phonons. Our results show the importance of frequency dependency in relaxation times and macroscopic moments to predict rarefaction effects. Good agreement with data reported in the literature is obtained. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0525-y Issue No:Vol. 29, No. 1 (2017)

Authors:C. Jiang; K. Davey; R. Prosser Pages: 145 - 186 Abstract: Abstract Tessellated continuum mechanics is an approach for the representation of thermo-mechanical behaviour of porous media on tessellated continua. It involves the application of iteration function schemes using affine contraction and expansion maps, respectively, for the creation of porous fractal materials and associated tessellated continua. Highly complex geometries can be produced using a modest number of contraction mappings. The associated tessellations form the mesh in a numerical procedure. This paper tests the hypothesis that thermal analysis of porous structures can be achieved using a discontinuous Galerkin finite element method on a tessellation. Discontinuous behaviour is identified at a discontinuity network in a tessellation; its use is shown to provide a good representation of the physics relating to cellular heat exchanger designs. Results for different cellular designs (with corresponding tessellations) are contrasted against those obtained from direct analysis and very high accuracy is observed. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0523-0 Issue No:Vol. 29, No. 1 (2017)

Authors:Raimondo Penta; Alf Gerisch Pages: 187 - 206 Abstract: Abstract 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. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0526-x Issue No:Vol. 29, No. 1 (2017)

Authors:Long-Qiao Zhou; Sergey V. Meleshko Pages: 207 - 224 Abstract: Abstract A linear thermoviscoelastic model for homogeneous, aging materials with memory is established. A system of integro-differential equations is obtained by using two motions (a one-dimensional motion and a shearing motion) for this model. Applying the group analysis method to the system of integro-differential equations, the admitted Lie group is determined. Using this admitted Lie group, invariant and partially invariant solutions are found. The present paper gives a first example of application of partially invariant solutions to integro-differential equations. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0524-z Issue No:Vol. 29, No. 1 (2017)

Authors:Andrzej Ziółkowski Pages: 225 - 249 Abstract: Abstract An apparatus of micromechanics is used to isolate the key ingredients entering macroscopic Gibbs free energy function of a shape memory alloy (SMA) material. A new self-equilibrated eigenstrains influence moduli (SEIM) method is developed for consistent estimation of effective (macroscopic) thermostatic properties of solid materials, which in microscale can be regarded as amalgams of n-phase linear thermoelastic component materials with eigenstrains. The SEIM satisfy the self-consistency conditions, following from elastic reciprocity (Betti) theorem. The method allowed expressing macroscopic coherency energy and elastic complementary energy terms present in the general form of macroscopic Gibbs free energy of SMA materials in the form of semilinear and semiquadratic functions of the phase composition. Consistent SEIM estimates of elastic complementary energy, coherency energy and phase transformation strains corresponding to classical Reuss and Voigt conjectures are explicitly specified. The Voigt explicit relations served as inspiration for working out an original engineering practice-oriented semiexperimental SEIM estimates. They are especially conveniently applicable for an isotropic aggregate (composite) composed of a mixture of n isotropic phases. Using experimental data for NiTi alloy and adopting conjecture that it can be treated as an isotropic aggregate of two isotropic phases, it is shown that the NiTi coherency energy and macroscopic phase strain are practically not influenced by the difference in values of austenite and martensite elastic constants. It is shown that existence of nonzero fluctuating part of phase microeigenstrains field is responsible for building up of so-called stored energy of coherency, which is accumulated in pure martensitic phase after full completion of phase transition. Experimental data for NiTi alloy show that the stored coherency energy cannot be neglected as it considerably influences the characteristic phase transition temperatures of SMA material. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0530-1 Issue No:Vol. 29, No. 1 (2017)

Authors:Nikolai D. Tutyshkin; Paul Lofink; Wolfgang H. Müller; Ralf Wille; Oliver Stahn Pages: 251 - 269 Abstract: Abstract On the basis of the physical concepts of void formation, nucleation, and growth, generalized constitutive equations are formulated for a tensorial model of plastic damage in metals based on three invariants. The multiplicative decomposition of the metric transformation tensor and a thermodynamically consistent formulation of constitutive relations leads to a symmetric second-order damage tensor with a clear physical meaning. Its first invariant determines the damage related to plastic dilatation of the material due to growth of the voids. The second invariant of the deviatoric damage tensor is related to the change in void shape. The third invariant of the deviatoric tensor describes the impact of the stress state on damage (Lode angle), including the effect of rotating the principal axes of the stress tensor (Lode angle change). The introduction of three measures with related physical meaning allows for the description of kinetic processes of strain-induced damage with an equivalent parameter in a three-dimensional vector space, including the critical condition of ductile failure. Calculations were performed by using experimentally determined material functions for plastic dilatation and deviatoric strain at the mesoscale, as well as three-dimensional graphs for plastic damage of steel DC01. The constitutive parameter was determined from tests in tension, compression, and shear by using scanning electron microscopy, which allowed to vary the Lode angle over the full range of its values . In order to construct the three-dimensional plastic damage curve for a range of triaxiality parameters \(-1 \le ST \le 1\) and of Lode angles , we used our own, as well as systematized published experimental data. A comparison of calculations shows a significant effect of the third invariant (Lode angle) on equivalent damage. The measure of plastic damage, based on three invariants, can be useful for assessing the quality of metal mesostructure produced during metal forming processes. In many processes of metal sheet forming the material experiences, a non-proportional loading accompanied by rotating the principal axes of the stress tensor and a corresponding change of Lode angle. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0529-7 Issue No:Vol. 29, No. 1 (2017)

Authors:S. Buchkremer; F. Klocke Pages: 271 - 289 Abstract: Abstract Performance and operational safety of many metal parts in engineering depend on their surface integrity. During metal cutting, large thermomechanical loads and high gradients of the loads concerning time and location act on the surfaces and may yield significant structural material modifications, which alter the surface integrity. In this work, the derivation and validation of a model of nanostructural surface modifications in metal cutting are presented. For the first time in process modeling, initiation and kinetics of these modifications are predicted using a thermodynamic potential, which considers the interdependent developments of plastic work, dissipation, heat conduction and interface energy as well as the associated productions and flows of entropy. The potential is expressed based on the free Helmholtz energy. The irreversible thermodynamic state changes in the workpiece surface are homogenized over the volume in order to bridge the gap between discrete phenomena involved with the initiation and kinetics of dynamic recrystallization and its macroscopic implications for surface integrity. The formulation of the thermodynamic potential is implemented into a finite element model of orthogonal cutting of steel AISI 4140. Close agreement is achieved between predicted nanostructures and those obtained in transmission electron microscopical investigations of specimen produced in cutting experiments. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0533-y Issue No:Vol. 29, No. 1 (2017)

Authors:Philipp Junker; Stephan Schwarz; Jerzy Makowski; Klaus Hackl Pages: 291 - 310 Abstract: Abstract Material models, including softening effects due to, for example, damage and localizations, share the problem of ill-posed boundary value problems that yield mesh-dependent finite element results. It is thus necessary to apply regularization techniques that couple local behavior described, for example, by internal variables, at a spatial level. This can take account of the gradient of the internal variable to yield mesh-independent finite element results. In this paper, we present a new approach to damage modeling that does not use common field functions, inclusion of gradients or complex integration techniques: Appropriate modifications of the relaxed (condensed) energy hold the same advantage as other methods, but with much less numerical effort. We start with the theoretical derivation and then discuss the numerical treatment. Finally, we present finite element results that prove empirically how the new approach works. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0528-8 Issue No:Vol. 29, No. 1 (2017)

Authors:Julian Heß; Yongqi Wang; Kolumban Hutter Pages: 311 - 343 Abstract: Abstract This paper presents a new, thermodynamically consistent model for granular-fluid mixtures, derived with the entropy principle of Müller and Liu. Including a pressure diffusion equation combined with the concept of extra pore pressure, and hypoplastic material behavior, thermodynamic restrictions are imposed on the constitutive quantities. The model is applied to a granular-fluid flow, using a closing assumption in conjunction with the fluid pressure. While the focal point of the work is the conceptional part, i.e. the thermodynamic consistent modeling, numerical simulations with physically reasonable results for simple shear flow are also carried out. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0535-9 Issue No:Vol. 29, No. 1 (2017)

Authors:Daria Andreeva; Wiktoria Miszuris Pages: 345 - 358 Abstract: Abstract In this paper, we consider the heat transfer problem in a cylindrical composite material with an adhesive interphase of closed curvilinear shape. The interphase exhibits nonlinear behaviour and at the same time has physical characteristics (size, thermal conductivity) that significantly vary from the properties of the surrounding components. As the latter complicates the direct use of FEM, we use asymptotic method to replace the interphase with an imperfect interface of zero thickness, preserving the essential features of the thermal behaviour of the interphase through the evaluated transmission conditions. Carefully designed numerical simulations verify their validity. We place special emphasis on the impact of geometric properties of the interphase, in particular, the curvature of its boundaries, on the accuracy of the conditions. PubDate: 2017-01-01 DOI: 10.1007/s00161-016-0532-z Issue No:Vol. 29, No. 1 (2017)

Authors:Tomáš Roubíček Abstract: Abstract A rather general model for fluid and heat transport in poro-elastic continua undergoing possibly also plastic-like deformation and damage is developed with the goal to cover various specific models of rock rheology used in geophysics of Earth’s crust. Nonconvex free energy at small elastic strains, gradient theories (in particular the concept of second-grade nonsimple continua), and Biot poro-elastic model are employed, together with possible large displacement due to large plastic-like strains evolving during long time periods. Also the additive splitting is justified in stratified situations which are of interest in modelling of lithospheric crust faults. Thermodynamically based formulation includes entropy balance (in particular the Clausius–Duhem inequality) and an explicit global energy balance. It is further outlined that the energy balance can be used to ensure, under suitable data qualification, existence of a weak solution and stability and convergence of suitable approximation schemes at least in some particular situations. PubDate: 2017-01-13 DOI: 10.1007/s00161-016-0547-5

Authors:A. Paglietti Abstract: Abstract In laminar flow, viscous fluids must exert appropriate elastic shear stresses normal to the flow direction. This is a direct consequence of the balance of angular momentum. There is a limit, however, to the maximum elastic shear stress that a fluid can exert. This is the ultimate shear stress, \(\tau _\mathrm{y}\) , of the fluid. If this limit is exceeded, laminar flow becomes dynamically incompatible. The ultimate shear stress of a fluid can be determined from experiments on plane Couette flow. For water at \(20\,^{\circ }\hbox {C}\) , the data available in the literature indicate a value of \(\tau _\mathrm{y}\) of about \(14.4\times 10^{-3}\, \hbox {Pa}\) . This study applies this value to determine the Reynolds numbers at which flowing water reaches its ultimate shear stress in the case of Taylor–Couette flow and circular pipe flow. The Reynolds numbers thus obtained turn out to be reasonably close to those corresponding to the onset of turbulence in the considered flows. This suggests a connection between the limit to laminar flow, on the one hand, and the occurrence of turbulence, on the other. PubDate: 2017-01-10 DOI: 10.1007/s00161-016-0549-3

Authors:P. Ván Abstract: Abstract Single-component nonrelativistic dissipative fluids are treated independently of reference frames and flow-frames. First the basic fields and their balances are derived, then the related thermodynamic relations and the entropy production are calculated and the linear constitutive relations are given. The usual basic fields of mass, momentum, energy and their current densities, the heat flux, pressure tensor and diffusion flux are the time- and spacelike components of the third-order mass–momentum–energy density-flux four-tensor. The corresponding Galilean transformation rules of the physical quantities are derived. It is proved that the non-equilibrium thermodynamic frame theory, including the thermostatic Gibbs relation and extensivity condition and also the entropy production, is independent of the reference frame and also the flow-frame of the fluid. The continuity-Fourier–Navier–Stokes equations are obtained almost in the traditional form if the flow of the fluid is fixed to the temperature. This choice of the flow-frame is the thermo-flow. A simple consequence of the theory is that the relation between the total, kinetic and internal energies is a Galilean transformation rule. PubDate: 2017-01-09 DOI: 10.1007/s00161-016-0545-7

Authors:Teodor M. Atanacković; Marko Janev; Sanja Konjik; Stevan Pilipović Abstract: Abstract We study waves in a viscoelastic rod whose constitutive equation is of generalized Zener type that contains fractional derivatives of complex order. The restrictions following from the Second Law of Thermodynamics are derived. The initial boundary value problem for such materials is formulated and solution is presented in the form of convolution. Two specific examples are analyzed. PubDate: 2017-01-06 DOI: 10.1007/s00161-016-0548-4