与温度相关的非均质圆筒的非线性耦合广义热弹性分析

基于能量守恒方程和L–S广义热弹性理论,借助状态空间技术和Newmark法求解了材料性质沿径向任意梯度分布同时又与温度相关的非均质圆筒非线性耦合广义热弹性问题。通过对材料性质与温度无关和相关功能梯度圆筒的算例分析,给出了在线性和非线性耦合下圆筒内温度和应力沿径向和随时间的变化关系,验证了本文解的正确性和有效性。数值结果表明,考虑材料性质是否与温度相关,能量守恒方程中耦合项是线性还是非线性,得到的温度与应力均存在不同程度的差异。本文解可方便地应用于不同边界条件和初始条件下圆筒的广义热弹性分析。

     

A state space formalism for anisotropic elasticity. Part I: Rectilinear anisotropy

A state space formalism for anisotropic elasticity including the thermal effect is developed. A salient feature of the formalism is that it bridges the compliance-based and stiffness-based formalisms in a natural way. The displacement and stress components and the thermoelastic constants of a general anisotropic elastic material appear explicitly in the formulation, yet it is simple and clear. This is achieved by using the matrix notation to express the basic equations and grouping the stress in such a way that it enables us to cast neatly the three-dimensional equations of anisotropic elasticity into a compact state equation and an output equation. The homogeneous solution to the state equation for the generalized plane problem leads naturally to the eigen relation and the sextic equation of Stroh. Extension, twisting, bending, temperature change and body forces are accounted for through the particular solution. Based on the formalism the general solution for generalized plane strain and generalized torsion of an anisotropic elastic body are determined in an elegant manner.     

DIFFERENTIAL-ALGEBRAIC APPROACH TO COUPLED PROBLEMS OF DYNAMIC THERMOELASTICITY

An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations (DAEs) where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula (BDF) using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.     

Ideal nonsymmetric 4D-medium as a model of invertible dynamic thermoelasticity

The space-time continuum (4D-medium) is considered, and a generalized model of reversible dynamic thermoelasticity is constructed as a theory of elasticity of an ideal (defect-free) nonsymmetric 4D-medium that is transversally-isotropic with respect to the time coordinate. The definitions of stresses and strains for the space-time continuum are introduced. The constitutive equations of the medium model relating the components of nonsymmetric stress and distortion 4D-tensors are stated. Physical interpretations of all tensor components of the thermomechanical properties are given. The Lagrangian of the generalized model of coupled dynamic thermoelasticity is presented, and the Euler equations are analyzed. It is shown that the three Euler equations are generalized equations of motion of the dynamic classical thermoelasticity, and the last, fourth, equation is a generalized heat equation which allows one to predict the wave properties of heat. An energy-consistent version of thermoelasticity is constructed where the Duhamel-Neumann and Maxwell-Cattaneo laws (a nonclassical generalization of the Fourier law for the heat flow) are direct consequences of the constitutive equations.     

Estimates of thermoelastic characteristics of composites reinforced by short anisotropic fibers

We construct a mathematical model describing thermomechanical interaction between composite structure elements (isotropic particles of the matrix and anisotropic short fibers) and the macroscopically isotropic elastic medium with desired thermoelastic characteristics. At the first stage of this model, the self-consistency method is used to obtain relations determining the elasticity moduli of the composite, and at the second stage, the model permits determining its linear thermal expansion coefficient. The dual variational statement of the linear thermoelasticity problem in an inhomogeneous solid permits obtaining two-sided estimates for the bulk elasticity modulus, shear modulus, and linear thermal expansion coefficient of the composite under study. The calculated dependencies presented in the paper permit predicting the thermoelastic characteristics of a composite reinforced by anisotropic short fibers (including those in the form of nanostructure elements).     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅴ)--POLAR THERMOMECHANICAL CONTINUA

The purpose is to reestablish rather complete basic balance equations and boundary conditions for polar thermomechanical continua based on the restudy of the traditional theories of micropolar thermoelasticity and thermopiezoelectricity . The equations of motion and the local balance equation of energy rate for micropolar thermoelasticity are derived from the rather complete principle of virtual power. The equations of motion, the balance equation of entropy and all boundary conditions are derived from the rather complete Hamilton principle . The new balance equations of momentum and energy rate which are essentially different from the existing results are presented. The corresponding results of micromorphic thermoelasticity and couple stress elastodynamics may be naturally obtained by the transition and the reduction from the micropolar case , respectively . Finally , the results of micropolar thermopiezoelectricity are directly given .     

Spatial Stability for the Quasi-static Problem of Thermoelasticity

This paper is concerned with the quasi-static problem of thermoelasticity. The classical system of equations of thermoelasticity is a coupling of an elliptic equation with a parabolic equation. It poses some new mathematical difficulties. Here we study the exponential spatial decay of solutions. An upper bound for the amplitude in terms of the boundary and initial conditions is obtained. The extension of the spatial stability results to thermoelasticity of type III is also treated. Dedicated to C.O. Horgan on the occasion of his 60th birthdayMathematics Subject Classifications (2000) 74F05, 74G50.     

Impossibility of localization in linear thermoelasticity with voids

In this note we prove the impossibility of the localization in time of the solutions of the linear thermoelasticity with voids. This means that the only solution for this problem that vanishes after a finite time is the null solution. From a thermomechanical point of view, this result says that the combination of the thermal and porous dissipation in the linear theory is not sufficiently strong to guarantee that the thermomechanical deformations will vanish after a finite time. The main idea to prove this result is to show the uniqueness of solutions for the backward in time problem.     

Multiple reciprocity boundary element analysis of two-dimensional anisotropic thermoelasticity involving an internal arbitrary non-uniform volume heat source

In the direct boundary element method (BEM) formulation of anisotropic thermoelasticity, thermal loads manifest themselves as additional volume integral terms in the boundary integral equation (BIE). Conventionally, this requires internal cell discretisation throughout the whole domain. In this paper, the multiple reciprocity method in BEM analysis is employed to treat the general 2D thermoelasticity problem when the thermal loading is due to an internal non-uniform volume heat source. By successively performing the “volume-to-surface” integral transformation, the general formulation of the associated BIE for the problem is derived. The successful implementation of such a scheme is illustrated by three numerical examples.     

Convergence of Anisotropically Decaying Solutions of a Supercritical Semilinear Heat Equation

We consider the Cauchy problem for a semilinear heat equation with a supercritical power nonlinearity. It is known that the asymptotic behavior of solutions in time is determined by the decay rate of their initial values in space. In particular, if an initial value decays like a radial steady state, then the corresponding solution converges to that steady state. In this paper we consider solutions whose initial values decay in an anisotropic way. We show that each such solution converges to a steady state which is explicitly determined by an average formula. For a proof, we first consider the linearized equation around a singular steady state, and find a self-similar solution with a specific asymptotic behavior. Then we construct suitable comparison functions by using the self-similar solution, and apply our previous results on global stability and quasi-convergence of solutions.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅴ)——POLAR THERMOMECHANICAL CONTINUA

The purpose is to reestablish rather complete basic balance equations and boundary conditions for polar thermomechanical continua based on the restudy of the traditional theories of micropolar thermoelasticity and thermopiezoelectricity. The equations of motion and the local balance equation of energy rate for micropolar thermoelasticity are derived from the rather complete principle of virtual power. The equations of motion, the balance equation of entropy and all boundary conditions are derived from the rather complete Hamilton principle. The new balance equations of momentum and energy rate which are essentially different from the existing results are presented. The corresponding results of micromorphic thermoelasticity and couple stress elastodynamics may be naturally obtained by the transition and the reduction from the micropolar case, respectively. Finally, the results of micropolar thermopiezoelectricity are directly given. Contributed by DAI Tian-min, Original Member of Editorial Committee, AMM Foundation items: the National Natural Science Foundation of China (10072024); the Research Foundation of Liaoning Education Committee (990111001) Biography: DAI Tian-min (1931-)     

A boundary element method for anisotropic coupled thermoelasticity

Summary A boundary element formulation is presented for the solution of the equations of fully coupled thermoelasticity for materials of arbitrary degree of anisotropy. By employing the fundamental solutions of anisotropic elastostatics and stationary heat conduction, a system of equations with time-independent matrices is obtained. Since the fundamental solutions are uncoupled and time-independent, a domain integral remains in the representation formula which contains the time-dependence as well as the thermoelastic coupling. This domain integral is transformed to the boundary by means of the dual reciprocity method. By taking this approach, the use of dynamic fundamental solutions is avoided, which enables an efficient calculation of system matrices. In addition, the solution of transient processes as well as, free and forced vibration analysis becomes straightforward and can be carried out with standard time-stepping schemes and eigensystem solvers. Another important advantage of the present formulation is its versatility, since it includes a number of simplified thermoelastic theories, viz. the theory of thermal stresses, coupled and uncoupled quasi-static thermoelasticity, and stationary thermoelasticity. The accuracy of the new thermoelastic boundary element method is demonstrated by a number of example problems. Support by the Deutsche Forschungsgemeinschaft (DFG) of the Graduate Collegium Modelling and discretization methods for continua and fluids (GKKS) at the University of Stuttgart is gratefully acknowledged.     

Evolutionary model of finite-strain thermoelasticity

The equation of state of finite-strain thermoelasticity is obtained using a formalized approach to constructing constitutive relations for complex media under the assumption of closeness of intermediate and current configurations. A variational formulation of the coupled thermoelastic problem is proposed. The constitutive equation, the heat-conduction equation, the relations for internal energy, free energy, and entropy, and the variational formulation of the coupled problem of finite-strain thermoelasticity are tested on the problem of uniaxial extension of a bar. The model adequately describes experimental data for elastomers, such as entropic elasticity, temperature inversion, and temperature variation during an adiabatic process. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 184–196, May–June, 2008.     

Closed solutions of boundary-value problems of coupled thermoelasticity

Coupled equations of thermoelasticity take into account the effect of nonuniform heating on the medium deformation and that of the dilatation rate on the temperature distribution. As a rule, the coupling coefficients are small and it is assumed, sometimes without proper justification, that the effect of the dilatation rate on the heat conduction process can be neglected. The aim of the present paper is to construct analytical solutions of some model boundary-value problems for a thermoelastic bounded body and to determine the body characteristic dimensions and the medium thermomechanical moduli forwhich it is necessary to take into account that the temperature and displacement fields are coupled. We consider some models constructed on the basis of the Fourier heat conduction law and the generalized Cattaneo-Jeffreys law in which the heat flux inertia is taken into account. The solution is constructed as an expansion in a biorthogonal system of eigenfunctions of the nonself-adjoint operator pencil generated by the coupled equations of motion and heat conduction. For the model problem, we choose a special class of boundary conditions that allows us to exactly determine the pencil eigenvalues.     

The effect of anisotropic damage evolution on the behavior of ductile and brittle matrix composites

Anisotropic damage evolution laws for ductile and brittle materials have been coupled to a micromechanical model for the prediction of the behavior of composite materials. As a result, it is possible to investigate the effect of anisotropic progressive damage on the macroscopic (global) response and the local spatial field distributions of ductile and brittle matrix composites. Two types of thermoinelastic micromechanics analyses have been employed. In the first one, a one-way thermomechanical coupling in the constituents is considered according to which the thermal field affects the mechanical deformations. In the second one, a full thermomechanical coupling exists such that there is a mutual interaction between the mechanical and thermal fields via the energy equations of the constituents. Results are presented that illustrate the effect of anisotropic progressive damage in the ductile and brittle matrix phases on the composite’s behavior by comparisons with the corresponding isotropic damage law and/or by tracking the components of the damage tensor.     

On the plane problem of orthotropic quasi-static thermoelasticity

The plane displacement boundary value problem of quasi-static linear orthotropic thermoelasticity is discussed. The thermoelastic system on a bounded simply-connected domain is decoupled. The decoupled temperature equation is investigated by using an accurate estimate and the contractive mapping principle. Representation of solution of the field equation is obtained, and some solvability results are proved. The results are of both theoretical and numerical interest.     

Scattering of plane SH-waves on cylindrical canyon of arbitrary shape in anisotropic media

The complex function method is used to solve problems of scattering of plane SH-waves on cylindrical canyon topography of arbitrary shape in anisotropic media. This paper gives the complete function series and general expressions with boundary condition to approach the solution of steady state scattering of plane SH-waves on two-dimensional canyon topography in anisotropic media. The problem to be solved can be reduced to a solution of infinite algebraic equation series by using Hermite function and it's orthogonal conditions. The solution can be obtained directly by using computers. Finally, as an example, computational results of scattering of plane SH-waves on a semi-cylindrical canyon topography are presented. The projects sponsored by The Joint Seismological Science Foundation.     

Semi-analytical analysis of thermally induced damage in thin ceramic coatings

An efficient procedure to analyze damage evolution in brittle coatings under influence of thermal loads is suggested. The approach is based on a general computational scheme to determine damage evolution parameters, which incorporates an analytical solution of the appropriate interim boundary-value thermoelasticity problem. For thin inhomogeneous coatings, the simplification in the analysis is achieved by application of the mathematical model with generalized boundary conditions of thermomechanical conjugation of the substrate with environment via the coating. Efficiency of the suggested approach is illustrated by an example of damage evolution in the alumina coating on the titanium-alloy and tungsten substrates under uniform heating.     

Fundamental solutions of the thermoelastic equations for transversely isotropic plates

The problem of thermoelasticity for transversely isotropic plates acted upon by concentrated heat sources is solved. The {1, 2}-order equations of thermoelasticity that incorporate the transverse shear and normal stresses are used. A bending heat source with symmetric heat transfer is considered. The dependence of thermal stress components on the thermal and thermomechanical parameters of transversely isotropic plates is studied