非线性粘弹理论中的单积分型本构关系

本文综述了非线性粘弹理论中的单积分本构表达,评述了多种有代表性的单积分型非线性粘弹理论,对几种本构方程加以分析比较,以揭示它们的内涵,明了其非线性表述原理。      

Burgers方程的小波精细积分算法

求解偏微分方程的常用方法包括有限差分法、有限元法等。近年来,小波分析在偏微分方程数值求解中的应用已引起很多学者的关注,例如采用Daubechies小波或shannon小波构造的小波配置方法已经取得较好的结果。钟万勰院士提出的偏微分方程的子域精细积分方法是一种半解析方法,方法简单,精度高。将小波方法和精细积分方法相结合应用于偏微分方程的数值求解中将有利于提高算法的精度和稳定性,为此本文以Burgers方程为例,提出了一种求解一维非线性抛物型偏微分方程的小波精积分方法。该方法用拟小波配点法对空间域进行离散,建立起对时间的常微分方程组,然后采用精细时程积分方法对该方程组求解。数值计算结果表明,该方法同其它方法相比,具有计算格式简单,数值稳定性和精度较高的优点。     

非粘滞阻尼系统时程响应分析的精细积分方法

考虑一个具有非粘滞阻尼特性的多自由度系统响应的时程分析问题.该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式体现为阻尼力等于质点速度与某一核函数的卷积.在利用状态空间方法将系统运动方程转换成一阶的状态方程的基础上,采用精细积分方法对状态方程进行数值求解,得到一种求解该阻尼系统时程响应的精确、高效的计算方法.通过两个数值算例表明,采用该方法得到几乎精确的数值计算结果,而且计算效率有成数量级的提高.     

率相关本构方程积分新算法

提出一种积分率相关本构方程的隐式积分新算法,引入0～1范围内的缩放因子λ对本构方程进行间接求解,可以避免直接求解等效塑性应变或等效塑性应变率时,由于其数值过大或过小而造成的收敛困难或收敛失败,实现对率相关本构方程的快速准确求解.以B-P统一本构方程及双曲正弦本构方程为例,验证了算法的可行性.结果表明,新算法对于准静态变形条件下的无硬化本构方程也可以得出准确的解.     

精细积分的非线性动力学积分方程及其解法

给出了非线性动力学积分方程的表达式,针对该方程提出了一个显式预测-校正的单步四阶精度的精细积分算法,适用于多自由度、强非线性,非保守系统.算例表明该方法精度高、计算量较少.     

结构非线性动力方程的精细积分算法

基于线性方程精细积分的思路,对具有惯性、阻尼、刚度非线性的动力方程及参变非线性动力方程提出了一种较高精度线性化精细积分迭代计算算法,算例表明该算法可用较大的步长取得满意的计算精度,并可在较大的线性化区间获得较高的计算精度。     

基于内变量理论的一种广义粘弹性本构方程

以Ｂｉｏｔ理论为基础，将弹性应变与非弹性尖变均取作状态为量，引入耗散势，建立了一种率相关非弹性本构关系及演变方程，进而导出一种广义的粘弹性本构方程。     

变温粘弹性蠕变型本构方程

文中给出了变温蠕变曲线的概念，给出了由一组恒温蠕变曲线确定沿某一温度历史的变温蠕变曲线的方法，给出了终态温度等效蠕变曲线的概念，最后给出了粘弹性材料在变温下的蠕变型本构方程的一般形式。文中证明了热流变简单材料理论是本文理论的特例，揭示了热流变简单材料理论数学方法的物理意义。     

有限元分析中硬化弹塑性本构方程积分算法

基于增量段上的应变关于塑性拉氏乘子的变化率为常矢量的假定,导出一种精确有效的、对有限元分析中线性混合硬化弹塑性本构方程的积分算法。以精确算例和等误差分布图的形式检验了本文方法以及其他两种常用方法的精度。     

瞬态热传导方程的子结构精细积分方法

对线性定常结构动力系统提出的精细积分方法,在数值精度等方面表现出极大优越性,但是当矩阵尺度很大时在数值计算与存储中将产生困难,对此,本文对瞬态热传导方程,根据结构的概念,将结构分为若干个子结构,对各子结构分别进行指数矩阵运算并通过了结构间界面的物理量相联系,从而提高精细积分方法的计算效率。     

矩形空腔内Stokes流的状态空间有限元法

基于Hellinger-Reissner二类变分原理,从平面Stokes流问题的平衡方程、连续性要求和边界条件出发,得到相应的Hamilton函数,建立Hamilton正则方程后,采用分离变量法对场变量进行离散求解:在x方向采用有限元插值,在y方向采用状态空间法给出控制坐标方向的解析解。计算过程中的指数矩阵均采用精细积分法求解,使得本文算法具有高效率、高精度、对步长不敏感的优点。通过对侧边自由液面边界条件的单板驱动矩形空腔Stokes流问题的求解,得到与文献相同的结果,从而验证了本文方法的有效性。本文旨在将弹性力学状态空间有限元法的思想引入到低雷诺数流体力学中,为Hamilton体系下研究复杂边界Stokes流问题提供新的途径。     

The finite element method of viscoelastic large deformation plane problem with Kirchhoff stress tensors-Green strain tensors constitutive relation

Based upon the updated Lagrangian approach, the principle of virtual work denoted by the updated Kirchhoff stress increment tensors and the updated Green strain increment tensors and the integral constitutive relation expressed by Kirchhoff stress tensors and Green strain tensors are used and the viscoelastic large deformation incremental variational equation is derived. By means of the 8-nodes isoparametric finite element the program of two-dimensional problem is written. Good agreement is found among the results obtained from this paper and other literatures.     

Three‐dimensional simulation of planar contraction viscoelastic flow by penalty finite element method

The planar contraction flow is a benchmark problem for the numerical investigation of viscoelastic flow. The mathematical model of three‐dimensional viscoelastic fluids flow is established and the numerical simulation of its planar contraction flow is conducted by using the penalty finite element method with a differential Phan‐Thien–Tanner constitutive model. The discrete elastic viscous split stress formulation in cooperating with the inconsistent streamline upwind scheme is employed to improve the computation stability. The distributions of velocity and stress obtained by simulation are compared with that of Quinzani's experimental results detected by laser–doppler velocimetry and flow‐induced birefringence technologies. It shows that the numerical results agree well with the experimental results. The numerical methods proposed in the study can be well used to predict complex flow patterns of viscoelastic fluids. Copyright © 2009 John Wiley & Sons, Ltd.     

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

基于实体造型的三维有限元网格自动剖分

本文叙述了用八叉树法对三维实体进行有限元网格自动剖分算法的设计和实现。本文的算法对传统的八叉树法进行了一系列的改进，包括新型数据结构的建立，对八叉树结点属性判断算法的改进，边界结点体处理方法的改进等。使得本文的算法既保持了传统算法中自动化程度高、层次结构分明、适于再进行网格加密等特点，又克服了其所需存贮空间大、执行速度慢、边界处理复杂、边界单元形状质量不好等缺点，使算法的实现取得了令人满意的结果     

An enriched finite element algorithm for the implicit simulation of the Stefan problem

An implicit enriched finite element algorithm is proposed to simulate heat transfer involving isothermal phase changes. This technique is based on a mixed variational formulation discretized by means of an enriched finite element approximation of the enthalpy in space. The interface is implicitly described without coupling with an interface-capturing technique. The time integration is carried out with an implicit (backward) Euler algorithm in time. Two examples in 1D and 2D clearly evidence the efficiency of the method developed.     

子域精细积分及偏微分方程数值解

对于偏微分方程半解析法的方程，精细时程积分虽然能求出高度准确的解，但往往面临矩阵尺度太大的困难；另一方面差分法虽然有带宽小的优点，但有稳定性及精度方面的问题．本文提出子域精细积分法，既可利用精细积分的数值优点，又有带宽小的好处．数值例题表明了子域精细积分法的效能．     

黏弹夹芯板的有限元建模及实验研究

从有限元分析和数值模拟及实验验证的角度研究了黏弹夹芯板的频率依赖振动特性。夹芯板中间层为黏弹性材料,其刚度和阻尼的频率依赖性行为直接影响系统的模态频率和阻尼,并导致振动模式求解的复杂化。采用三阶七参数Biot模型描述黏弹性材料频率相关的黏弹性行为。开发了三层四节点28自由度的夹芯板单元,基于经典板理论和哈密顿原理建立了黏弹夹芯板的有限元动力学方程。通过引入辅助耗散坐标,将Biot模型和黏弹夹芯板的有限元动力学模型结合起来,并将其转化为常规二阶线性系统形式,极大简化了求解非线性振动特性的过程。对一边固定、另三边自由的黏弹夹芯板进行了前三阶固有频率和损耗因子的预测,并与实验结果对比。数值模拟结果和实验结果吻合良好,说明所提有限元方法是正确有效的。     

Criteria for finite element algorithm of generalized heat conduction equation

To eliminate oscillation and overbounding of finite element solutions of classical heat conduction equation, the author and Xiao have put forward two new concepts of monotonies and have derived and proved several criteria. This idea is borrowed here to deal with generalized heat conduction equation and finite element criteria for eliminating oscillation and overbounding are also presented. Some new and useful conclusions are drawn.     

判别二维有限元网格图的简便算法

根据图论中欧拉公式推得二维有限元网格图中单元数、节点数、边数及边界数之间的关系式,获得判别网格图中是否存在节点号重复和单元遗漏的简便算法。