分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

Theory of fractional order in electro-thermoelasticity

A new mathematical model of electro-thermoelasticity has been constructed in the context of a new consideration of heat conduction with fractional order. The state space approach developed in Ezzat (2008) is adopted for the solution of one-dimensional problem in the presence of heat sources. The Laplace transform technique is used. A numerical method is employed for the inversion of the Laplace transforms. According to the numerical results and its graphs, conclusion about the new theory has been constructed. Some comparisons have been shown in figures to estimate the effects of the fractional order parameter on all the studied fields.     

Synthesis of thermal effects in misaligned hydrodynamic journal bearings

The analysis presented herein deals with the evaluation of pressure and temperature fields which are generated in thin fluid films of varying thickness. The particular problem of a misaligned journal bearing has been studied by solving simultaneously the Reynolds and energy equations, which also include the effects of viscous dissipation and the variation of fluid viscosity with temperature. The method has been used to predict pressure and temperature fields as well as global performance parameters for a typical journal bearing operation.     

Localization of elastic deformation in strongly anisotropic, porous, linear materials with periodic microstructures: Exact solutions and dilute expansions

Exact solutions are derived for the problem of a two-dimensional, infinitely anisotropic, linear-elastic medium containing a periodic lattice of voids. The matrix material possesses either one infinitely soft, or one infinitely hard loading direction, which induces localized (singular) field configurations. The effective elastic moduli are computed as functions of the porosity in each case. Their dilute expansions feature half-integer powers of the porosity, which can be correlated to the localized field patterns. Statistical characterizations of the fields, such as their first moments and their histograms are provided, with particular emphasis on the singularities of the latter. The behavior of the system near the void close-packing fraction is also investigated. The results of this work shed light on corresponding results for strongly non-linear porous media, which have been obtained recently by means of the “second-order” homogenization method, and where the dilute estimates also exhibit fractional powers of the porosity.     

短脉冲激光加热分数阶导热及其热应力研究

短脉冲激光加热引起材料内部复杂的传热过程及热变形,现有的以Fourier定律或Cattaneo-Vernotte松弛方程结合弹性理论为框架建立起来热应力理论在刻画其热物理过程存在严重缺陷. 本文基于分数阶微积分理论, 以半空间为研究对象, 建立了分数阶Cattaneo热传导方程和相应的热应力方程, 给出了问题的初始条件和边界条件, 采用拉普拉斯变换方法, 给出了非高斯时间分布激光热源辐射下温度场和热应力场的解析解, 研究了短脉冲激光加热的温度场及热应力场的热物理行为. 数值计算中, 首先对理论解进行数值验证, 然后取分数阶变量$p=0.5$研究温度场和热应力场的变化特点及激光参数对温度和热应力的影响,最后数值计算分数阶参数对温度和热应力场的影响. 计算结果表明, 分数阶Cattaneo传热方程和热应力方程描述的温度和热应力任然具有波动特性,与经典的Fourier传热模型和标准的Cattaneo传热模型相比, 分数阶阶次越大, 热波波速越小, 热波波动性越明显; 反之, 则热波波速越大, 热扩散性越强.激光加热和冷却的速度越快, 温度上升和下降的速度越快, 压应力和拉应力交替变化越快, 温度变化幅值越小, 热应力幅值影响不明显.      

A unified treatment of the elastic elliptical inclusion under antiplane shear

Summary A generalized and unified treatment is presented for the antiplane problem of an elastic elliptical inclusion undergoing uniform eigenstrains and subjected to arbitrary loading in the surrounding matrix. The general solution to the problem is obtained through the use of conformal mapping technique and Laurent series expansion of the associated complex potentials. The resulting elastic fields are derived explicitly in both transformed and physical planes for the inclusion and the surrounding matrix. These relations are universal in the sense of being independent of any particular loading as well as the geometry of the matrix. The complete field solutions are provided for an elliptical inclusion under uniform loading at inifinity, and for a screw dislocation interacting with the elastic elliptical inclusion.     

Construction of two dimensional temperature field from first derivative fields

Schlieren and its variants such as differential interferometry and schlieren microscopy, have been used extensively for flow visualization where first derivative fields are captured. The derivative fields obtained from the schlieren like methods can be further processed to estimate the temperature field when the first derivative fields relate to those of temperature. Temperature construction from first derivative field is an ill-posed problem owing to the experimental noise and a few discrete points where measured temperatures may be available. A new approach has been proposed where the domain is discretized into a large number of triangular elements and least-squares based finite-element analysis is performed over the discretized domain. The domain and boundaries are identified manually based on prior knowledge. Temperature fields have been constructed for experimentally obtained first derivative fields from a Differential Interferometer (DI) for different cases. The performance of the new methodology is found to be good.     

An analysis of the effect of the solid layer for the modified chemical vapor deposition process

The heat transfer problem relative to the modified chemical vapor deposition process has been analyzed and the effects of solid layer thickness, torch speed and tube rotation are studied. The quasi-steady three-dimensional energy equations have been solved for the temperature fields in the gas and the solid layer with a Gaussian heat flux boundary condition on the outer surface. Of particular interest is the effect of the solid layer thickness and the torch speed on inner surface temperature, gas temperature and thermophoretic velocity. The large change of the axial temperature distribution of the surface occurs for different solid layer thicknesses or torch speeds. The presence of the solid layer and tube rotation reduce the effects of nonuniform torch heating in the circumferential direction and the resulting surface temperatures are very uniform in this direction.     

Interactions due to mechanical/thermal sources in a micropolar thermoelastic medium possessing cubic symmetry

The present problem is the deformation of micropolar thermoelastic solids with cubic symmetry under the influence of various sources acting on the plane surface. Analytic expressions for displacement components, microrotation, force stress, couple stress, and temperature distribution are obtained in the physical domain for Lord–Shulman (L–S) and Green–Lindsay (G–L) theories of thermoelasticity by applying integral transforms. A numerical inversion technique has been applied to obtain the solution in the physical domain. The numerical results are presented graphically for a particular model.     

Problem of generalized thermoelastic infinite medium with cylindrical cavity subjected to a ramp-type heating and loading

This paper presents the problem of thermoelastic interactions in an elastic infinite medium with cylindrical cavity at an elevated temperature field arising out of a ramp-type heating and loading bounding surface of the cavity, and the surface is assumed initially quiescent. The governing equations are taken in a unified system from which the field equations for coupled thermoelasticity as well as for generalized thermoelasticity can be easily obtained as particular cases. Due attention has been paid to the finite time of rise of temperature, stress, displacement, and strain. The problem has been solved analytically using a direct approach. The derived analytical expressions have been computed for a specific situation. Numerical results for the temperature distribution, thermal stress, displacement, and strain are represented graphically. A comparison is made with the results predicted by the three theories.     

Axisymmetric Dynamic Problem of Coupled Thermoviscoplasticity

An axisymmetric dynamic thermoviscoelastic problem is formulated with allowance for the coupling of mechanical and thermal fields. The behavior of the material is described by the Bodner–Partom model. A technique for numerical solution of the problem is developed. The laws governing the stress–strain state and the temperature field of a circular disk under forced flexural vibrations are studied     

饱和多孔黏弹地基热-水-力耦合动力响应分析

天然土体由于沉积条件和应力状态不同, 往往会表现出一定的流变性. 本文研究地基上表面受外载荷作用时, 渗透系数和孔隙率变化对饱和多孔黏弹性地基热-水-力耦合动力响应问题的影响. 基于Biot波动方程、达西定律和Lord-Shulman广义热弹性理论, 并引入了考虑黏弹性松弛时间因子的Kelvin-Voigt黏弹性模型研究地基上表面受热/力源作用时, 孔隙率和渗透系数变化对均质各向同性饱和多孔黏弹性地基中所考虑的各无量纲量的影响. 根据不同的边界条件采用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度的解析表达式, 结合算例分析了不同变量对各物理量的影响. 正则模态法是一种加权残差法, 可不经正、反积分变换将方程快速解耦并消除数值反变换的局限性. 结果表明: 无论何种载荷作用时, 载荷频率变化对所有考虑的物理量均有明显的影响; 孔隙率和渗透系数均对无量纲超孔隙水压力有明显的影响, 当仅考虑热载荷作用时, 孔隙率和渗透系数变化对无量纲温度均无影响. 正则模态法可广泛应用于岩土工程领域, 尤其适用于商业建筑、高速铁路和公路能源基础的热、力学特性研究中. 该研究结果可为工程施工奠定一定的理论基础, 具有一定的指导性意义.      

Chebyshev collocation method and multidomain decomposition for the incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations in primitive variables have been solved by a pseudospectral Chebyshev method using a semi-implicit fractional step scheme. The latter has been adapted to the particular features of spectral collocation methods to develop the monodomain algorithm. In particular, pressure and velocity collocated on the same nodes are sought in a polynomial space of the same order; the cascade of scalar elliptic problems arising after the spatial collocation is solved using finite difference preconditioning. With the present procedure spurious pressure modes do not pollute the pressure field. As a natural development of the present work a multidomain extent was devised and tested. The original domain is divided into a union of patching sub-rectangles. Each scalar problem obtained after spatial collocation is solved by iterating by subdomains. For steady problems a C¹ solution is recovered at the interfaces upon convergence, ensuring a spectrally accurate solution. A number of test cases have been solved to validate the algorithm in both its single-block and multidomain configurations. The preliminary results achieved indicate that collocation methods in multidomain configurations might become a viable alternative to the spectral element technique for accurate flow prediction.     

Transient thermopiezoelectric response of a one-dimensional functionally graded piezoelectric medium to a moving heat source

The thermopiezoelectricity problem of a one-dimensional (1-D), finite length, functionally graded medium excited by a moving heat source is investigated in this paper. The Lord and Shulman theory of generalized coupled thermoelasticity is employed to account for both the finite speed of thermal waves and coupling of temperature field with displacement and electric fields. Except thermal relaxation time and specific heat, which are taken to be constant for simplicity, all other properties are assumed to vary exponentially along the length through an arbitrary non-homogeneity index. Laplace transform has been used to eliminate the time effect, and three coupled fields, namely, displacement, temperature, and electric fields are obtained analytically in the Laplace domain. The solutions are then inverted to time domain using a numerical Laplace inversion method. Numerical examples are displayed to illustrate the effects of non-homogeneity index, length and thermal relaxation time on the results. When the medium is homogeneous, the results of the current paper are reduced to exactly the same results available in the literature.     

Quasi-similar solutions for blast waves with internal heat transfer effects

The point explosion problem with internal heat transfer effects taken into account is analysed. The classical inviscid solutions to this problem yield a non-physical phenomenon of infinite temperature and zero density at the center of explosion for all times. With heat transfer fluxes considered, the solution near the center of symmetry is improved and finite values for temperature are obtained. The non-self-similar solution of the problem is based on the quasi-similar approximate technique which reduces the non-linear partial differential equations governing the problem to ordinary differential ones. However, this formulation yields a two-point boundary-value problem. To facilitate the integration, the flow field is first divided into two regions: an outer inviscid region and an inner region where dissipation effects are manifested. This results in two sets of ordinary differential equations expressing the conservation equations in the inner and outer regions which are then solved and matched together to yield the composite solution. Secondly, the problem is then transformed into an initial-value one. Using the results of the composite solution, the governing equations can be integrated directly from the center until the shock front. The structure of the non-self-similar flow fields with internal heat transfer effects is then fully determined for specified values of the heat transfer parameters.     

Homotopy Coiflets wavelet solution of electrohydrodynamic flows in a circular cylindrical conduit

In previous studies, the nonlinear problem of electrohydrodynamic(EHD)ion drag flows in a circular cylindrical conduit has been studied by several authors. However, those studies seldom involve the computation for large physical parameters such as the electrical Hartmann number and the magnitude parameter for the strength of the nonlinearity due to the existence of strong nonlinearity in these extreme cases. To overcome this faultiness, the newly-developed homotopy Coiflets wavelet method is extended to solve this EHD flow problem with strong nonlinearity. The validity and reliability of the proposed technique are verified. Particularly, the highly accurate homotopy-wavelet solution is obtained for extreme large parameters, which seems to be overlooked before.Discussion about the effects of related physical parameters on the axial velocity field is presented.     

Stress and pressure fields around two wellbores in a poroelastic medium

The problem of two circular wellbores of different size in a poroelastic medium is considered in the present work. The constitutive behaviour of the poroelastic medium is assumed to comply with the classical Biot model for isotropic porous materials infiltrated by compressible fluid. The wellbores are assumed infinitely long and the fluid flow is taken stationary, thus making it possible to perform a plane strain analysis. Owing to the geometrical layout of the system, bipolar cylindrical coordinates have been adopted. Three different sets of BCs on the pressure field and on the fluid flux have been considered, founding the corresponding forms of the pressure field. Based on Helmholtz representation, a displacement potential has been introduced, and the corresponding stress field in the poroelastic medium has been assessed. However, such a solution does not satisfy the BCs at the edges of the wells. Then, an auxiliary stress function, which allows accomplishing the BCs, is introduced, leading to the complete solution of the problem. The cases of two coaxial wellbores (eccentric annulus), a single hole bored in a poroelastic half plane and two intersecting holes have been considered also. The proposed approach allows evaluating the pore pressure and the stress and strain fields in the system varying the amplitude of the wells and the physical parameters of the porous material. In particular, the evaluation of the peak values of the stress components around the circular boreholes plays a key role in a variety of engineering contexts, with particular reference to the stability analysis of wellbores and tunnels and failure of vascular vessels in biological tissues.