关于薄膜大挠度方程的极限解及其特性

关于薄膜大挠度方程的极限解及其特性袁驷（清华大学土木系）ＯＮＴＨＥＬＩＭＩＴＳＯＬＵＴＩＯＮＯＦＴＨＥＬＡＲＧＥＤＥＦＬＥＣＴＩＯＮＯＦＭＥＭＢＲＡＮＥＳ￥ＹｕａｎＳｉ（ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ．Ｂｅｉｊｉｎｇ）Ａｂｓｔｒａｃｔ：Ｔｈｉｓ...     

线性元时程积分按最大模自适应步长公式的证明

该文对运动方程时程积分解法曾提出一种线性有限元自适应步长求解的算法,给出了按最大模自适应步长的计算公式,其中含有步长h的5/2阶的分数阶次。该文对该公式给出数学证明,并通过单自由度和多自由度的数值算例验证了其5/2分数阶次是最优的。     

一维变分不等式问题的自适应有限元分析新探

结构工程中的弹性薄膜接触和杆件弹塑性扭转等问题是典型的变分不等式问题,对其高效精确求解,特别是满足给定精度要求下的自适应求解,是挑战性课题。该文作者新近成功实现了一维变分不等式问题的自适应有限元分析,该文对此进展作一报道。对于变分不等式的有限元求解,该文提出区域二分法和C检验技术,极大提升了松弛迭代的收敛速度,一般4次~5次线性解即可得到收敛的有限元解答,进而采用作者提出的EEP(单元能量投影)超收敛公式计算超收敛解答,用其检验误差并指导网格细分,逐步得到堪称为数值精确解的解答,亦即得到按照最大模度量逐点满足精度要求的解答。该文给出的数值算例表明所提出的算法具有高效、可靠、精确的优良特性。     

一维Ritz有限元超收敛计算的EEP法简约格式的误差估计

该文对一维问题Ritz有限元后处理超收敛计算的EEP(单元能量投影)法简约格式给出误差估计的数学证明,即对足够光滑问题的(>1)次单元的有限元解答,采用EEP法简约格式计算得到的单元内任一点位移和应力(导数)超收敛解均可以达到的收敛阶,即位移比常规有限元解的收敛阶至少高一阶,而应力则至少高二阶。     

二维自由振动的有限元线法自适应分析新进展

有限元线法(FEMOL)是一种优良的半解析、半离散方法,可将其比拟为广义一维问题,进而将一维有限元中单元能量投影(EEP)法及相应的自适应求解技术引入,使FEMOL由半解析方法变为完全解析、数值精确的方法。在对二维线性问题成功地实现了自适应FEMOL分析的基础上,该文进一步报道FEMOL自适应方法在二维自由振动问题中的成功应用和最新进展。该文简要介绍了FEMOL自适应分析二维振动问题的求解策略和技术,整套方法思路清晰、算法严谨、高效可靠,可以得到满足精度要求的自振频率和按最大模度量满足用户事先给定误差限的振型,均为数值精确解。该文给出的数值算例表明所提出的算法具有高效、稳定、通用、可靠的优良特性。     

具有最佳超收敛阶的EEP法计算格式:I算法公式

对一维C0问题的高次有限元后处理中超收敛计算的EEP(单元能量投影)法提出改进的最佳超收敛计算格式,即用m次单元对足够光滑问题的有限元解答,采用该格式计算的任意一点的位移和应力都可以达到h2m阶的最佳超收敛结果.整个工作分为3个部分,分别给出算法公式、数值算例和数学证明.该文是系列工作的第一部分,针对高次单元提出了凝聚形函数的概念,并证明了相关的逼近定理和等价定理,在此基础上给出了具体的算法公式.     

有限元线法求解非线性模型问题——Ⅱ.扭转杆的最优截面

作为有限元线法(FEMOL)求解非线性问题的系列工作之二,本文将该法应用于形状优化问题,对扭转杆的截面优化这一模型问题作了分析求解。文中首先对双连域截面的扭转问题作了FEMOL推导,然后允许结线的长度改变以描述不同的截面形状,再利用若干变换技巧将形状变量及优化条件引入常微分方程(ODE)体系中,从而将问题转换成标准的非线性ODE问题,并由ODE求解器进行求解。文中算例显示了本法对形状优化问题的求解具有方法简洁、实施方便、效率显著等优点。     

索结构找形分析的精确单元方法

索找形问题是索结构分析中所要解决的首要问题,在给定边界条件下,所施加的预张力和外部荷载通过调节索的 形状来平衡。本文研究索结构初始形状确定的精确有限单元法,对于常见的荷载形式,构造了线性和非线性两大类共5种 单元,适用于一般的索结构找形计算,并且可以给出精确的解答。本文通过将水平方向和竖直方向的平衡方程解藕,得出 了索单元的精确描述格式,并保证了索结构形状的唯一性。文中以索结构内部结点坐标作为基本未知量,以结点平衡方程 为基本方程,通过直接求解单元的平衡微分方程得到单元刚度矩阵的解析表达式。对于由线性单元组成的索结构,其基本 方程为线性,可直接求解;对于含有非线性单元的索结构,其基本方程为非线性,需通过迭代求解,文中构造了相应的 Newton法迭代格式。本文方法所需单元数目少,计算量小,所得到的解答为数值精确解。数值算例表明本法稳定可靠。     

运动方程自适应步长求解的一个新进展——基于EEP超收敛计算的线性有限元法

该文采用最简单的Galerkin型线性单元,对运动方程构建了简捷高效的单步法递推公式;进而基于EEP超收敛计算技术,开发了单元步长自动优化和结点位移精度修正两项关键技术,可在整个时域上得到误差分布均匀且逐点满足给定的误差限的解答——堪称数值解析解。该文给出了单自由度和多自由度的数值算例以验证本法的有效性。     

旋转梁振动分析的ODE求解器法

本文提出用标准的常微分方程（ＯＤＥ）求解器直接求解旋转梁的挥舞振动固有频率和振型，主要作法是通过建立旋转梁的挥舞振动方程综合运用ＯＤＥ变换技巧，将其化为标准的非线性问题，再输入求解器求解，本法简单明了且直接可靠，可由小到大求得各阶频率和振型，并使计算结果几乎可达到所需的任意精度。