Thermoconductivity and thermoelasticity of thin isotropic shells heated by a moving heat pulse

The paper presents a solution to the problem of thermal conduction and thermoelasticity for a thin shallow spherical shell heated by a concentrated or local impulsive heat source moving over the shell surface. It is assumed that temperature is linearly distributed across the shell thickness and that the shell, on its sides, exchanges heat with the environment in accordance with Newton’s law of cooling. The Fourier and Laplace transforms are used to find an analytic solution. The dependence of the temperature field and stress/strain components on the type of heating and the form of heat source is studied __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 11, pp. 85–92, November 2006.     

On the theory of elastic shells made from a material with voids

In this paper we present a theory for porous elastic shells using the model of Cosserat surfaces. We employ the Nunziato–Cowin theory of elastic materials with voids and introduce two scalar fields to describe the porosity of the shell: one field characterizes the volume fraction variations along the middle surface, while the other accounts for the changes in volume fraction along the shell thickness. Starting from the basic principles, we first deduce the equations of the nonlinear theory of Cosserat shells with voids. Then, in the context of the linear theory, we prove the uniqueness of solution for the boundary initial value problem. In the case of an isotropic and homogeneous material, we determine the constitutive coefficients for Cosserat shells, by comparison with the results derived from the three-dimensional theory of elastic media with voids. To this aim, we solve two elastostatic problems concerning rectangular plates with voids: the pure bending problem and the extensional deformation under hydrostatic pressure.     

Analytical solution to the unsteady one-dimensional conduction problem with two time-varying boundary conditions: Duhamel’s theorem and separation of variables

An analytical resolution of the time-dependent one-dimensional heat conduction problem with time-dependent boundary conditions using the method of separation of variables and Duhamel’s theorem is presented. The two boundary conditions used are a time-dependent heat flux at one end and a varying temperature at the other end of the one-dimensional domain. It is put forth because the author found that the prescribed resolution method using separation of variables and Duhamel’s theorem presented in heat conduction textbooks is not directly applicable to problems with more than one time-dependent boundary condition. The analytical method presented in this paper makes use of one of the property of the heat conduction equation: the apparent linearity of the solutions. For that reason, in order to solve a problem with two time-dependent boundary conditions, the author first separates the initial problem into two independent but complementary problems, each with only one time-dependent boundary condition. Doing that, both simpler problems can be solved independently using a prescribed method that is known to work and the final solution can be obtained by joining the two independent solutions from the simpler separated problems. Every step of the resolution method is presented in this paper, along with a numerical validation of the final solution of three test case problems.     

机械荷载与热荷载共同作用下变厚度梁的热弹性力学解

本文研究了两端简支变厚度梁受机械荷载与热荷载共同作用下的热弹性力学解.温度场用调和级数展开,通过求解热传导方程可首先确定温度沿梁厚度方向的非线性分布情况.从二维热弹性力学理论的基本方程出发,导出满足控制微分方程和两端简支边界条件的位移函数的一般解,对上下表面的边界方程作傅立叶正弦级数展开确定待定系数,数值结果与商业有限元软件ANSYS进行了比较,显示出很高的精度.本文方法可直接应用于对应力和位移分析要求较高的工程问题,如航空航天和微型机械的设计.     

基于热-磁-弹性耦合的二维导电薄柱壳数值分析

基于几何方程、物理方程、运动方程和电动力学方程,建立了二维导电薄柱壳热磁弹性基本方程．考虑到Joule热效应,引入热平衡方程及广义Ohm定律,得出导电薄柱壳的温度场．利用变量代换方法,整理成具有10个基本未知量的标准型方程组．采用差分法和准线性化方法,给出准线性微分方程组．对于二维导电薄柱壳,导出了Lorentz力表达式、温度场积分特征值．通过实例计算,得到了二维导电薄柱壳应力、位移、温度随外加电磁参量的变化规律．研究结果可为二维导电薄壳热磁弹性问题研究提供理论参考．     

Transient hyperbolic heat conduction in thick-walled FGM cylinders and spheres with exponentially-varying properties

This paper focuses on non-Fourier hyperbolic heat conduction analysis for heterogeneous hollow cylinders and spheres made of functionally graded material (FGM). All the material properties vary exponentially across the thickness, except for the thermal relaxation parameter which is taken to be constant. The cylinder and sphere are considered to be cylindrically and spherically symmetric, respectively, leading to one-dimensional heat conduction problems. The problems are solved analytically in the Laplace domain, and the results obtained are transformed to the real-time space using the modified Durbin’s numerical inversion method. The transient responses of temperature and heat flux are investigated for different inhomogeneity parameters and relative temperature change values. The comparisons of temperature distribution and heat flux between various time and material properties are presented in the form of graphs.     

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method. The solution accuracy is further improved by implementing an adaptive meshing technique to generate finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems are: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to a moving heat source. The English text was polished by Yunming Chen.     

An exact solution for the steady-state thermoelastic response of functionally graded orthotropic cylindrical shells

We analyze the steady-state response of a functionally graded thick cylindrical shell subjected to thermal and mechanical loads. The functionally graded shell is simply supported at the edges and it is assumed to have an arbitrary variation of material properties in the radial direction. The three-dimensional steady-state heat conduction and thermoelasticity equations, simplified to the case of generalized plane strain deformations in the axial direction, are solved analytically. Suitable temperature and displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the thermoelastic equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which are then solved by the power series method. In the present formulation, the cylindrical shell is assumed to be made of an orthotropic material, although the analytical solution is also valid for isotropic materials. Results are presented for two-constituent isotropic and fiber-reinforced functionally graded shells that have a smooth variation of material volume fractions, and/or in-plane fiber orientations, through the radial direction. The cylindrical shells are also analyzed using the Flügge and the Donnell shell theories. Displacements and stresses from the shell theories are compared with the three-dimensional exact solution to delineate the effects of transverse shear deformation, shell thickness and angular span.     

A boundary identification method for an inverse heat conduction problem with an application in ironmaking

A boundary identification problem in inverse heat conduction is studied, based on data from internal measurement of temperature and heat flux. Formulated as a sideways heat conduction equation, a spatial continuation technique is applied to extend the solution to a known boundary condition at the desired boundary position. Recording the positions traversed in the continuation for each time instant yields the boundary position trajectory and hence the solution of the identification problem. A prospective application of the method can be found in the ironmaking blast furnace, where it is desired to monitor the thickness of the accreted refractory wall based on measurement of its internal state. Simulations featuring noisy measurement data demonstrate the feasibility of the identification method for blast furnace wall thickness estimation.     

基于遗传算法的混凝土三维非稳态温度场反分析

大体积混凝土结构施工期间的合理温度控制问题非常重要，而精确进行温度计算所需的一些材料参数往往不易直接测得，需要根据一些易测得的量进行反求。本文基于三维瞬态温度场有限元求解理论与反问题理论，建立了混凝土三维瞬态温度场反问题求解数值模型。运用遗传算法寻求非线性反演问题全局最优解，只需要若干点温度实测值便可实现混凝土多个热学参数如绝热温升、导温、导热系数及热交换系鼓等的同时反演，算例对本文反演方法的反演精度及数值稳定性给出了满意的证明。     

Pulsed laser heating – Fourier and electron kinetic theory approaches

Laser surface pulse heating of engineering metals is in demand in the metal industry and investigation into laser pulse heating becomes fruitful in this regard. Application of Fourier theory to heat conduction due to high power laser irradiation may give closed form solution to the problem. On the other hand, the heat flux through a given plane depends on the electron energy distribution through the material and at the scale of distance required to examine the problem, the material can no longer be considered as being homogeneous continuum, therefore, errors may occur when considering the Fourier theory in laser heating process. The problem requires to be examined in the quantum field. The present study examines the pulse laser heating process when considering both Fourier conduction and electron-kinetic theory approaches. Analytical solution to Fourier conduction equation is obtained for intensity exponential pulses while numerical scheme is introduced to solve the heat transfer equation resulted from kinetic theory approach. It is found that both Fourier and electron kinetic theory approaches result in similar temperature profiles for the pulses having the same energy content. In the case of electron kinetic theory approach the rise time for surface temperature to reach the melting point is shorter than that obtained from the analytical solution. Received on 23 February 1998     

Homogenization of the equations of the Cosserat theory of elasticity of inhomogeneous bodies

The paper deals with the homogenization of a boundary value problem for an inhomogeneous body with Cosserat properties, which is referred to as the original problem. The homogenization process is understood as a method for representing the solution of the original problem in terms of the solution of precisely the same problem for a body with homogeneous properties. The problem for a body with homogeneous properties is called the accompanying problem, and the body itself, the accompanying homogeneous body. As a rule, a constructive homogenization procedure includes the following three stages: at the first stage, the properties of the inhomogeneous body are used to find the properties of the accompanying homogeneous body (efficient properties); at the second stage, the boundary value problem is solved for the accompanying body; at the third stage, the solution of the accompanying problem is used to find the solution of the original problem. This approach was implemented in mechanics of composite materials constructed of numerous representative elements. A significant contribution to the development of mechanics of composites is due to Rabotnov [1–3] and his students. Recently, the homogenization method has been widely used to solve problems for composites of regular structure by expanding the solution of the original problem in a power series in a small geometric parameter equal to the ratio of the characteristic dimension of the periodicity cell to the characteristic dimension of the entire body. The papers by Bakhvalov [4–6] and Pobedrya [7] were the first in the field. At present, there are numerous monographs partially or completely dealing with the method of a small geometric parameter [8–14]. Isolated problems for inhomogeneous bodies with nonperiodic dependence of their properties on the coordinates were considered by many authors. Most of such papers published before 1973 are collected in two vast bibliographic indices [15, 16]. General methods were considered, and many specific problems of the theory of elasticity of continuously inhomogeneous bodies were solved in Lomakin’s papers and his monograph [17]. The theory of torsion of inhomogeneous anisotropic rods was considered in [18]. In 1991, in his Doctoral dissertation, one of the authors of this paper proposed a version of the homogenization method based on an integral formula representing the solution of the original static problem of inhomogeneous elasticity via the solution of the accompanying problem [19, 20]. An integral formula for the dynamic problem of elasticity was published somewhat later [21]. This integral formula was used to develop a constructive method for the homogenization of the dynamic problem of inhomogeneous elasticity, which can be used in the case of both periodic and nonperiodic inhomogeneity of the properties [22]. The integral formula in the case of the Cosserat theory of elasticity was published in [23]. The present paper briefly presents constructive methods for homogenizing the problems of the Cosserat theory of elasticity based on the integral formula.     

Transient thermal cracking associated with non-classical heat conduction in cylindrical coordinate system

This paper studies the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction.The heat conduction equation is solved by separation of variable technique.Closed form solution for the temperature field and the associated thermal stress are established.The critical parameter governing the level of the transient thermal stress is identified.Exact expression for the transient stress intensity factor is obtained for a crack in the cylinder.The difference between the non-classical solutions and the classical solution are discussed.It is found that the traditional classical heat conduction considerably underestimates the transient thermal stress and thermal stress intensity factor.     

Approximate analysis for transient heat conduction with radiation in a slab

In an attempt to minimize the numerical computations associated with the solution of transient heat conduction with radiation in a slab, a perturbation type of analysis is being applied to the temperature field and radiation heat flux simultaneously. The resulting partial differential equations for the perturbation functions for the temperature are solved in explicit forms by use of the energy integral methods, while the radiation heat flux is determined by an appropriate scheme of approximating the temperature distribution in the slab. Included in the analysis are the effects of the parameters: the optical thickness, the ratio of conduction transport to radiation and the wall emissivity. It is found that, in a wide range of these governing parameters, the results compare very favorably with those obtained by the numerical solution of the formulated integro-differential equation. With the present analysis, the temperature, conduction and radiation heat fluxes can be predicted without resorting to lengthy numerical analysis.     

AN ANALYSIS OF TRANSIENT HEAT CONDUCTION FOR FUNCTIONALLY GRADED SPHERES BASED ON STATE SPACE THEORY1)

基于状态空间理论研究功能梯度圆球的球对称瞬态热传导问题。根据热传导方程和热流密度的定义,取温度场和热流密度为系统的状态向量,通过将圆球分层和在时域内应用差分格式对控制方程进行离散,建立了系统的状态方程,给出了功能梯度圆球瞬态热传导问题的半解析解。算例分析表明：本文解不但结果正确、计算效率高,而且适用于材料参数沿径向任意梯度变化的圆球瞬态热传导分析。     

Thermomechanical formulation of ductile damage coupled to nonlinear isotropic hardening and multiplicative viscoplasticity

In this paper, we present a thermomechanical framework which makes use of the internal variable theory of thermodynamics for damage-coupled finite viscoplasticity with nonlinear isotropic hardening. Damage evolution, being an irreversible process, generates heat. In addition to its direct effect on material's strength and stiffness, it causes deterioration of the heat conduction. The formulation, following the footsteps of Simó and Miehe (1992), introduces inelastic entropy as an additional state variable. Given a temperature dependent damage dissipation potential, we show that the evolution of inelastic entropy assumes a split form relating to plastic and damage parts, respectively. The solution of the thermomechanical problem is based on the so-called isothermal split. This allows the use of the model in 2D and 3D example problems involving geometrical imperfection triggered necking in an axisymmetric bar and thermally triggered necking of a 3D rectangular bar.     

功能梯度板条的热弹性力学解

本文利用推广后的Mian 和Spencer 功能梯度板理论,研究了功能梯度板条在非均布温度场作用下的热弹性问题．采用该理论中的位移展开公式,在板厚度方向上考虑热传导引起的稳态温度场,材料常数沿板厚方向可以任意连续变化,从而得到了基于弹性理论的功能梯度板条在温度场作用下的解析解．通过数值算例分析,验证了本文理论的正确性并讨论了边界条件和梯度变化程度对功能梯度板条热弹性响应的影响．     

Thermal stresses factors of orthotropic graded layers with a finite line bond

In the paper, the problem of a finite line bond between two orthotropic functionally graded strips under thermal loading is considered. Considering some new boundary conditions, it is assumed that the temperature drop across the finite line bond is the result of the thermal conductivity index controlling heat conduction through the bond region. Using Fourier transforms technique, the therm-elastic mixed boundary value problems are reduced to a system of singular integral equations which can be solved approximately by applying the Chebyshev polynomials. The numerical results for the temperature and displacement field as well as the thermal stress intensity factors (TSIFs) are presented. The influence of the thickness of the layers and the thermo-elastic nonhomogeneity parameters on the temperature distribution and the TSIFs is discussed in detail. These results can be expected to be used for the purpose of gaining better understanding of the thermo-mechanical behavior of layered structures.     

Transient heat conduction in functionally graded thick hollow cylinders by analytical method

In this article, transient heat conduction in a cylindrical shell of functionally graded material is studied by using analytical method. The shell is assumed to be in axisymmetry conditions. The material properties are considered to be nonlinear with a power law distribution through the thickness. The temperature distribution is derived analytically by using the Bessel functions. To verify the proposed method the obtained numerical results are compared with the published results. The comparisons of temperature distribution between various time and material properties are presented.     

THE SCHRÖDINGER EQUATION OF THIN SHELL THEORIES

This work is the continuation of the discussions of [50] and [51]. In this paper: (A) The Love-Kirchhoff equation of small deflection problem for elastic thin shell with constant curvature are classified as the same several solutions of Schrodinger equation, and we show clearly that its form in axisymmetric problem;(B) For example for the small deflection problem, we extract me general solution of the vibration problem of thin spherical shell with equal thickness by the force in central surface and axisymmetric external field, that this is distinct from ref. [50] in variable. Today the variable is a space-place, and is not time;(C) The von Kármán-Vlasov equation of large deflection problem for shallow shell are classified as the solutions of AKNS equations and in it the one-dimensional problem is classified as the solution of simple Schrodinger equation for eigenvalues problem, and we transform the large deflection of shallow shell from nonlinear problem into soluble linear problem.