On thermodynamic potentials in linear thermoelasticity

The four thermodynamic potentials, the internal energy u=u(ε_ij,s), the Helmholtz free energy f=f(ε_ij,T), the Gibbs energy g=g(σ_ij,T) and the enthalpy h=h(σ_ij,s) are derived, independently of each other, by using the Duhamel–Neumann extension of Hooke's law and an assumed linear dependence of the specific heat on temperature. A systematic procedure is then presented to express all thermodynamic potentials in terms of four possible pairs of independent state variables. This procedure circumvents a tedious transition from one potential to another, based on the formal change of variables, and inversions of the stress–strain and entropy–temperature relations. The general results are applied to uniaxial loading paths under isothermal, adiabatic, constant stress, and constant strain conditions. An interplay of adiabatic and isothermal elastic constants in the expressions for exchanged heat along certain thermodynamic paths is indicated.     

Thermomechanical material modelling based on a hybrid free energy density depending on pressure,isochoric deformation and temperature

In order to represent temperature-dependent mechanical material properties in a thermomechanical consistent manner it is common practice to start with the definition of a model for the specific Helmholtz free energy. Its canonical independent variables are the Green strain tensor and the temperature. But to represent calorimetric material properties under isobaric conditions, for example the exothermal behaviour of a curing process or the dependence of the specific heat on the temperature history, the temperature and the pressure should be taken as independent variables. Thus, in the field of calorimetry the Gibbs free energy is usually used as thermodynamic potential whereas in continuum mechanics the Helmholtz free energy is normally applied. In order to simplify the representation of calorimetric phenomena in continuum mechanics a hybrid free energy density is introduced. Its canonical independent variables are the isochoric Green strain tensor, the pressure and the temperature. It is related to the Helmholtz free energy density by a Legendre transformation. In combination with the additive split of the stress power into the sum of isochoric and volumetric terms this approach leads to thermomechanical consistent constitutive models for large deformations. The article closes with applications of this approach to finite thermoelasticity, curing adhesives and the glass transition.     

Il'iushin's postulate and resulting thermodynamic conditions on elasto-plastic coupling

Il'iushin's postulate is restated within a general thermodynamic strain space formulation of rate independent plasticity by means of plastic internal variables. This yields a general expression in terms of appropriate thermodynamic potentials. A combination of a thermodynamic condition, derived from the general development, with the results of Il'iushin's postulate, furnishes explicit conditions on elasto-plastic coupling. A specific example is presented, with the plastic work being the only plastic internal variable. Necessary and sufficient consitbns on the elastic moduli and their change with plastic deformation are derived, for the thermodynamic condition to be satisfied.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

Mean-field homogenization of elasto-viscoplastic composites based on a general incrementally affine linearization method

We propose a general formulation – which we believe to be new – for the mean-field homogenization of inclusion-reinforced elasto-viscoplastic composites assuming small strains. Our proposal is based on an interplay between constitutive equations and numerical algorithms, and the key ideas behind it are the following. The evolution equations for inelastic strain and internal variables at the beginning of each time interval are linearized around the ending time of the same interval. The linearized equations are then numerically integrated using a fully implicit backward Euler scheme. The obtained algebraic equations lead to an incrementally affine stress–strain relation which involves two important terms. The first one is the algorithmic tangent operator, obtained by consistent linearization of the time discretized constitutive equations. The second term is a new one and called an affine strain increment. The proposal leads to thermoelastic-like relations directly in the time domain, and not in the Laplace–Carson (L–C) one. There is no need for viscoelastic-type integral rewriting of the evolution equations, for L–C transforms, or for numerical inversion back from L–C to time domains. The proposed method can be readily applied to sophisticated elasto-viscoplastic models with an arbitrary set of scalar or tensor internal variables, and is valid for multi-axial, non-monotonic and non-proportional loading histories. The theory is applied in detail to a well-known constitutive model, and verified against finite element simulations of representative volume elements or unit cells, for a number of composite materials.     

Multiple cracks in thermoelectroelastic bimaterials

The formulation for thermal stress and electric displacement in an infinite thermopiezoelectric plate with an interface and multiple cracks is presented. Using Green's function approach and the principle of superposition, a system of singular integral equations for the unknown temperature discontinuity defined on each crack face is developed and solved numerically. The formulation can then be used to calculate some fracture parameters such as the stress–electric displacement and strain energy density factor. The direction of crack growth for many cracks in thermopiezoelectric bimaterials is predicted by way of the strain energy density theory. Numerical results for stress–electric displacement factors and crack growth direction at a particular crack tip in two crack system of bimaterials are presented to illustrate the application of the proposed formulation.     

Micromechanics of Short-Term Thermal Damageability

The structural theory of microdamageability of a homogeneous material is generalized to the case of a thermal action. The theory is based on the stochastic thermoelastic equations of a medium with micropores, hollow or filled with particles of a damaged material. This medium models a material with dispersed microdamages. The Schleicher–Nadai fracture criterion is used as the condition of origin of a micropore in a microvolume of an undestroyed material. It is assumed that the particles of the damaged material in the micropores do not resist shear and triaxial tension and behave as the undamaged material under triaxial compression. The porosity balance equation is corrected for the thermal component and together with the relations between macrostresses, macrostrains, and temperature forms a closed system describing the concurrent action of deformation and microdamage. Nonlinear stress–strain diagrams and dependences of microdamage on macrostrain and temperature are constructed     

First-order gradient damage theory

Taking the strain tensor, the scalar damage variable, and the damage gradient as the state variables of the Helmholtz free energy, the general expressions of the firstorder gradient damage constitutive equations are derived directly from the basic law of irreversible thermodynamics with the constitutive functional expansion method at the natural state. When the damage variable is equal to zero, the expressions can be simplified to the linear elastic constitutive equations. When the damage gradient vanishes, the expressions can be simplified to the classical damage constitutive equations based on the strain equivalence hypothesis. A one-dimensional problem is presented to indicate that the damage field changes from the non-periodic solutions to the spatial periodic-like solutions with stress increment. The peak value region develops a localization band. The onset mechanism of strain localization is proposed. Damage localization emerges after damage occurs for a short time. The width of the localization band is proportional to the internal characteristic length.     

热冲击下物性与温度相关性对热弹性响应的渐进分析

计及材料物性与温度的相关性,基于Green-Naghdi能量无耗散广义热弹性理论(G-N II理论),对热冲击下具有变物性特征材料的热弹性响应进行了求解分析。借助Laplace正、反变换技术以及Krichhoff变换,在热物性参数随真实温度呈线性规律的前提下,推导了半无限大体受热冲击作用时热弹性响应的解析表达式,通过求解分析,得到了热冲击下热波、热弹性波的传播规律,位移场、温度场以及应力场的分布情况,以及物性随温度相关性对热弹性响应的影响效果。结果表明：当考虑材料物性随温度的变化时,热波、热弹性波的传播以及各物理场的分布均受到不同程度的影响,且物性随温度相关性对热弹性响应的作用效果将受到材料热-力耦合特性的影响。     

On tensorial forms of thermodynamic potentials in mixtures theory

Four thermodynamic tensorial quantities, equivalents of the thermodynamic potentials, namely the chemical potential tensor, the tensor of enthalpy, the tensor of free energy and the tensor of internal energy are presented. The last three of them are proposed in this paper and connections between all of these tensors are derived. The tensorial forms of the thermodynamic potentials are expressed in terms of four possible pairs of independent state variables. A set of sixteen alternative expressions, four for each tensorial form of the thermodynamic potential are derived and their importance discussed.     

On the Gibbs conditions of stable equilibrium,convexity and the second-order variations of thermodynamic potentials in nonlinear thermoelasticity

The Gibbs conditions of stable thermodynamic equilibrium are formulated for nonlinear thermoelastic materials, based on the constrained minimization of four fundamental thermodynamic potentials. Sufficient conditions for incremental stability are stated in each case. The previously unexplored connections between the second-order variations of thermodynamic potentials are used to establish the convexity or concavity properties of all thermodynamic potentials in relation to each other, and to derive the relationships between the specific heats at constant stress and deformation, and between the isentropic and isothermal elastic moduli and compliances. The comparison with the derivation based on the classical thermodynamic approach is also given.     

On thermoelastic dielectrics

Constitutive equations for a linear thermoelastic dielectric are derived from the energy balance equation assuming dependence of the stored energy function on the strain tensor, the polarization vector, the polarization gradient tensor and entropy. A method is indicated for constructing a hierarchy of constitutive equations for materials with arbitrary symmetry by introducing various thermodynamic potentials. Maxwell's relations are constructed for the thermodynamic potential W^L. The entropy inequality is used to obtain stability conditions for an elastic dielectric in equilibrium under prescribed boundary constraints. Frequencies are explicitly determined for a plane wave propagating along the x₁-axis in an infinite centro-symmetric isotropic thermoelastic dielectric.     

A Note on the Theory of Short-Term Microdamageability of Granular Composites under Thermal Actions

The theory of microdamageability of granular composites is stated with allowance made for the thermal effect. Microdamages in the components are modeled by pores, hollow or filled with particles of the destroyed material that resist compression. The fracture criterion is assumed to have the Schleicher–Nadai form, which takes into account the difference between the tensile and compressive ultimate loads, with the ultimate strength being a random function with a power or Weibull distribution. The stress–strain state and effective properties of the material are determined from the stochastic thermoelastic equations for granular composites with porous components. The equations of deformation and microdamage are closed by the equation of porosity balance corrected for the thermal effect. Nonlinear diagrams are plotted for the concurrent processes of deformation of a granular material and microdamage of the matrix as functions of macrostrains and temperature. The influence of the physical and geometrical parameters on the processes is analyzed.     

饱和土中的任意形状孔洞对弹性波的散射

根据Biot波动理论建立了求解饱和土中任意形状孔洞对弹性波散射问题的复变函数方法.首先通过引入位移势函数把稳态条件下的Biot波动方程解耦为势函数所满足的Helmholtz方程.利用分离变量方法即得到Helmholtz方程完备的通解.根据所得位移势函数的通解,可得骨架位移、流体相对骨架的位移、应力和孔压的表达式.通过保角变换方法,把物理平面上的孔洞映射到像平面上单位圆.利用土骨架和流体的边界条件,即可确定波函数展开式中的未知系数.给出了一些数值结果.     

Constraints, Constitutive Limits, and Instability in Finite Thermoelasticity

This paper examines all the possible types of thermomechanical constraints in finite-deformational elasticity. By a thermomechanical constraint we mean a functional relationship between a mechanical variable, either the deformation gradient or the stress, and a thermal variable, temperature, entropy or one of the energy potentials; internal energy, Helmholtz free energy, Gibbs free energy or enthalpy. It is shown that for the temperature-deformation, entropy-stress, enthalpy-deformation, and Helmholtz free energy-stress constraints equilibrium states are unstable, in the sense that certain perturbations of the equilibrium state grow exponentially. By considering the constrained materials as constitutive limits of unconstrained materials, it is shown that the instability is associated with the violation of the Legendre–Hadamard condition on the internal energy. The entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints do not exhibit this instability. It is proposed that stability of the rest state (or equivalently convexity of internal energy) is a necessary requirement for a model to be physically valid, and hence entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints are physical, whereas temperature-deformation constraints (including the customary notion of thermal expansion that density is a function of temperature only), entropy-stress constraints, enthalpy-deformation constraints, and Helmholtz free energy-stress constraints are not.     

Elastodynamic analysis at an interface of viscous fluid/thermoelastic micropolar honeycomb medium due to inclined load

In this paper, the effect of angle inclination at the interface of a viscous fluid and thermoelastic micropolar honeycomb solid due to inclined load is investigated. The inclined load is assumed to be a linear combination of normal load and tangential load. Laplace transform with respect to time variable and Fourier transform with respect to space variable are applied to solve the problem. Expressions of stresses, temperature distribution, and pressures in the transformed domain are obtained by introducing potential functions. The numerical inversion technique is used to obtain the solution in the physical domain. The frequency domain expressions for steady state are also obtained with appropriate change of variables. Graphic representations due to the response of different sources and changes of angle inclination are shown. Some particular cases are also discussed.     

EFFECT OF FRACTIONAL ORDER PARAMETER ON THERMOELASTIC BEHAVIORS OF ELASTIC MEDIUM WITH VARIABLE PROPERTIES

This paper is concerned with the thermoelastic behaviors of an elastic medium with variable thermal material properties. The problem is in the context of fractional order heat conduction. The governing equations with variable thermal properties were established by means of the fractional order calculus. The problem of a half-space formed of an elastic medium with variable thermal material properties was solved, and asymptotic solutions induced by a sudden temperature rise on the boundary were obtained by applying an asymptotic approach. The propagations of thermoelastic wave and thermal wave, as well as the distributions of displacement,temperature and stresses were obtained and plotted. Variations in the distributions with different values of fractional order parameter were discussed. The results were compared with those obtained from the case of constant material properties to evaluate the effects of variable material properties on thermoelastic behaviors.     

Representation of the glass-transition in mechanical and thermal properties of glass-forming materials: A three-dimensional theory based on thermodynamics with internal state variables

In the vicinity of the glass transition, glass-forming materials exhibit pronounced frequency-dependent changes in the mechanical material properties, the thermal expansion behaviour and the specific heat. The frequency dependence becomes apparent under harmonic stress, strain or temperature excitations. The Prigogine-Defay ratio is a characteristic number which connects the changes in magnitude of these quantities at the glass transition. In order to represent the thermoviscoelastic properties of glass-forming materials in continuum mechanics, a three-dimensional approach which is based on the Gibbs free energy as thermodynamic potential is developed in this article. The Gibbs free energy depends on the stress tensor, the temperature and a set of internal variables which is introduced to take history-dependent phenomena into account. In the vicinity of an equilibrium reference state, the specific Gibbs free energy is approximated up to second order terms. Evaluating the Clausius-Duhem inequality, the constitutive relations for the strain tensor, the entropy and the internal variables are derived. In comparison with other approaches, the entropy, the strain tensor and the internal variables are functionals not only of the stress tensor but also of the temperature. Applying harmonic temperature- or stress-controlled excitations, complex frequency-dependent relations for the specific heat under constant stress, for the thermal expansion coefficients as well as for the dynamic mechanical compliance are obtained. The frequency-dependence of these quantities depicts the experimentally observed behaviour of glass-forming materials as published in literature. Under the assumption of isotropic material behaviour, it is shown that the developed theory is compatible with the Prigogine-Defay inequality for arbitrary values of the material parameters.     

Mechanical behavior of amorphous polymers in shear

Based on the non-equilibrium thermodynamic theory, a new thermo-viscoelastic constitutive model for an incompressible material is proposed. This model can be considered as a kind of generalization of the non-Gaussian network theory in rubber elasticity to include the viscous and the thermal effects. A set of second rank tensorial internal variables was introduced, and in order to adequately describe the evolution of these internal variables, a new expression of the Helmholtz free energy was suggested. The mechanical behavior of the thermo-viscoelastic material under simple shear deformation was studied, and the “ viscous dissipation induced“ anisotropy due to the change of orientation distribution of molecular chains was examined. Influences of strain rate and thermal softening produced by the viscous dissipation on the shear stress were also discussed. Finally, the model predictions were compared with the experimental results performed by G‘ Sell et al. , thus the validity of the proposed model is verified.