ANALYTICAL SOLUTION OF THERMAL CONSOLIDATION OF SATURATED SOILS WITH TWO PARALLEL CYLINDRICAL HEAT SOURCES1）

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

两平行圆柱热源作用下饱和土体热固结解析解

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

衰变热源作用下饱和多孔介质热固结问题的扩展精细积分法

热源作用下饱和多孔介质热固结效应是土木及能源工程领域的一个重要课题.由于问题的复杂性,已有的研究大多将介质假定为均匀各向同性,且将热源假定为恒定强度.实际工程中,天然饱和多孔介质常表现出明显的分层特性,热源强度也存在衰变性,为此本工作采用扩展精细积分法对衰变热源作用下层状饱和多孔介质的热固结问题进行研究.借助于积分变换,将饱和多孔介质热固结问题的偏微分方程转化为变换域内的常微分方程;然后对饱和多孔介质微层元进行合并消元,并结合边界条件,推导出衰变热源作用下层状饱和多孔介质热固结问题在积分变换域内的扩展精细积分解;对所得解答进行相应的数值积分逆变换,可获得所求温度、超静孔压及竖向位移在物理域内的解答.基于上述求解过程,编制相应的计算程序进行数值计算,通过与已有文献对比,验证本文扩展精细积分法在求解层状饱和多孔介质热固结问题中的适应性和正确性;最后通过几组算例,分析热源衰变周期、热源埋深及介质的成层性对热固结效应的影响.结果表明:热源衰变周期对温度和超静孔压的峰值、以及达到峰值的时间均有明显影响,衰变周期越长,二者峰值均越大,且达到峰值所需时间越长;热源埋深对超静孔压及竖向位移变化影响显著,深埋热源作用时热源两侧竖向位移呈对称分布,而浅埋热源两侧则无此现象;饱和多孔介质的分层特性对热固结效应影响明显.     

变边界粘弹性轴对称问题的复变函数法

对边界几何形状、位置随时间变化的变边界结构,给出了用复变函数求解粘弹问题的解析方法。文中用拉普拉斯变换结合平面弹性复变方法,对内外边界变化时粘弹性轴对称问题进行求解。引入两个与时间、空间相关的解析函数,给出了变边界情况下应力、位移以及边界条件与解析函数的关系。当解析函数形式部分确定,则可用边界条件求解其中与时间相关的待定函数。求解待定函数的方程一般情况下为一系列积分方程,特殊情况可求得解析解。对轴对称问题中应力边值问题、位移边值问题以及混合边值问题,分别利用边界条件求得相关系数,从而得到了应力与位移的解析表达。当取Boltzmann粘弹模型时,进行不同边值问题的分析。分析显示,应力、位移的形态与大小均与边界变化过程相关,与固定边界粘弹性问题有较大不同。本文解答可用于粘弹性轴对称问题内外边界任意变化及各种边值问题的力学分析。此外,该法可进一步进行荷载非对称、复杂孔型变边界问题的求解。     

分数阶广义热弹性理论下中空柱热弹性分析

基于Laplace变换及特征值法,推导并给出了分数阶广义热弹性理论下中空柱内表面作用有热冲击情况的解析解,通过Laplace数值逆变换法求解得到了位移场、温度场、应力场的分布规律。结果表明:特征值法能准确给出Laplace域内方程组的解;分数阶参数对温度场和应力场有较大影响,对位移场影响较小。作为广义热弹性理论的一种推广,在处理热传导问题时,通过分数阶广义热弹性理论进行研究更科学、全面。     

半无限长压电杆的瞬态热冲击问题

采用具有一个热松弛时间的L－S广义热弹性理论，研究了一维半无限长压电杆在一端受到热冲击时的边值问题。借助拉普拉斯正，反变换技术，在所考虑时间非常短的情况下，对问题进行了求解，得到了压电杆上的位移，压力及温度分布的近似解析解，发现应力及温度分布中分别存在两个阶跃点，并给出了算例。     

半无限粘弹杆的广义磁热粘弹问题

应用Lord-Shulman(L-S)和Green-Lindsay(G-L)广义热弹性理论,研究了在磁场中受移动热源作用的半无限长均质各向同性粘弹杆的磁热粘弹动态响应,并与经典耦合理论进行了对比.给出了杆的广义磁热粘弹耦合的控制方程,借助拉普拉斯积分变换及其数值反变换对控制方程进行了求解,计算得到了杆内温度、应力和位移的分布规律.研究结果表明:时间、热源移动速度和磁场大小对以上分布规律都有一定的影响.     

考虑热渗效应的有限热源固结近似解1）

通过建立考虑热渗效应和热流固耦合效应的饱和土体固结方程,研究了无限长圆柱热固结问题. 利用Fourier 和Laplace 变换及其逆变换,给出了热固结问题的解析解;然后对空间内无限长圆柱形热源问题进行研究,得到非等温条件下柱体周围饱和土体温度、孔隙水压力的近似解,并总结其规律,分析了热渗系数、固结系数对温度作用下土体固结的影响.     

考虑热渗效应的有限热源固结近似解

通过建立考虑热渗效应和热流固耦合效应的饱和土体固结方程,研究了无限长圆柱热固结问题.利用Fourier和Laplace变换及其逆变换,给出了热固结问题的解析解;然后对空间内无限长圆柱形热源问题进行研究,得到非等温条件下柱体周围饱和土体温度、孔隙水压力的近似解,并总结其规律,分析了热渗系数、固结系数对温度作用下土体固结的影响.     

无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.     

Response of saturated porous media subjected to local thermal loading on the surface of semi-infinite space

Heat source function method is adopted in the present paper to derive elementary solutions of coupled thermo-hydro-mechanical consolidation for saturated porous media under conjunct actions of instantaneous point heat source, instantaneous point fluid source and constant volume force. By using the so-called fictitious heat source method and images method, the solutions of a semi-infinite saturated porous medium subjected to a local heat source with time-varied intensity on its free surface are developed from elementary solutions. The numerical integral methods for calculating the unsteady temperature, pore pressure and displacement fields are given. The thermomechanical response are analyzed for the case of a circular planar heat source. Besides, the thermal consolidation characteristics of a saturated porous medium subjected to a harmonic thermal loading are also given, and the fluctuation processes of the field variables located below the center of heat source are analyzed. The project supported by the National Natural Science Foundation of China (50578008) The English text was polished byYunming Chen.     

THERMAL CONSOLIDATION OF LAYERED POROUS HALF-SPACE TO VARIABLE THERMAL LOADING

An analytical method was derived for the thermal consolidation of layered, saturated porous half-space to variable thermal loading with time. In the coupled governing equations of linear thermoelastic media, the influences of thermo-osmosis effect and thermal filtration effect were introduced. Solutions in Laplace transform space were first obtained and then numerically inverted. The responses of a double-layered porous space subjected to exponential decaying thermal loading were studied. The influences of the differences between the properties of the two layers (e.g., the coefficient of thermal consolidation, elastic modulus) on thermal consolidation were discussed. The studies show that the coupling effects of displacement and stress fields on temperature field can be completely neglected, however, the thermo-osmosis effect has an obvious influence on thermal responses.     

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.     

DYNAMIC RESPONSE OF TWO-DIMENSIONAL GENERALIZED THERMOELASTIC COUPLING PROBLEM SUBJECTED TO A MOVING HEAT SOURCE

Based on the Lord and Shulman generalized thermoelasticity theory with one relaxation time, an isotropic semi-infinite plate subjected to a moving heat source has been studied by employing the finite element method directly in time domain, whose distributions of nora dimensional temperature, displacement and stress are illustrated graphically. The results show that the present method is an effective and exact numerical one for solving the thermoelastic coupling problem and is capable of overcoming the defects of traditional integrated transformation and inverse integrated transformation methods. At the same time, the temperature step of the thermal wave front is obtained exactly in contrast with conventional numerical transformation methods.     

采用包络-准则法优化考虑瞬态传热的薄板结构

刚度和强度是薄板结构的两个主要性能。在瞬态传热过程中,考虑热-力耦合,随时间和空间变化的非均匀温度场在结构中会引起热变形和热应力,温度场随时间变化的规律和空间分布依赖于板的厚度变化,进而影响板的刚度和强度。因此,考虑瞬态传热的薄板优化问题具有更强的非线性,更加难以求解。本文给出一种包络-准则方法处理这类结构优化问题。首先,针对外力荷载,进行一个结构柔顺性的优化设计;以这一设计为基础,通过瞬态热-力耦合分析及优化准则,计算多个时刻的优化设计变量并取其包络,对上述优化结果进行迭代修正,以消除瞬态温度场作用下较高的局部应力。优化算例表明,该方法对于考虑瞬态传热薄板优化问题有效。     

考虑热渗效应的有限长圆柱热固结解

在作者已建立的饱和多孔介质耦合非线性热弹性理论基础上,考虑热渗效应,建立了饱和多孔介质耦合热弹性固结方程,并推导了有限长圆柱热固结问题的解析解,进而以温控三轴试验的试样为例进行了算例分析,同时利用COMSOL软件进行了数值模拟,并将解析结果和数值结果进行对比．结果表明：在不排水条件下,影响试样最终孔压大小的参数是：土的泊松比>弹性模量>水的体膨胀系数,渗透系数对孔压变化影响不大,考虑与不考虑水土压缩性不影响孔压的计算;相对线弹性情况,考虑非线弹性膨胀系数随温度变化时的孔压有所下降,轴向应变变化很小．