黏弹性地基上梁的力学行为分析

 考虑地基土(或基础)的黏弹性性态，基于黏弹性力学的基本原 理，采用Fourier级数，建立了梁的力学分析模型，给出了Maxwell与 Kelvin地基(或基础)上梁的变形和内力表达式. 文末的工程实例表明， 所建立的力学模型正确，可用于对各种工程中的基础梁分析.     

多孔有限大弹性薄板弯曲应力集中问题

应用弹性力学的复变函数理论,采用多保角变换的方法,推出了含有任意多孔有限大弹性薄板弯曲的多复变量应力函数的表达式.在内边界上进行复Fourier级数展开,在外边界采用配点法来确定应力函数的未知系数,从而计算有限大弹性薄板的应力场.本文以外边界为矩形,内边界为任意多椭圆孔的有限薄板为例,编制了相应的计算程序,进行了算例分析.结果表明本方法对处理多孔有限大弹性平面问题是简单且行之有效的.     

黏弹性地基上功能梯度材料板的振动分析

为研究黏弹性地基上功能梯度材料板的自由和强迫振动特性,基于Reddy高阶剪切变形理论以及由Shen导得的广义Karman型方程,用双重Fourier级数法推导了三参数黏弹性地基上四边简支功能梯度材料板自由振动和动力响应的解析解,计算了各模态自振频率和半波冲击载荷作用下的动力响应,讨论了材料组分指数、黏弹性地基参数、边厚比等因素对自由振动和动力响应的影响.结果表明,黏弹性地基的剪切和压缩刚度显著提升了功能梯度材料板的振动频率,减小了动力响应;另外,地基的黏性对振动频率和动力响应也有一定的影响.     

不确定结构时域响应分析的多项式维数分解法

针对具有不确定参数结构,提出时域不确定性传播和量化的多项式维数分解法,确定了结构响应统计量的演变过程.首先,采用参数概率模型来描述结构参数的不确定性,建立结构动力学方程,将结构响应表达为不确定参数的函数;进一步,将所关心的结构响应采用成员函数进行维数分解,并利用正交多项式基底对成员函数进行Fourier展开;最后,应用降维积分方法进行展开系数的求解,给出了响应均值和标准差的计算表达式.在数值算例中,将本文方法与蒙特卡洛方法进行对比,结果表明所建立方法具有较高的求解效率和计算精度.     

黏弹性管道水柱分离弥合水锤有限体积法模型

基于有限体积法二阶Godunov求解格式对黏弹性输水管道中水柱分离弥合现象进行建模和模拟研究. 在传统的弹性管道模型基础上考虑管道黏弹性效应的影响. 在瞬变流控制方程中引入管道黏弹性项和动态摩阻项, 采用有限体积法进行求解, 考虑压力修正系数来模拟自由气体对计算单元的影响, 同时为避免数值模拟结果产生虚假震荡引入斜率限制器MINMOD函数; 通过虚拟单元法进行边界构建, 实现了计算区域的统一计算. 将所建模型计算结果与已有模型结果、试验结果进行对比, 并对影响模型的各参数进行敏感性分析. 结果表明, 本文模型能够准确模拟出纯水锤、水柱分离弥合水锤两种情况下的瞬态压力变化, 均能与试验数据高度吻合; 与传统的特征线方法相比, 当库朗数C_r小于1时, 有限体积法二阶Godunov格式计算结果更准确、稳定; 在压力波动的衰减过程中, 黏弹性效应相比于管道摩阻起主导作用; 与弹性管道模型相比, 考虑管道黏弹性效应后可显著提高模拟结果的准确性, 尤其是压力波峰值的相对误差明显降低.      

弹性地基上梁的力学行为分析

考虑地基土（或基础）的弹性性态，基于黏弹性力学的基本原理，采用Fourier级数，建立了梁的力学分析模型，给出了Maxwell与Kelvin地基（或基础）上梁的变形和内力表达式。末的工程实例我表明，所建立的力学模型正确，可用于对各种工程中的基础梁分析。     

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型,介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展.通过线性稳定性理论,分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, DarcyBrinkman-Oldroyd以及Darcy-Brinkman-Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响.利用弱非线性分析方法,揭示失稳临界点附近热对流流动的分叉情况,以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式.采用数值模拟方法,研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的,而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞,最后发展为混沌状态.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

尼龙6/蒙脱土复合材料蠕变行为研究

针对经典黏弹性模型不能很好分析黏弹性材料的蠕变行为问题,运用分数阶导数的类标准线性体模型与Prony级数模型研究了黏弹性材料尼龙6/蒙脱土复合材料的蠕变行为.采用原位聚合法制备了尼龙6/蒙脱土复合材料,在室温环境下对其进行蠕变实验.然后运用分数阶导数的类标准线性体模型和Prony级数模型对复合材料的蠕变实验数据进行分析...     

Parameter determination and response analysis of viscoelastic material

The parameter determination of viscoelastic material is a multi-variable, multi-aim nonlinear optimization problem, which made the optimization process very complicated. In this paper a hybrid optimal algorithm was proposed to determine the viscoelastic parameters in the constitutive relation according to the experimentally obtained mechanical properties. This algorithm merges the Broydon–Fletcher–Goldfarb–Shanno search into a genetic algorithm framework as a basic operator in order to enhance the local search capability. The proposed hybrid algorithm not only can reduce the iterative times greatly but can abolish the limitation of initial parameter values. Nonlinear material characteristic curve-fitting was carried out using the proposed algorithm and other existing approaches. And the comparison results show this algorithm is accurate and effective. The numerical simulation and experimental study of viscoelastic cantilever beam also indicates that the finite element formulation and the calculative viscoelastic model parameters are reliable. The proposed optimization method can be extended to further complex parameter estimation researches.     

The Flow Analysis of Viscoelastic Fluid with Fractional Order Derivative in Horizontal Well

The fractional calculus approach is introduced into the seepage mechanics. A three-dimensional relaxation model of viscoelastic fluid is built. The models based on four boundary conditions of exact solution in Laplace space for some unsteady flows in an infinite reservoir is obtained by using the Laplace transform and Fourier sine and cosine integral transform. The pressure transient behavior of non-Newtonian viscoelastic fluid is studied by using Stehfest method of the numerical Laplace transform inversion and Gauss–Laguerre numerical integral formulae. The viscoelastic fluid is very sensitive to the order of the fractional derivative. The change rules of pressure are discussed when the parameters of the models change. The plots of type pressure curves are given, and the results can be provided to theoretical basis and well-test method for oil field.     

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

State space approach to unsteady magnetohydrodynamics natural convection heat and mass transfer through a porous medium saturated with a viscoelastic fluid

A technique of the state space approach and the inversion of the Laplace transform method are applied to dimensionless equations of an unsteady one-dimensional boundary-layer flow due to heat and mass transfer through a porous medium saturated with a viscoelastic fluid bounded by an infinite vertical plate in the presence of a uniform magnetic field is described. Complete analytical solutions for the temperature, concentration, velocity, and induced magnetic and electric fields are presented. The inversion of the Laplace transforms is carried out by using a numerical approach. The proposed method is used to solve two problems: boundary-layer flow in a viscoelastic fluid near a vertical wall subjected to the initial conditions of a stepwise temperature and concentration and viscoelastic fluid flow between two vertical walls. The solutions are found to be dependent on the governing parameters including the Prandtl number, the Schmidt number, the Grashof number, reaction rate coefficient, viscoelastic parameter, and permeability of the porous medium. Effects of these major parameters on the transport behavior are investigated methodically, and typical results are illustrated to reveal the tendency of the solutions. Representative results are presented for the velocity, temperature, concentration, and induced magnetic and electric field distributions, as well as the local skin-friction coefficient and the local Nusselt and Sherwood numbers.     

IMPROVEMENT OF EXPLICIT ALGORITHMS STABILITY WITH VISCO-ELASTIC ARTIFICIAL BOUNDARY ELEMENTS 1)

黏弹性人工边界单元是目前常用的处理半无限空间波动问题的数值模拟方法,可有效吸收计算区域内产生的外行波动.黏弹性人工边界单元具有与内部介质不同的质量密度、刚度和阻尼,受其影响,对整体模型进行显式时域逐步积分时,在边界区域易发生失稳现象,影响整体系统显式积分的计算效率. 针对该问题目前尚无行之有效的解决方法.本文针对二维黏弹性人工边界单元,建立可代表整体系统典型特征的侧边子系统和角点子系统,利用传递矩阵谱半径分析方法,基于传统中心差分格式,推导得到局部子系统稳定性条件的解析解.在此基础上通过研究解析解中各物理参数对稳定性条件的影响,给出通过增加人工边界单元的质量密度,以改善采用黏弹性人工边界单元时显式算法稳定性的方法.均匀和成层半空间波动问题算例分析表明,将内部单元质量密度设置为人工边界单元质量密度的上限,可以在保证黏弹性人工边界计算精度的前提下,有效改善整体系统显式时域逐步积分的数值稳定性,大幅提高计算效率.     

采用黏弹性人工边界时显式算法稳定性条件

黏弹性人工边界是处理无限域波动问题常用的数值模拟方法。采用显式时域逐步积分算法进行计算时,受黏弹性人工边界的阻尼、刚度等影响,人工边界区的稳定性比内部计算域的更严格,尚无明确、实用的稳定性判别准则,这限制了黏弹性人工边界在显式动力分析中的应用。针对二维黏弹性人工边界,利用基于局部子系统的稳定性分析方法和基于传递矩阵谱半径的稳定性判别准则,给出了可代表整体模型局部特征的不同边界子系统的稳定性条件解析解。通过对比分析不同计算区域的稳定性条件及其影响因素,证明了整体模型的稳定性由角点子系统控制。在此基础上,获得了含黏弹性人工边界的整体模型在显示动力计算中的统一稳定性判别准则和简化实用计算方法。在实际应用中,令积分时间步长满足稳定性条件,即可顺利完成整体模型的动力计算。以上研究可为将黏弹性人工边界应用于显式动力计算时积分时间步长的合理选取提供参考。     

Convection development in a liquid layer with a free boundary

In certain calculations of the critical Rayleigh number for a liquid layer with free boundary which is heated from below, the linearization method has been used and it has been assumed that the temperature perturbations disappear at the undisturbed free boundary.Proper linearization shows that the temperature perturbation is proportional to the free surface perturbation, and the latter is proportional to the normal stress perturbation with the proportionality factor F=²/gh³ (g is the free-fall acceleration, is the kinematic viscosity, h is the liquid layer thickness). In §1 we present a formulation of the problem with account for the parameter F; in §2 we consider the linearized equations and the existence of a stability threshold is proved-a positive eigenvalue-and it is established that with an increase in the parameter F/P (P is the Prandtl number) the value of the critical Rayleigh number Ra_* decreases; §3 presents the results of a numerical calculation of Ra as a function of the parameter F/P.Convection development in a liquid layer with a free surface on which a given temperature is maintained was studied in [1, 2]. The value R_*=1100 found for the critical Rayleigh number agrees well with the experimental value. In the calculations made in [1, 2] the linearization method is used, and it is assumed that the temperature perturbations disappear at the undisturbed free boundary. Strictly speaking, this assumption is not correct.Correct linearization shows that the temperature perturbation is proportional to the perturbation of the free boundary, and the latter is proportional to the normal stress perturbation (see below (2.3)).The problem formulation is presented in §1; §2 deals with the linearized equations and the existence (Theorem 2.1) is demonstrated of a stability threshold—which is a simple positive eigenvalue; §3 presents the results of a numerical calculation of R_* as a function of the parameter =F/P.     

Extrudate swell through an orifice die

The extrudate swell of a viscoelastic fluid through an orifice die is investigated by using a mixed finite element and a streamline integration method (FESIM), using a version of the K-BKZ model. The free surface calculation is based on a local mass conservation scheme and an approximate numerical treatment for the contact point movement of the free surface. The numerical results show a vortex growth and an increasing swelling ratio with the Weissenberg number. Convergence with mesh refinement is demonstrated, even at a high Weissenberg number of O(587), where the swelling ratio reaches a value of about 360%. In addition, it is found that the effective flow channel at the entrance region next to the orifice die is reduced due to the enhanced vortex growth, which may be a source of flow instability.