耦合热弹塑性问题的泛函及其变分原理*

为用有限元法求解耦合热弹塑性问题,建立适当的泛函是有必要的.本文试图在K.N.Rysinko和E.I.Blinov提出的热弹塑性本构方程和热传导方程[1]的基础上,用泛函分析理论导出耦合热弹塑性问题的泛函.     

非线性热传导方程的解和热像法

文中证明了非线性热传导方程只能有连续解,不可能存在间断解。列出了解非线性热传导方程的一些方法。用热像法给出了平面热波对平面分界面的反射和透射,并证明了热像法的一阶近似性。最后,对激光聚爆氘氚小球,应用高原子序数介质的球壳作为挡过热电子的传播问题作了解释。     

金属材料热处理过程中的瞬态温度场与解态相变的数值计算方法

本文首先将Kirchhoff变换推广到导热系数为温度的多项式的非定常非线性热传导问题.并用分析方法确定热传导问题的边界条件.其次提出以孕育期叠加法并引用线性混合法则来模拟金属热处理过程的多相瞬态相变,较为简便地确定相变的开始时间、相变的种类及相变组织的数量.最后利用三维双重边界元法分析工件多种形式的热处理全过程.算例的数值计算结果表明本文方法行之有效.     

地下结构和岩土介质弹塑性耦合分析的摄动半解析法*

本文研究了一个用于物理非线性相互作用分析的有效的数值方法。结构和介质耦合分析的弹塑性问题可用摄动法转化为几个线性问题,然后对相应的线性问题分别用有限条和有限层法分析地下结构和岩土介质以达到简化计算的目的。这种方法用了两次半解析技术——摄动和半解析解函数——将三维非线性耦合问题化为一维的数值问题。此外,本法是半解析法结合解析的摄动法应用于非线性问题的新进展,同时也是近年来发展的摄动数值法的一个分支。     

多孔结构热传导问题的并行多尺度分析

多孔材料在航空航天、汽车、机械领域应用广泛,由于其微结构的多孔性,需要发展多尺度算法用于其性能预测.本文先应用Fourier变换将多孔区域的热传导问题转换为频域空间的复值问题,然后对频域问题做多尺度渐近分析,通过构造边界层证明了频域方程的多尺度截断误差估计.进一步,本文利用孔洞填充思想提出了一套在无孔区域上研究多孔区域的统一的多尺度方法,构造了一套预测多孔材料热性能的新的并行多尺度算法,结合逆积分变换给出了整个算法的误差估计.     

非线性热传导方程的几类精确解

通过引入恰当的试探函数,将非线性热传导方程化为易于求解的常微分方程组并对其求解,进而得到非线性热传导方程的孤波解、奇异行波解、三角函数周期波解等一些不同形式的行波解.     

基于分数阶热传导方程激光加热瞬态温度场研究

基于分数阶Taylor(泰勒)级数展开原理,建立单相延迟一阶分数阶近似方程,获得分数阶热传导方程.针对短脉冲激光加热问题建立分数阶热传导方程组,并运用Laplace(拉普拉斯)变换方法进行求解,给出非Gauss(高斯)时间分布的激光内热源温度场解析解.针对具体算例数值研究温度波传播特性.结果表明热传播速度与分数阶阶次有关,分数阶阶次增加,热传播速度减小,温度变化幅度增加.分数阶方程可以用于描述介于扩散方程和热波方程间的热传输过程,且对热传播机制与分数阶热传导方程中分数阶项的关系做了深入剖析.     

在材料非线性问题中的摄动有限元法

摄动法是解决非线性连续介质力学问题的一种有效方法.这种方法是建立在该问题的线性解析解的基础上的,因此,若得不到一个简单的解析解,应用这种方法去解决一些复杂的非线性问题将遇到困难.有限元法对解非线性问题也是一种十分有用的工具,然而一般来说,它需要相当长的计算时间. 本文介绍摄动有限元法.这种方法吸取上述两种方法的优点,能够解决更复杂的非线性问题,而且也能大量节省计算机的计算时间. 本文讨论了比例加载下的弹塑性力学问题,并提出一个带孔拉板的数值解.     

裂纹尖端场弹塑性分析的加权残数法及塑性应力强度因子的计算

本文以幂强化材料,平面应变情形为例,系统地提出了裂纹尖端场弹塑性分析的加权残数法,并根据此法,得出了裂纹尖端场的解析式弹塑性近似解.在此基础上.对整个裂纹区域,构造了弹塑性解叠加非线性有限元计算塑性应力强度因子的方法,从而为裂纹尖端场和整个裂纹体的分析和计算,提供了一个方法.     

半无界非线性热传导方程的Laguerre拟谱方法

以Laguerre-Gauss-Radau节点为配置点,利用拟谱方法求数值解,逼近半无界非线性热传导方程非齐次Neumann边界条件的正确解.给出算法格式和相应的数值例子,表明所提算法格式的有效性和高精度.这里所用方法也可用于求解其他非线性问题.     

A local iteration scheme for nonlinear two-dimensional steady-state heat conduction: a BEM approach

An efficient algorithm is proposed to solve the steady-state nonlinear heat conduction equation using the boundary element method (BEM). Nonlinearity of the heat conduction equation arises from nonlinear boundary conditions and temperature dependence of thermal conductivity. Using Kirchhoff's transformation, the case of temperature dependence of thermal conductivity can be transformed to the nonlinear boundary conditions case. Applying the BEM technique, the resulting matrix equation becomes nonlinear. The nonlinearity, however, only involves the boundary nodes that have nonlinearboundary conditions. The proposed local iterative scheme reduces the entire BEM matrix equation to a smaller matrix equation whose rank is the same as the number of boundary nodes with nonlinear boundary conditions. The Newton-Raphson iteration scheme is used to solve the reduced nonlinear matrix equation. The local iterative scheme is first applied to two one-dimensional problems (analytical solutions are possible) with different nonlinear boundary conditions. It is then applied to a two-region problem. Finally, the local iterative scheme is applied to two cavity problems in which radiation plays a role in the heat transfer.     

Nonlinear thermoelastic analysis of layered structures

The nonlinear thermoelastic behavior of orthotropic layered slabs and cylinders including radiation boundaries, temperature-dependent material properties, and stress-dependent layer interface conditions is investigated. A one-dimensional finite element formulation employing quadratic layer and linear interface elements is used to perform the analyses. The transient heat conduction portion of the program is temporally discretized via an implicit linear time interpolation algorithm which includes Crank-Nicolson, Galerkin, and Euler backward differencing. The nonlinear heat conduction equations are iteratively evaluated using a modified Newton-Raphson scheme. Direct iteration between heat conduction and stress analysis is employed when stress-dependent interface thermal resistance coefficients are utilized. Verification problems are presented to demonstrate the accuracy of the finite element code.     

四段式有限元法分析热耦合的流体-固体相互作用问题

给出了一种流(体)-热-结构综合的分析方法,固体中的热传导耦合了粘性流体中的热对流,因而在固体中产生热应力.应用四段式有限元法和流线逆风Petrov-Galerkin法分析热粘性流动,应用Galerkin法分析固体中的热传导和热应力.应用二阶半隐式Crank-Nicolson格式对时间积分,提高了非线性方程线性化后的计算效率.为了简化所有有限元公式,采用3节点的三角形单元,对所有的变量:流体的速度分量、压力、温度和固体的位移,使用同阶次的插值函数.这样做的主要优点是,使流体-固体介面处的热传导连接成一体.数个测试问题的结果表明,这种有限元法是有效的,且能加深对流(体)-热-结构相互作用现象的理解.     

Numerical analysis of quenching – Heat conduction in metallic materials

Metallic materials present a complex behavior during heat treatment processes. In a certain temperature range, change of temperature induces a phase transformation of metallic structure, which alters physical properties of the material. Indeed, measurements of specific heat and conductivity show strong temperature-dependence during processes such as quenching of steel. Several mathematical models, as solid mixtures and thermal–mechanical coupling, for problems of heat conduction in metallic materials, have been proposed. In this work, we take a simpler approach without thermal–mechanical coupling of deformation, by considering the nonlinear temperature-dependence of thermal parameters as the sole effect due to those complex behaviors. The above discussion of phase transformation of metallic materials serves only as a motivation for the strong temperature-dependence as material properties. In general, thermal properties of materials do depend on the temperature, and the present formulation of heat conduction problem may be served as a mathematical model when the temperature-dependence of material parameters becomes important. For this mathematical model we present the error estimate using the finite element method for the continuous-time case.     

Nonlinear heat equation for nonhomogeneous anisotropic materials: A dual‐reciprocity boundary element solution

A dual‐reciprocity boundary element method is presented for the numerical solution of initial‐boundary value problems governed by a nonlinear partial differential equation for heat conduction in nonhomogeneous anisotropic materials. To assess the validity and accuracy of the method, some specific problems are solved. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Parametric study of simultaneous transient conduction and radiation in a two-dimensional participating medium

The lattice Boltzmann method (LBM) has been used to solve the energy equation of a transient conduction–radiation heat transfer problem in a two dimensional Cartesian enclosure filled with an emitting, absorbing and scattering media. The control volume finite element method (CVFEM) is used to obtain the radiative information. Based on this new nonlinear hybrid algorithm, the effects of various influencing parameters on the transient thermal response such as the scattering albedo, the conduction–radiation parameter, and the wall emissivity are studied on the distributions of temperature, radiative heat fluxes. Numerical results are presented as compared with other published works and they are found to be in satisfactory agreement. The advantages of the proposed numerical approach include, among others, simple implementation on a computer, accurate CPU time and capability of stable simulation. Therfore, the method can capture fundamental behaviours in thermal flows of engineering interest, in addition it will have computaional advantages when the geometry is more complex.     

Numerical and experimental study of mitigation of welding distortion

Welding stresses and deformations are closely related phenomena. During the heating and cooling cycles thermal strains may occur in the weld and adjacent area. The strains produced during the heating stage of welding are always accompanied by plastic deformation of the metal. The stresses resulting from these strains combine and react to produce internal forces that cause a variety of welding distortions. Welding deformation needs to be minimized and also the designer should know before hand the extent of deformation so that it can be accounted for in the design as well as in the construction stages.In this paper, heat sinking as a method of distortion mitigation has been studied. Heat sinking has been affected by circulating water through channel clamped at the bottom surface of the plates undergoing welding. The pseudolinear equivalent constant rigidity concept has been used in this investigation for thermo-mechanical analysis of plates undergoing welding with simultaneous heat sinking. The initial nonlinear problem with varying modulus dependent on temperature is transformed into a pseudolinear equivalent system of constant rigidity that is solved by linear analysis.The numerical results compared very well with those of the experimental ones. The proposed concept is found to be computationally more efficient and simpler to model compared to FEM for solving similar thermo-elasto-plastic nonlinear problems. The procedure presented in this work and the results thus obtained, holds a great promise for determining the heat sinking parameters for effectively controlling welding distortion.     

Use of geometry in finite element thermal radiation combined with ray tracing

In heat transfer for space applications, the exchanges of energy by radiation play a significant role. In this paper, we present a method which combines the geometrical definition of the model with a finite element mesh. The geometrical representation is advantageous for the radiative component of the thermal problem while the finite element mesh is more adapted to the conductive part. Our method naturally combines these two representations of the model. The geometrical primitives are decomposed into cells. The finite element mesh is then projected onto these cells. This results in a ray tracing acceleration technique. Moreover, the ray tracing can be performed on the exact geometry, which is necessary if specular reflectors are present in the model. We explain how the geometrical method can be used with a finite element formulation in order to solve thermal situation including conduction and radiation. We illustrate the method with the model of a satellite.     

三维热传导型半导体问题的交替方向有限元方法

1　引　　言三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述[1 ,2 ] ,记 Ω为 Ω=[0 ,1 ] 3的边界 ,三维问题-Δψ =α( p -e+ N( x) ) ,　　 ( x,t)∈Ω× [0 ,T] ,( 1 .1 ) e t= . ( De( x) e-μe( x) e ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .2 ) p t= . ( Dp( x) p +μp( x) p ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .3 )ρ( x) T t-ΔT =[( Dp( x) p +μp( x) p ψ) -( De( x) e-μe( x) e ψ) ] . ψ,　　　　　　 ( x,t)∈Ω× ( 0 ,T] . ( 1 .4 )ψ( x,t) =e( x,t) =p( …     

Hybrid differential transform and finite difference method to solve the nonlinear heat conduction problem

In this paper, the differential transform is employed to discuss the behaviors of nonlinear heat conduction problem. A hybrid method of differential transform and finite difference approach is proposed to solve the transient responses of a nonlinear heat conduction problem. Different parameters of the equation and boundary conditions are considered to verify the feasibility of the proposed method to such problems. Simulation results are illustrated and discussed in comparison with the linear case. The results show that the hybrid method can achieve good results for such problems.