Thermodynamic and statistical principles of the theory of elastoviscoplastic deformation and strengthening of materials

We propose a thermodynamic method and a statistical one for constructing the constitutive equations of elastoviscoplastic deformation and strengthening of materials. The thermodynamic method is based on the energy conservation law as well as the equations of entropy balance and entropy generation in the presence of self-equilibrated internal microstresses, which are characterized by coupled strengthening parameters. The general constitutive equations consist of the relations between thermodynamic flows and forces, which follow from nonnegativity of entropy generation and satisfy the generalized Onsager principle, as well as the thermoelasticity relations and the expression for entropy, which follow from the energy conservation law. The specific constitutive equations are obtained on the basis of representation of the energy dissipation rate as a sum of two constituents that describe translational and isotropic strengthening and are approximated by power and hyperbolic sine laws. Starting from the stochastic microstructural concepts, we construct the constitutive equations of elastoviscoplastic deformation and strengthening on the basis of the linear model of thermoelasticity and the nonlinear Maxwell model for spherical and deviatoric components of microstresses and microstrains, respectively. The solution of the problem of the effective properties and stress-strain state of a three-component material is constructed with the use of the combined Voigt–Reuss scheme and leads to constitutive equations coinciding, as to their form, with similar equations constructed by the thermodynamic method.     

On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework

We provide a thermodynamic basis for the development of models that are usually referred to as ??phase-field models?? for compressible, incompressible, and quasi-incompressible fluids. Using the theory of mixtures as a starting point, we develop a framework within which we can derive ??phase-field models?? both for mixtures of two constituents and for mixtures of arbitrarily many fluids. In order to obtain the constitutive equations, we appeal to the requirement that among all admissible constitutive relations that which is appropriate maximizes the rate of entropy production (see Rajagopal and Srinivasa in Proc R Soc Lond A 460:631?C651, 2004). The procedure has the advantage that the theory is based on prescribing the constitutive equations for only two scalars: the entropy and the entropy production. Unlike the assumption made in the case of the Navier?CStokes?CFourier fluids, we suppose that the entropy is not only a function of the internal energy and the density but also of gradients of the partial densities or the concentration gradients. The form for the rate of entropy production is the same as that for the Navier?CStokes?CFourier fluid. As observed earlier in Heida and Málek (Int J Eng Sci 48(11):1313?C1324, 2010), it turns out that the dependence of the rate of entropy production on the thermodynamical fluxes is crucial. The resulting equations are of the Cahn?CHilliard?CNavier?CStokes type and can be expressed both in terms of density gradients or concentration gradients. As particular cases, we will obtain the Cahn?CHilliard?CNavier?CStokes system as well as the Korteweg equation. Compared to earlier approaches, our methodology has the advantage that it directly takes into account the rate of entropy production and can take into consideration any constitutive assumption for the internal energy (or entropy).     

An approach to constructing models of multiphase elastic porous media

A novel approach to describing the behaviour of multiphase elastic porous media is proposed. The average values of the physical quantities needed to describe the motions of porous media are formulated using an integral relation. The validity of this relation is taken as the fundamental hypothesis. The integral definition of the average values enables integral relations to be devised for the average values from the integral laws of conservation of mass, momentum and energy and the increase in entropy. Along with the average values, the integral relations contain new variables that can be identified with generalized thermodynamic forces, which can be used to take into account the phase interaction in a porous medium. The integral relations are used to derive differential equations for the rate of entropy change and Gibbs relations for a porous medium as a basis for obtaining the constitutive relations. Relationships between the thermomechanical parameters of the model are established from the Gibbs relations under additional assumptions. The equation for the rate of entropy change can be used to establish relations between the generalized thermodynamic forces and fluxes. A complete system of differential equations in the defining parameters, which describes the motion of multiphase elastic porous media, is finally obtained.     

Global energy‐tracking identities and global energy‐tracking splitting FDTD schemes for the Drude Models of Maxwell's equations in three‐dimensional metamaterials

In this article, we study the Drude models of Maxwell's equations in three‐dimensional metamaterials. We derive new global energy‐tracking identities for the three dimensional electromagnetic problems in the Drude metamaterials, which describe the invariance of global electromagnetic energy in variation forms. We propose the time second‐order global energy‐tracking splitting FDTD schemes for the Drude model in three dimensions. The significant feature is that the developed schemes are global energy‐preserving, unconditionally stable, second‐order accurate both in time and space, and computationally efficient. We rigorously prove that the new schemes satisfy these energy‐tracking identities in the discrete form and the discrete variation form and are unconditionally stable. We prove that the schemes in metamaterials are second order both in time and space. The superconvergence of the schemes in the discrete H¹ norm is further obtained to be second order both in time and space. Their approximations of divergence‐free are also analyzed to have second‐order accuracy both in time and space. Numerical experiments confirm our theoretical analysis results. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 763–785, 2017     

A thermodynamical framework for chemically reacting systems

In this paper, we develop a thermodynamic framework that is capable of describing the response of viscoelastic materials that are undergoing chemical reactions that takes into account stoichiometry. Of course, as a special sub-case, we can also describe the response of elastic materials that undergo chemical reactions. The study generalizes the framework developed by Rajagopal and co-workers to study the response of a disparate class of bodies undergoing entropy producing processes. One of the quintessential feature of this framework is that the second law of thermodynamics is formulated by introducing Gibbs?? potential, which is the natural way to study problems involving chemical reactions. The Gibbs potential?Cbased formulation also naturally leads to implicit constitutive equations for the stress tensor. Another feature of the framework is that the constraints due to stoichiometry can also be taken into account in a consistent manner. The assumption of maximization of the rate of entropy production due to dissipation, heat conduction, and chemical reactions is invoked to determine an equation for the evolution of the natural configuration ?? _p(t)(B), the heat flux vector and a novel set of equations for the evolution of the concentration of the chemical constituents. To determine the efficacy of the framework with regard to chemical reactions, those occurring during vulcanization, a challenging set of chemical reactions, are chosen. More than one type of reaction mechanism is considered and the theoretically predicted distribution of mono, di and polysulfidic cross-links agree reasonably well with available experimental data.     

Relaxation‐time limit of the three‐dimensional hydrodynamic model with boundary effects

In this paper, we study three‐dimensional (3D) unipolar and bipolar hydrodynamic models and corresponding drift‐diffusion models from semiconductor devices on bounded domain. Based on the asymptotic behavior of the solutions to the initial boundary value problems with slip boundary condition, we investigate the relation between the 3D hydrodynamic semiconductor models and the corresponding drift‐diffusion models. That is, we discuss the relation‐time limit from the 3D hydrodynamic semiconductor models to the corresponding drift‐diffusion models by comparing the large‐time behavior of these two models. These results can be showed by energy arguments. Copyrightcopyright 2011 John Wiley & Sons, Ltd.     

Three‐dimensional blockmodels

In the classic blockmodel formulation, a social network among members of a population with n actors and k relations (types of tie) is arrayed as k n X n matrices. Though this is a three‐dimensional data structure, it is typically reduced to a two‐way analysis. In this paper, a three‐way procedure for analyzing multigraph data is developed. Specifically, in addition to applying the rule of structural equivalence to collapse actors, it is also applied to the relations (the third dimension), and structurally equivalent sets of relations are collapsed. The result is a three‐dimensional blockmodel (image) of social structure that is a more parsimonious representation of social structure than the classic two‐dimensional blockmodel images. The three‐dimensional approach is illustrated by application to three case studies: Homan's Bank Wiring Room, Sampson's monastery, and a local economy of hospital services. The structural equivalence approach to relations is further explored by applying it to the individual‐level Liking and Antagonism relations and their compounds (of length four or less) in the Bank Wiring Room. This application demonstrates that the structural equivalence approach can be used to identify equality equations for primitive and compound relations.     

An experimentally validated rod model for soft continuum robots

Soft robots are bio-inspired, highly deformable robots with the ability to interact with workpieces in a manner that complements their hard robot counterparts. To develop practical applications and reproducible designs of soft robots, new models are necessary to describe their kinematics and dynamics. In the present work, we describe experimental and numerical investigations of a popular pneumatically-actuated soft continuum arm. These works enable us to derive constitutive relations and develop a rod model for large deformations of the arm that faithfully describes its bending behavior. We show how the resulting non-classical constitutive relation can be defined either through experiments or through quasi-static finite element simulations. With the help of this relation, the resulting rod model can be used to study the dynamics of the soft robot arm in a fast and tractable manner. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three dimensional linear motion transformation in a higher order dimensional space

The objective of the work presented in this paper is an attempt at solving and transforming of the known from the classical mechanics three dimensional – single mass mathematical and mechanical vibration models in a higher order dimensional space with any virtual sectional curvature – positive or negative, constant or variable. The object of the investigation is a class of three dimensional surfaces. The aims of the work presented in this paper are to illustrate the performance of the common algorithm in three dimensional linear motion transformation, that means to transform 3D space in a higher order dimensional space and a comparison is derived on the behavior of the common algorithm depending on the surface properties. A characterization of the Riemannian Manifolds is performed by means of curvature operators in the three dimensional solution. The computer codes Mathematica and MATLAB are used in the numerical simulation. The system motion is investigated in a 3-D qualitative aspect in time and frequency domain. The application can be in topology when geodesists make snap shots of the surface profile, then the curved lines can be analyzed and transformed in the desired space dimension. Any kind of a trajectory of motion can be transformed successfully in a higher order dimensional space and vice verse by means of applying of the common algorithm.     

The distance-decay function of geographical gravity model: Power law or exponential law?

The distance-decay function of the geographical gravity model is originally an inverse power law, which suggests a scaling process in spatial interaction. However, the distance exponent of the model cannot be reasonably explained with the ideas from Euclidean geometry. This results in a dimension dilemma in geographical analysis. Consequently, a negative exponential function was used to replace the inverse power function to serve for a distance-decay function. But a new puzzle arose that the exponential-based gravity model goes against the first law of geography. This paper is devoted for solving these kinds of problems by mathematical reasoning and empirical analysis. New findings are as follows. First, the distance exponent of the gravity model is demonstrated to be a fractal dimension using the geometric measure relation. Second, the similarities and differences between the gravity models and spatial interaction models are revealed using allometric relations. Third, a four-parameter gravity model possesses a symmetrical expression, and we need dual gravity models to describe spatial flows. The observational data of China's cities and regions (29 elements indicative of 841 data points) in 2010 are employed to verify the theoretical inferences. A conclusion can be reached that the geographical gravity model based on power-law decay is more suitable for analyzing large, complex, and scale-free regional and urban systems. This study lends further support to the suggestion that the underlying rationale of fractal structure is entropy maximization. Moreover, it suggests that many dimensional dilemmas of spatial modeling can be solved using the concepts from fractal geometry.     

Mixing Closures for Conservation Laws in Stratified Flows

A closure for shocks involving the mixing of the fluids in two-layer stratified flows is proposed. The closure maximizes the rate of mixing, treating the dynamical hydraulic equations and entropy conditions as constraints. This closure may also be viewed as yielding an upper bound on the mixing rate by internal shocks. It is shown that the maximal mixing rate is accomplished by a shock moving at the fastest allowable speed against the upstream flow. Depending on whether the active constraint limiting this speed is the Lax entropy condition or the positive dissipation of energy, we distinguish precisely between internal hydraulic jumps and bores. Maximizing entrainment is shown to be equivalent to maximizing a suitable entropy associated to mixing. By using the latter, one can describe the flow globally by an optimization procedure, without treating the shocks separately. A general mathematical framework is formulated that can be applied whenever an insufficient number of conservation laws is supplemented by a maximization principle.     

Theoretical study and numerical simulation of textiles

We propose a new approach for developing continuum models fit to describe the mechanical behavior of textiles. We develop a physically motivated model, based on the properties of the yarns, which can predict and simulate the textile behavior. The approach relies on the selection of a suitable topological model for the patch of the textile, coupled with constitutive models for the yarn behavior. The textile structural configuration is related to the deformation through an energy functional, which depends on both the macroscopic deformation and the distribution of internal nodes. We determine the equilibrium positions of these latter, constrained to an assigned macroscopic deformation. As a result, we derive a macroscopic strain energy function, which reflects the possibly nonlinear character of the yarns as well as the anisotropy induced by the microscopic topological pattern. By means of both analytical estimates and numerical experiments, we show that our model is well suited for both academic test cases and real industrial textiles, with particular emphasis on the tricot textile.     

Compressible generalized Newtonian fluids

We systematically derive models that would be suitable to describe flows of compressible fluids with the material moduli depending on the symmetric part of the velocity gradient and temperature, within the context of a thermodynamic framework that has been quite successful in developing models to describe the response of bodies that produce entropy while undergoing processes.     

A variety of (3 + 1)‐dimensional Burgers equations derived by using the Burgers recursion operator

A new variety of (3 + 1)‐dimensional Burgers equations is presented. The recursion operator of the Burgers equation is employed to establish these higher‐dimensional integrable models. A generalized dispersion relation and a generalized form for the one kink solutions is developed. The new equations generate distinct solitons structures and distinct dispersion relations as well. Copyright © 2014 John Wiley & Sons, Ltd.     

The homogenization of the three dimensional skew polynomial algebras of type I

This paper studies a special case of graded central extensions of three dimen-sional Artin-Schelter regular algebras, see [9, §3]. The algebras are homoge-nizations of two classes of three dimensional skew polynomial algebras. We refer to these algebras as Type I and Type II algebras. We describe the non-commutative projective geometry and compute the finite dimensional simple modules for the homogenization of Type I algebras in the case that α is not a primitive root of unity. In this case, all finite dimensional simple modules are quotients of line modules that are homogenizations of Verma modules.     

On the constitutive relations for isotropic and transversely isotropic materials

In this work the strain and stress spaces constitutive relations for isotropic and transversely isotropic softening materials are developed. The loading surface is considered in the strain space and the normality rule; the stress relaxation is proportional to the gradient of the loading surface, is adopted. It is found that the strain space plasticity theory allows us to describe the hardening, perfectly plastic and softening materials more accurately. The validity of the strain space constitutive relation for transversely isotropic materials are confirmed by comparing with the experimental data for fiber reinforced composite materials. Some numerical examples in two and three dimensional elasto-plastic problems for various loading–unloading conditions are presented, and give a very good agreement with the existing results.     

On the stability of the upper airways system

Simple models for upper pharyngeal obstruction, describing the sleep apnea syndrome are proposed. Stability is discussed, of two and three individualized elements, with and without elastic connections, interacting with the steady flow. Considering the flow as the controlling parameter, critical steady state flows are located and their post-critical behavior is discussed for various models. It is pointed out that non-linear constitutive elastic laws are necessary contrary to the linear ones introduced by Fodil (1998) [15]. Finally the three element model will be presented and studied with non-linear constitutive relations and side connections. Applications of the theory will be performed and discussion for the three models will be presented. It is pointed out that the sleep apnea syndrome is due to the instability of the upper pharyngeal region.     

An unbounded Sum-of-Squares hierarchy integrality gap for a polynomially solvable problem

In this expository article, we study optimization problems specified via linear and relative entropy inequalities. Such relative entropy programs (REPs) are convex optimization problems as the relative entropy function is jointly convex with respect to both its arguments. Prominent families of convex programs such as geometric programs (GPs), second-order cone programs, and entropy maximization problems are special cases of REPs, although REPs are more general than these classes of problems. We provide solutions based on REPs to a range of problems such as permanent maximization, robust optimization formulations of GPs, and hitting-time estimation in dynamical systems. We survey previous approaches to some of these problems and the limitations of those methods, and we highlight the more powerful generalizations afforded by REPs. We conclude with a discussion of quantum analogs of the relative entropy function, including a review of the similarities and distinctions with respect to the classical case. We also describe a stylized application of quantum relative entropy optimization that exploits the joint convexity of the quantum relative entropy function.     

Modeling of nonlinear viscoelastic contact problems with large deformations

Contact problems are one of the most important engineering problems. These problems become much more tedious when one of the contacting bodies behaves nonlinear viscoelasticity and large deformations. This paper presents an incremental-iterative finite element model for the analysis of two dimensional quasistatic frictionless contact problems. Nonlinear viscoelastic behavior and large deformations are considered. The Schapery’s single-integral creep model with stress-dependent properties is used for nonlinear viscoelasticity. The constitutive equations are transformed into an incremental form resulting in a recursive relationship. Thereby, the need to store the entire strain histories is eliminated, except that from the previous time increment. The updated Lagrangian formulation is used to model the material and geometrical nonlinearities. Also, the Lagrange multiplier method is adopted to enforce the contact constraints. The converged solution is obtained using the Newton–Raphson iterative technique. The developed model has been verified with the previously published works and found a good agreement with them. To demonstrate the efficient capability of the developed computational model, three contact problems with different nature are analyzed.     

Infinite‐dimensional exponential attractors for fourth‐order nonlinear parabolic equations in unbounded domains

We study in this article the long‐time behavior of solutions of fourth‐order parabolic equations in bfR³. In particular, we prove that under appropriate assumptions on the nonlinear interaction function and on the external forces, these equations possess infinite‐dimensional exponential attractors whose Kolmogorov's ε‐entropy satisfies an estimate of the same type as that obtained previously for the ε‐entropy of the global attractor. Copyright © 2011 John Wiley & Sons, Ltd.